HDU-4196 Remoteland（数论）

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4196

题意：给出一个正整数n，求出一个最大的数，这个数满足是由[1,n]之间不同的数相乘得到，并且这个数是完全平方数。

分析：乘积最大的数是n!，因为n!可以拆成一堆素数的乘积，只需要剔除n!里面多余的素数即可。如果n!里面包含的素数因子的指数是奇数，则需要剔除一个这样的素数才能构成平方数，是偶数就不用剔除，找素数因子的幂要用到数论里面的勒让德定理。

  因为最终得到的结果是模之后的结果，所以如果先求阶乘再除以需要剔除的素数实现比较麻烦（要用到逆元），这里采用一个巧妙的处理：在筛选素数时预处理阶乘，当这个数是素数先不用乘到阶乘中去，最后判断幂时如果是偶次幂再乘到结果中去。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<string>
#include<algorithm>
using namespace std;
typedef long long ll;
const int mod=1e9+7;
const int maxn=1e7+5;
bool ispri[maxn];
int primes[maxn];
ll fac[maxn];
int k;

void pre()//预处理操作，筛素数和预处理阶乘
{
    memset(ispri,true,sizeof ispri);
    ispri[0]=ispri[1]=false;
    fac[0]=fac[1]=1;
    for(int i=2;i<maxn;i++)
    {
        fac[i]=fac[i-1];
        if(ispri[i])
        {
            primes[k++]=i;
            for(int j=i+i;j<maxn;j+=i)
                ispri[j]=false;
        }
        else fac[i]=(fac[i]*i)%mod; //素数暂时不乘到阶乘中去，fac[i]只保存i!的合数部分
    }
}
int main()
{
    pre();
    int n;
    while(~scanf("%d",&n)&&n)
    {
        ll ans=fac[n];
        for(int i=0;i<k&&primes[i]<=n;i++)
        {
            int cnt=0,t=n;
            while(t)
            {
                t/=primes[i];
                cnt+=t;
            }             //cnt所求的是primes[i]的幂指数，这里用到的是勒让德定理
            if(cnt%2==0)   ans=(ans*primes[i])%mod; //如果素数i的幂指数是偶数就乘到结果中去
        }
        printf("%d\n",(int)ans);
    }
    return 0;
}