hdu4196 remoteland

此题开始有点小纠结，一直不知道为什么要打合数表。。。后面想了一下觉得是：

一个数用n表示最大就是N! 那么N!可以用素数表示为a1^k1*a2^k2*....*an^kn其中ak为质数

要是完全平方数那么ki必为偶数，如果k为奇数-1就可以了。

如果ki为奇数那么统计的时候-1采用，如果是偶数就采用。

现在统计合数的意义，是因为把质数幂次-1，所以就相当于N！去掉质数 = a1^（k1-1）*a2^（k2-1）*....*an^（kn-1）

在main中的while统计n中每个质数的幂次数，那么如果是偶数则表明我k-1的时候少乘了一个质数，那么乘上就可以了。

如果是奇数表明我是需要-1的幂次不需要统计了

#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
#define N 10000002
const int mod = 1000000007;
 
bool a[10000002]={1,1,0};
int p[680000],num=0;
int sum[10000002]={0,1};
void prime(){
    int i,j;
    for(i=2;i<3164;i++)
        if(!a[i]) 
            for(j=i*i;j<N;j+=i)
                a[j]=1;
    for(__int64 i=2;i<N;i++)
        if(!a[i])  p[num++]=i, sum[i]=sum[i-1];
        else sum[i]=sum[i-1]*i%mod;  //保存i之前所有合数的积
}
int main(){
    prime();
    int i,j,k,temp;
    __int64 n,ans;
    while(scanf("%I64d",&n)&&n){
        ans=sum[n];
        for(i=0;p[i]*2<=n;i++){
            int k=0;
            for(temp=n;temp;)
                k+=(temp/=p[i]);
            if(!(k&1)) ans=p[i]*ans%mod;
        }
        printf("%I64d\n",ans);
    }
}