View Single Post
Old 2021-11-05, 08:02   #8
axn
 
axn's Avatar
 
Jun 2003

2·2,683 Posts
Default

Does this look right?

Code:
merit(p, n) = log(p) * (p-1) * p^(n-1);
find_count(p, t)=my(m); for(i=2, oo, m=merit(p, i); if(m > t, return(i-1)));
old_find_count(p, B1)=floor(log(B1)/log(p));

B1 = 1000000; last_prime = precprime(B1); target_merit = merit(last_prime, 1);
forprime(p=2, last_prime, c = find_count(p, target_merit); print(p, ":", c, ":", old_find_count(p, B1)); if(c == 1, break()))
Code:
2:25:19
3:15:12
5:10:8
7:8:7
11:6:5
13:6:5
17:5:4
19:5:4
23:4:4
29:4:4
31:4:4
37:4:3
41:4:3
43:4:3
47:3:3
53:3:3
59:3:3
61:3:3
67:3:3
71:3:3
73:3:3
79:3:3
83:3:3
89:3:3
97:3:3
101:3:2
103:3:2
107:3:2
109:3:2
113:3:2
127:3:2
131:3:2
137:3:2
139:3:2
149:2:2
151:2:2
157:2:2
163:2:2
167:2:2
173:2:2
179:2:2
181:2:2
191:2:2
193:2:2
197:2:2
199:2:2
211:2:2
223:2:2
227:2:2
229:2:2
233:2:2
239:2:2
241:2:2
251:2:2
257:2:2
263:2:2
269:2:2
271:2:2
277:2:2
281:2:2
283:2:2
293:2:2
307:2:2
311:2:2
313:2:2
317:2:2
331:2:2
337:2:2
347:2:2
349:2:2
353:2:2
359:2:2
367:2:2
373:2:2
379:2:2
383:2:2
389:2:2
397:2:2
401:2:2
409:2:2
419:2:2
421:2:2
431:2:2
433:2:2
439:2:2
443:2:2
449:2:2
457:2:2
461:2:2
463:2:2
467:2:2
479:2:2
487:2:2
491:2:2
499:2:2
503:2:2
509:2:2
521:2:2
523:2:2
541:2:2
547:2:2
557:2:2
563:2:2
569:2:2
571:2:2
577:2:2
587:2:2
593:2:2
599:2:2
601:2:2
607:2:2
613:2:2
617:2:2
619:2:2
631:2:2
641:2:2
643:2:2
647:2:2
653:2:2
659:2:2
661:2:2
673:2:2
677:2:2
683:2:2
691:2:2
701:2:2
709:2:2
719:2:2
727:2:2
733:2:2
739:2:2
743:2:2
751:2:2
757:2:2
761:2:2
769:2:2
773:2:2
787:2:2
797:2:2
809:2:2
811:2:2
821:2:2
823:2:2
827:2:2
829:2:2
839:2:2
853:2:2
857:2:2
859:2:2
863:2:2
877:2:2
881:2:2
883:2:2
887:2:2
907:2:2
911:2:2
919:2:2
929:2:2
937:2:2
941:2:2
947:2:2
953:2:2
967:2:2
971:2:2
977:2:2
983:2:2
991:2:2
997:2:2
1009:2:1
1013:2:1
1019:2:1
1021:2:1
1031:2:1
1033:2:1
1039:2:1
1049:2:1
1051:2:1
1061:2:1
1063:2:1
1069:2:1
1087:2:1
1091:2:1
1093:2:1
1097:2:1
1103:2:1
1109:2:1
1117:2:1
1123:2:1
1129:2:1
1151:2:1
1153:2:1
1163:2:1
1171:2:1
1181:2:1
1187:2:1
1193:2:1
1201:2:1
1213:2:1
1217:2:1
1223:2:1
1229:2:1
1231:2:1
1237:2:1
1249:2:1
1259:2:1
1277:2:1
1279:2:1
1283:2:1
1289:2:1
1291:2:1
1297:2:1
1301:2:1
1303:2:1
1307:2:1
1319:2:1
1321:2:1
1327:2:1
1361:2:1
1367:2:1
1373:2:1
1381:2:1
1399:1:1
axn is offline   Reply With Quote