[来源未知]完全平方数

内存128M,时间1s
完全平方数
【问题描述】
一个数如果是另一个整数的完全平方，那么我们就称这个数为完全平方数（Pefect Sqaure），也称平方数。
小A认为所有的平方数都是很perfect的。
于是他给了小B一个任务：用任意个不大于n的不同的正整数相乘得到完全平方数，并且小A希望这个平方数越大越好。
请你帮助小B告诉小A满足题意的最大的完全平方数。
【输入格式】
输入仅 1行，一个数n。
【输出格式】
输出仅1行，一个数表示答案。由于答案可以很大，所以请输出答案对100000007取模后的结果。
【输入输出样例1】
number.in
7
number.out
144
【输入输出样例解释1】
144=2×3×4×6，是12的完全平方。
【输入输出样例2】
number.in
9
number.out
5184
【输入输出样例解释2】
5184=3×4×6×8×9，是72的完全平方。
【数据范围】
对于20%的数据，0 < n≤100；
对于50%的数据，0 < n≤5,000；
对于70%的数据，0 < n≤100,000；
对于100%的数据，0 < n≤5,000,000。

sol:
对于一个完全平方数显然其分解质因数后每个质因数的次幂为偶数，我们将1到n的数分解质因数后，将质因数个数为奇数的-1(相当于去掉了若干个互不相同的质数)。显而易见的这可以得到一个最大的完全平方数。
因为将每个数质因数分解的复杂度比较大，考虑对于每个质，有多少个数包含他，显然1-n中有n/prime[j]个数有prime[j]这个质数，包含2次的为n/prime[j]/prime[j]个。这个可以用分配律来考虑。

#include<cstdio>
#include<algorithm>
#include<string>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<iostream>
using namespace std;
typedef long long ll;
const int N=5000007;
const int pyz=1e8+7;
int n,m;
inline int read()
{
    char c;
    int res,flag=0;
    while((c=getchar())>'9'||c<'0') if(c=='-')flag=1;
    res=c-'0';
    while((c=getchar())>='0'&&c<='9') res=(res<<3)+(res<<1)+c-'0';
    return flag?-res:res;
}
bool is[N];
int ans,tmp,cnt[N],prime[N];
inline int ksm(int s,int t)
{
    int res=1;
    while(t)
    {
        if(t&1) res=(ll)res*s%pyz;
        s=(ll)s*s%pyz;
        t>>=1;
    }
    return res;
}
int main()
{
    freopen("number.in","r",stdin);
    freopen("number.out","w",stdout);
    n=read();
    for(int i=2;i<=n;++i)
    {
        if(!is[i]) prime[++prime[0]]=i;
        for(int j=1;j<=prime[0]&&i*prime[j]<=n;++j)
        {
            is[i*prime[j]]=1;
            if(!(i%prime[j])) break;
        }
    }
    for(int i=1;i<=prime[0];++i)
    {
        tmp=n;
        while(tmp/prime[i]>0)
        {
            cnt[i]+=tmp/prime[i];
            tmp/=prime[i];
        }
    }
    ans=1;
    for(int i=1;i<=prime[0];++i)
    {
        if(cnt[i]&1) cnt[i]--;
        ans=(ll)ans*ksm(prime[i],cnt[i])%pyz;
    }
    printf("%d",ans);
}