蒲丰投针,其实这算法本质是个求定积分的算法
[C] 纯文本查看 复制代码
#include <windows.h>
#define N 10000
#define C 20
#define R 19
__int64 a=1000000, b=0, c=N*R, d, e, f[N*R+1], g;
main()
{
for (;b-c;) f[b++]=C*N;
for(;d=0,g=c*2;c-=R,printf("%.6d",e+d/a),e=d%a)
for(b=c;d+=f[b]*a,f[b]=d%--g,d/=g--,--b;d*=b);
system("pause");
return 0;
} |