USACO算法系列三十四——prime3

    题目：http://www.nocow.cn/index.php/Translate:USACO/prime3

    我没打算做这道题，太恐怖了。看看漂亮的暴力解法

#include<stdio.h>#include<string.h>int a[6][6],ss[100000],n,sum;char ans[201][27];int main(){FILE *fi,*fo;fi=fopen("prime3.in","r");fo=fopen("prime3.out","w");int i,j,k,tmp;fscanf(fi,"%d%d",&n,&a[1][1]);ss[1]=1;for(i=2;i<=99999;i++){if(ss[i]==1)continue;for(j=i*2;j<=99999;j+=i)ss[j]=1;}for(a[2][2]=0;a[2][2]<=9;a[2][2]++)for(a[3][3]=0;a[3][3]<=9;a[3][3]++)for(a[4][4]=0;a[4][4]<=9;a[4][4]++){a[5][5]=n-a[4][4]-a[3][3]-a[2][2]-a[1][1];if(a[5][5]<0||a[5][5]>9)continue;if(a[5][5]%2==0)continue;tmp=a[1][1]*10000+a[2][2]*1000+a[3][3]*100+a[4][4]*10+a[5][5];if(ss[tmp]==1)continue;for(a[5][1]=1;a[5][1]<=9;a[5][1]+=2)for(a[4][2]=0;a[4][2]<=9;a[4][2]++)for(a[2][4]=0;a[2][4]<=9;a[2][4]++){a[1][5]=n-a[5][1]-a[4][2]-a[3][3]-a[2][4];if(a[1][5]<0||a[1][5]>9)continue;if(a[1][5]%2==0)continue;tmp=a[5][1]*10000+a[4][2]*1000+a[3][3]*100+a[2][4]*10+a[1][5];if(ss[tmp]==1)continue;for(a[3][1]=1;a[3][1]<=9;a[3][1]++)for(a[3][2]=0;a[3][2]<=9;a[3][2]++)for(a[3][4]=0;a[3][4]<=9;a[3][4]++){a[3][5]=n-a[3][1]-a[3][2]-a[3][3]-a[3][4];if(a[3][5]<0||a[3][5]>9)continue;if(a[3][5]%2==0)continue;tmp=a[3][1]*10000+a[3][2]*1000+a[3][3]*100+a[3][4]*10+a[3][5];if(ss[tmp]==1)continue;for(a[2][1]=1;a[2][1]<=9;a[2][1]++){a[4][1]=n-a[1][1]-a[2][1]-a[3][1]-a[5][1];if(a[4][1]<=0||a[4][1]>9)continue;tmp=a[1][1]*10000+a[2][1]*1000+a[3][1]*100+a[4][1]*10+a[5][1];if(ss[tmp]==1)continue;for(a[2][5]=1;a[2][5]<=9;a[2][5]+=2){a[2][3]=n-a[2][1]-a[2][2]-a[2][4]-a[2][5];if(a[2][3]<0||a[2][3]>9)continue;tmp=a[2][1]*10000+a[2][2]*1000+a[2][3]*100+a[2][4]*10+a[2][5];if(ss[tmp]==1)continue;a[4][5]=n-a[1][5]-a[2][5]-a[3][5]-a[5][5];if(a[4][5]<0||a[4][5]>9)continue;if(a[4][5]%2==0)continue;tmp=a[1][5]*10000+a[2][5]*1000+a[3][5]*100+a[4][5]*10+a[5][5];if(ss[tmp]==1)continue;a[4][3]=n-a[4][1]-a[4][2]-a[4][4]-a[4][5];if(a[4][3]<0||a[4][3]>9)continue;tmp=a[4][1]*10000+a[4][2]*1000+a[4][3]*100+a[4][4]*10+a[4][5];if(ss[tmp]==1)continue;for(a[5][2]=1;a[5][2]<=9;a[5][2]+=2){a[1][2]=n-a[2][2]-a[3][2]-a[4][2]-a[5][2];if(a[1][2]<=0||a[1][2]>9)continue;tmp=a[1][2]*10000+a[2][2]*1000+a[3][2]*100+a[4][2]*10+a[5][2];if(ss[tmp]==1)continue;for(a[5][3]=1;a[5][3]<=9;a[5][3]+=2){a[1][3]=n-a[2][3]-a[3][3]-a[4][3]-a[5][3];if(a[1][3]<=0||a[1][3]>9)continue;tmp=a[1][3]*10000+a[2][3]*1000+a[3][3]*100+a[4][3]*10+a[5][3];if(ss[tmp]==1)continue;a[1][4]=n-a[1][1]-a[1][2]-a[1][3]-a[1][5];if(a[1][4]<=0||a[1][4]>9)continue;tmp=a[1][1]*10000+a[1][2]*1000+a[1][3]*100+a[1][4]*10+a[1][5];if(ss[tmp]==1)continue;a[5][4]=n-a[1][4]-a[2][4]-a[3][4]-a[4][4];if(a[5][4]<0||a[5][4]>9)continue;if(a[5][4]%2==0)continue;tmp=a[1][4]*10000+a[2][4]*1000+a[3][4]*100+a[4][4]*10+a[5][4];if(ss[tmp]==1)continue;tmp=a[5][1]*10000+a[5][2]*1000+a[5][3]*100+a[5][4]*10+a[5][5];if(ss[tmp]==1)continue;sum++;for(i=1;i<=5;i++)for(j=1;j<=5;j++)ans[sum][i*5-6+j]=a[i][j]+'0';}}}}}}}for(i=1;i<=sum-1;i++)for(j=i+1;j<=sum;j++)if(strcmp(ans[i],ans[j])>0){strcpy(ans[0],ans[i]);strcpy(ans[i],ans[j]);strcpy(ans[j],ans[0]);}for(i=1;i<=sum;i++){for(j=1;j<=5;j++){for(k=1;k<=5;k++)fprintf(fo,"%c",ans[i][j*5-6+k]);fprintf(fo,"/n");}if(i!=sum)fprintf(fo,"/n");}fclose(fi);fclose(fo);return 0;}