cpp.rar_tyvj_tyvj 1268 c

#include "stdio.h" #include "math.h" int a[10000],c[10000];//9592个质数，a是质数组，c计次数。 int r[5000000];//结果 long long time=1; int main() { int i,j,f=1,sq,k=1; int b,d,n,m,p,w=0; a[0]=2; scanf("%d",&n); b=2*n; for(i=3;i<=b;i+=2) { sq=sqrt(i); for(j=3;j<=sq;j+=2) if(i%j==0) goto next; a[k++]=i; next:; }//产生需要的质数 for(i=0;a[i]<b&&a[i];i++) { m=b; while(m!=0) { m/=a[i]; c[i]+=m; } } p=i+1; b=n; for(i=0;a[i]<b;i++) { m=b; while(m!=0) { m/=a[i]; c[i]-=m; } } if(a[i]==b) c[i]-=1; b=n+1; for(i=0;a[i]<b;i++) { m=b; while(m!=0) { m/=a[i]; c[i]-=m; } } if(a[i]==b) c[i]-=1; r[0]=1; //分解完成 for(i=0;i<10000;i++) { for(j=0;j<c[i];j++) { if(time*a[i]<200000) time*=a[i]; else { for(k=0;k<=w||r[k]!=0;k++) r[k]*=time; for(k=0;k<=w||r[k]!=0;k++) { r[k+1]+=r[k]/1000; r[k]%=1000; } w+=2; time=a[i]; } } } for(k=0;k<=w||r[k]!=0;k++) r[k]*=time; for(k=0;k<=w||r[k]!=0;k++) { r[k+1]+=r[k]/1000; r[k]%=1000; } w+=2; //连乘 while(r[w]==0) w--; printf("%d",r[w]); w--; for(i=w;i>=0;i--) printf("%03d",r[i]); while(getchar()!=EOF); //输出 return 0; }