1日放置した完全数プログラム（少し修正したりなど）

前回の記事の完全数計算プログラム、少し修正しました．
関数の仮引数にメモリを裂くのがもったいなく感じたのでメルセンヌ素数判定を関数に入れ込んでみました
実際に効果あるのかは知らないけど・・・
あと5スタートにして、2→4→2→4と偶数と3の倍数を回避するようにしました

#include <iostream>
#include <boost/multiprecision/cpp_int.hpp>
#include "MT.h"
#include
#define TEST_MAX 100

namespace mp = boost::multiprecision;

unsigned long long Rand(unsigned long long A, unsigned long long B) { return genrand_int32() % (B - A + 1) + A; }

unsigned long long gcd(unsigned long long a, unsigned long long b) { return (a % b == 0) ? b : gcd(b, a % b); }

unsigned long long FelmatL(unsigned long long a, unsigned long long num) {
unsigned long long result = 1, tmp = a, tnum = num;
while (tnum > 0) {
if (tnum % 2 == 1)result = (result * tmp) % num;
tmp = (tmp * tmp) % num;
tnum >>= 1;
}
return result;
}

bool FermatTest(unsigned long long num) {
if (num == 2 || num == 3)return true;
for (unsigned long long cnt = 0; cnt < TEST_MAX; cnt++) {
unsigned long long a = Rand(2, num - 1);
if (gcd(a, num) != 1)return false;
if (FelmatL(a, num) != a)return false;
}
return true;
}

void PerfectNum() {
mp::cpp_int num=2;
std::ofstream outputfile("Perfect Number.txt");
std::cout << "2 : 6\n" << std::endl;
outputfile << "2 : 6\n" << std::endl;
std::cout << "3 : 28\n" << std::endl;
outputfile << "3 : 28\n" << std::endl;
char f = 0;
for(unsigned long long i = 5; i < ULLONG_MAX; i += (f==1)?2:4) {
num *= (f==1)? 4:16;
if (!FermatTest(i))continue;
mp::cpp_int s = 4;
for (unsigned long long n = 2; n < i; n++) {
s = (s * s - 2) % (num-1);
}
if (s==0) {
mp::cpp_int result = (num / 2) * (num - 1);
std::cout << i << " : " << result << "\n" << std::endl;
outputfile << i << " : " << result << "\n" << std::endl;
}
f = 1 - f;
}
outputfile.close();
}

int main(){
init_genrand( (unsigned)time(NULL));
PerfectNum();
}

---

以下計算結果

2 : 6

3 : 28

5 : 496

7 : 8128

13 : 33550336

17 : 8589869056

19 : 137438691328

107 : 13164036458569648337239753460458722910223472318386943117783728128

127 : 14474011154664524427946373126085988481573677491474835889066354349131199152128

607 : 141053783706712069063207958086063189881486743514715667838838675999954867742652380114104193329037690251561950568709829327164087724366370087116731268159313652487450652439805877296207297446723295166658228846926807786652870188920867879451478364569313922060370695064736073572378695176473055266826253284886383715072974324463835300053138429460296575143368065570759537328128

1279 : 54162526284365847412654465374391316140856490539031695784603920818387206994158534859198999921056719921919057390080263646159280013827605439746262788903057303445505827028395139475207769044924431494861729435113126280837904930462740681717960465867348720992572190569465545299629919823431031092624244463547789635441481391719816441605586788092147886677321398756661624714551726964302217554281784254817319611951659855553573937788923405146222324506715979193757372820860878214322052227584537552897476256179395176624426314480313446935085203657584798247536021172880403783048602873621259313789994900336673941503747224966984028240806042108690077670395259231894666273615212775603535764707952250173858305171028603021234896647851363949928904973292145107505979911456221519899345764984291328

2281 : 99497054337086473442435202604522816989643863571126408511774020575773849326355529178686629498151336416502516645641699516813140394897940636561646545947753232301453603583223268085613647233768081645727669037394385696522820301535888041815559513408036145123870584325525813950487109647770743827362571822870567643040184723115825645590386313377067112638149253171843914780065137373446222406322953569124771480101363180966448099882292453452395428270875732536311539266115116490704940164192417744919250000894727407937229829300578253427884494358459949535231819781361449649779252948099909821642207485514805768288115583409148969875790523961878753124972681179944234641016960011815788847436610192704551637034472552319820336532014561412028820492176940418377074274389149924303484945446105121267538061583299291707972378807395016030765440655601759109370564522647989156121804273012266011783451102230081380401951383582987149578229940818181514046314819313206321375973336785023565443101305633127610230549588655605951332351485641757542611227108073263889434409595976835137412187025349639504404061654653755349162680629290551644153382760681862294677414989047491922795707210920437811136712794483496437355980833463329592838140157803182055197821702739206310971006260383262542900044072533196137796552746439051760940430082375641150129817960183028081010978780902441733680977714813543438752546136375675139915776

3217 : 33570832131986724437010877211080384841138028499879725454996241573482158450444042882048778809437690388449535774260849885573694759906173841157438424730130807047623655942236174850509108537827658590642325482494761473196579074656099918600764404702181660294469121778737965822199901663478093006075022359223201849985636144177185925402078185073015045097727084859464743635537781500284915880244886306461785982956072060013474955617851481680185988557136609224841817877083608951191123174885226416130683197710667392351007374503755403352531476227943590071651702697594241031955529898971218001214641774673134944471562560957179657881556419122102935450299751813340515170956167951095453649485576150660101689160658011770193274226308280507786835049549112576654510119670456745939890194205255175384484489909328967646988163155982471564998196261632751283127879509198074253193409580454562488664383465379885002735506153988851506645137759275553988219425439764732399824712438125054117523837438256744437055019441051006489972341609117978404563794992004873057518455748701444951238377139620494287982489529827233140637014837408856156199515457669607964052126908149265601786094447595560440059050091763547114092255371397425807867554352112542194784815494784276201170845949274674632985210421075531784918358926690395463649721452265405713484388043911634485432358638806645313826206591131266232422007835577345584225720310518698143376736219283021119287617896146885584860065048876315701088796219593640826311622273328035603309475642390804499460156797855361018246696101253922254567240908315385468240931846166962495983407607141601251889544407008815874744654769507268678051757746956891212485456261121386667407711139619071530923355823178662705374393035049022603882479742334799407130280148769298597743778193050348749740786928096033906295910199238181338557856978191860647256209708168229116156300978059197026855726877649767072684960463452763160384093838292277544911857859658328888332628525056

4253 : 18201749040140430273835167914751015293362889506937596103357359437740048438463485836303232258902400875568993809975036246337556768238366233707757573415903278080066115951092522543145952013119010340408099881835975590975835499003463227280543140756802691525621894771543951403656868707748606069885095036335283528581813088651653833089052785834280723280067200484903706704970113521040604064538274154252165373175595605248628998615687866235395077520314640188634353318831397392517095413127961059169783805545093120277974213161792712750189234973622853507624924636687588461310469170436099238232492598937308477793932361698480685664770932715430491007662682405163044765107996164343261202355969912750334459238591875176292279879181327972920090211444842306766331109269166235125509668704478958734853068636529397153443016761158813447481567868262367893981171020212189792514086148090215189947704212407282687964267041783873543080261582448824858431838554183153705837057341487700727731550997896239330477471263230726484877117521254868468174015562068490713644719481693089433418532636835808942200716440819313366362403191355576523173825131875335647759016207232990842875902548737480801898289658229210850221695663421649611729577245208490877897945867171915873804061345544603919338281884397538634259920738193259330950837923823615411304543599290641330796754637229721314074245135865275394292301532835414540035740253567232435673868968242806554916467387945199026147468808031055150281915733287931898181708779241724409277593094799901224641825850183281480975988801751016901958988814849430087926938683408662084176896164178170688871544229415824450055526787148715365214142372692896788902761679810950828268208180845309986581369910322372338718437892376924042425319226250711104694125177110946650435829229609212441927473293781744486775943173110295432226780186422343937384580708855852269417012387614583884866613236875101300845603785945237925447338784922452682337026775686191693583794243538309813874806083439394116832472101884074936138517168612519077402832074575032255369216417037979614393919019507244388546924316787128440250098924928769352393698587810869831127888084696751795845351130410199628821737636610886633120714458623846784948362409208330416268144841565505869827619451490634144640754203543507421045031995754480556197763855213987476437990952173006055503017305481394836591954388794314558786433824434673257753459825605748638411145315395326160048592155490496363095590187403880545069860151237473584639382355987779291886392614452838326884468500853448420623482240841078202705004101172745586672975106848437133377536

9941 : 598885496387336159215158337179296643572869874198415438721517031956294391267467138052833975210822790379645435252765184809790907879059390832431456890763826317589211758237474574018159857925570634144282954314842906400063553663884242875986410271598361502228854072762831378553289935007636618509623086614399709248699419865215841100064058520556986140758791562598688902501114759998860081837041267261910926840834156939147901883683592934604445887790093281004039491643639906988228730211742465708290649815465091229137126198956630386416170359815427838179479070235472610399015035810959640653610781106772716202393062410044423415928785961279533870544231701973513163229825784566847945851733257855753088163409148988722084210091625175105625266415505647159128966267046695766062982729521155903897799050637426474617110612526598459278585486968447249024042786308910850051017106038889061687767151283091964070893823182959298891383459811079403729066740772169275606698939480096480948919175320023826313542796993224153325663400369376016129480939312138469723199723534110595657777687480880412793207072338171511658351904166668048242956453135022989001809464768832188595267188382036606053289907501689359295334416673910623695996312840110689651484585283440480471698222497719355503194168560950281935485655294739301609663777113373659535718967303547951623526956338396985018174631435834656785781957446464286463809460246212247286223440883955123461676820718682834656934175017941968544782841849057554072567204274097863751252547881116514783088210584399664096224612954055771939895263749729051585416405834890522726931499793407634795102451963136189675604110386882667666700202021296939123142592598861407877348129582919695159135624222759532153769870798883632159023890309632619900847829779211195760910227705676555741607808745978228036024945439485997483773117861761846816143038325765449577516576845228775618794694822323427772717649125027086076998811187338852893643099277612863434755567049300968668608075953921793767085814197717905978584997781578440077423937447925873669288761103924157758615611230568564058097818416825907409237840644980565685003936400472616338363054823832113925515656141515345511440781430776107400081182066238156903296613221289208392191798325657897497026838673102803471094612120455368011519554557205053376663671919907929622429762923746679903431626862986005514917878931929165684488842314522696190750751104130758772492742171221506476388227098197462790072396811258045100484996907101399477804665561319641934728133590562936267652166532237603284227934840610439816138264438446708643682085053962432046459648769845201613113808872277523981275235642270948095440191476800941546005369324739646459243063300015084138888650243268999554932706194400103202052002256475547454921546756210434441774249495366559167821730580815611872299157633380987381812207714590674498242069538447774482969637038762378605688864932175145427786356322290946332387532834726569678238370085937344737751948866649747361373588969170028812987868238296955348481127935674306231597404181875425238311701462864830821158701402077328211001889853890086949823322215668362320216048221349456542834993635806413226788683589993095956742107311836481298366437924635728768234856559136863642247989697063646989127406234640433936180524078096916562215885382221221440174447087595474342702340083128927752717896538651110872000256594041548014611801466523195753961690032322786290650070510846917869315870454207312856770297593768844524696425727849441292728482259863151304805682511045957634104047520304681087570488096199426249870968715151352899262711414859028400306247511682097755259210044273287907069828213053162082804915911440204049504420189851777238231126675540195073124882510073151638213199509097027588754146115313488501895618592704144939889978242021363965347241635421115065464033561832225847238709861244250411633809994450569265409560194147735449585298329207794798511236425587495469494467899179662053326961034916517867933367808057191175385336962971973530250601332617559991043382608515546850921036964488304328204714534588455104569762031196764380573914497030120132291308170538348397837511623996297933672121009901941694829665518167940396094144081711936184081318439929556609569648834030333939030167533881530753927743675683048377118347128650493153271095307571548267716602161565111731981963404253792831860675668529002291340501102920781703654993775479142326056554305406934385461270207799502316565734128418269704729154185772068761081919140269479613439718429584336562493991766156958042518706487778614666958547617944366057917035046081816932122930397831885899367394704668767909522089490980411104133713860758251156103716971790705966662212818116148565527822567454545606791772821700331972738940950104057875162524624791534716722495993330597176066518166716680800409152545336552754751419653578617676888721974359792782967881648099394891147826325257661113821525698043547390457514434463473999530885287210000009489867014184718369266887805926383484231605051335786390632303781597551018628117508352980030059724480548533949284681845788829364189309471523884413954194859086869472657365044985433141123019472681576182090016085717562012755665853499189921886618259400375906348395090310900142148715431411955215491914127573351385908146653482483997520485329568230705228757993906938515288896844165579094197758186574987120643678489659211515541853732261407728812733452051953588927628461384781414660582371641376089400661867134510390162683271326617270588318894292606033372701978317622740978994527704380727470767368180702967071072288812051097001504045715988812189767358074825262336006037202238159764792224045818165003646975285929365924426658451465248929517024070447106402584830973235253084098352978893906920352652116340904008935681654751699242318452472319485281464039177405304692021677989652462657737972297496705447616703675014511838739629530737098712639069034012066077188936467469996574123242443286789792602919789894555982048116453916717116130063058668784438394777624864833158156184468926035001575570475329794750896833182478324073496576

11213 : 3959613212817942196968287884854607138249695010687462454364841708851250756886730143010200497312491404781485112128545293038267410557687432146950851043426362765816215434907069232646845830697880715114498237144171635503366409257826346755285153686414093613591743288984395713590371609988997426858436873871662377222014131535012470582972340263924545236916165111578915024271639042423405991498025090677985246615888708814389494097187288049287415979235786510926548495671041495178109433955976169433598463074752794749833261671894416779587164022745726731816558059180128410052463097878416881346421940585665688499556244932031405619167193309877502125598274799692861848222186789839237161415950715076488613323777606157414531516260757170134150987782139353548730742670290744779603220618058769008513945985936788405184681857401721860333491468615517160265447106635779775116189694725965858433646555605352303369006895827491495266728584625923369382805582715543398299534252478300474804990013583056030504517708745604400254795470610884055712298069745897434773516459162650061100638644799141135564026144547699549924290389689515702725711923690519702231947789561422688360442715880502500729718659351717498594995066205790473119554476042449582563590579444721556498539231662047211644081609762161430568424342342242854439576719470899833805336316533137178302986916911543659790598077303875904224315053242517735730246856298337892312316528117634037166026440422177134130244659685698851719028656908543491242029580830020429397597945878491134064323290454637270114479800986458574381828790795066222590999367076594501777434267126916310434035574778547532869917547234512828881393561556657175848151136191220333591579083059191456438448403687051180171942755839609237451068823648275304071278515020701325144649797077281242875887293908713148087425084416129687389477450919839336384319643456353169764604635371056121820998845612577368007974752727989955080372815379334341844407691376688259658161238974567126603085075726800913322215326148218783626818232287223320959841024144071595668287525943886827137995088926789416420742895781629112626832021473618920626418326431458859284960751788990850345628853589391195089799880857621991457010020868383595460091757038998473321366148098536042605377671105942810444355441422008803094913405247372082842185535759207340245260120031482268289664155256889927160530153161598395632054016672558073391690755434186837560699746115376647075646124446313120539536864397378512528982338685436961057988195294019498594823503504893403210126005985201863912422945220587524487839774578151169639124967362823293256513202024803706423568550925173949161614257797912393511359242591869863619478797919472369000047141181208965931111421137288360760624748454875972501545483496405906288261654093770960025270251861037644320567270181065051970090049460704195505843123733155220934945908153070015156209041273961692602971492185076141511781001655030687032702373890902215036807381975613011048299714755167226115128307272101251476288442863767624217431718205364113124570399549481236095990457409373076315965718684534535322865287093771568758027809525357907282500979011293443231165575569445732992370476000920942922538926667095513699704161340592161159746236219773071356763754840639957122750858985185806119754237880099327946072326916962603744275213064634260928453104075010608297854587467237575467613798511046548517920545266623844188959929873419513930739752663887282896737602143712362102641121459357895275428584164006135336277168927794561201868749212479904416098658145097189359149052467301494653771651765986740908064319954228532907359916469374457144363856976480527707515526676064592447766901755037113256345750431361771536700237369409228138098929435715361717105374968249666768874716435458200086730293793319001160279518097064968936679995306918201968957611690366686286176964386475078710733881610365681670783276733231189516158495404203758005340538520187971152745328446019029100477501136956786385655556910835844321697477173648683328687159029719670236620302879815424433038397696239748990986244362183429736060133537688767006408936924928832453938314048346985862079684109578836282050219667623205694961246466059800792378418863583553057111154634775615626877637572769606110362768130814746824813196424289379204437786177275989321295838925014870509289623514139031100726371141651614681759490693712759141745584399038806017019558016549316152380251283298159184556076419304690117554103795728169472045336685854608739179828126567603445118239120699047381406429484442798343983663772265584848128132389264433157905449690280520607189331645267901161622386987161693221619602785888171129338673595852827689093839707616062594377604840077015747288019189892043058477610088173932710132968306471813414899036286602298340723365613904239299983917104952795547250211759660947922790165051457907929340870440808147913806352433854516953162243223960155861517020505099646669991472038855323639652044194991473889260921390120261070269987629039025773816966793417585417122418599005880133843186348612916355868235875600906126702646688542306262184400426360242713095741516009443419776144343751414966901561389451915222341230446882811797103693651699066520715449528551184687804313384261458874649260709607073211038923829085522743537368449086435360105601379194853966106312934064278617268327046200022923910857473355560421960130292197229185669933476428166642771389528948327200871438484659718235903183972267526825777356633446483538626285143021212714440333388352004712871330413326013160289574383811033720860915757855552976096653360492054042681357248970182573621851163508326287628494591608099752379466269094909310662446283468714450702831071516608374697328651174027719486158467818669370866178096050879261993434167401298795985980707851123427900202340473958560450518585690111010553544771878460143855957897625860681757637110872011856062933099161018703113954043675184893355122174208255701575495043843977938871685082192922187720143207515798671187372505602843788344758500881234814579017069183547924250517037035205307303196533903103753462639048730900319270404200704747783372674233117363014608397115419207799450124137593286730107373233455619658143491195747254705698916589686716679282080947654819934175511108443756527314614392729337477755314615118341084990438815138836433800030341206135764693071159103512059488720937401309490936269796380778420747225153504370370779582878038454534715791401735105081399755601073803578039237439792052053490239583704450989419267388307026202215022125780896032731025741953389880626859531073411113775249989905677420885901013463986839439831074089587489357533219454039264412034803623183594322135894548339913390308331649495108148299631153809459746028790153181652286533316810217288835890967080578311471722124410314477370827358400434067739029355888063072530341428376185805388255702691086336


一日放置してみた結果．
次はだから12003桁・・・うーん次の完全数ちゃんと出てくるのかな・・・