【脑洞系列】胡策思维题

【问题描述】
从1 − N中找一些数乘起来使得答案是一个完全平方数，求这个完全平方数最大可能是多少。
【输入格式】
第一行一个数字N。
【输出格式】
一行一个整数代表答案对100000007取模之后的答案。
【样例输入】
7
【样例输出】
144
【数据规模与约定】
对于20%的数据，1 ≤ N ≤ 100。
对于50%的数据，1 ≤ N≤ 5000。
对于70%的数据，1 ≤ N ≤ 10^5
对于100%的数据，1 ≤N ≤ 5 × 10^6 。

这题，说实话真不好想。。。反正考试的时候我没想出来。。
思路：因为是完全平方数，所以这个数的因子在1~n中一定是出现了偶数次，道理显然，否则一定不会出现这样一个数。因此我们可以对于N分解质因数，用线性筛O(N)求出1~n的素数。求出素数后用N除以每个素数，开一个数组记录一下出现的次数，对于除出来的商，我们在判断一下能否在继续除以当前的素数，能继续除就继续加，道理显然。
因此，如果数组中记录的数为奇数，就-1，这样就能保证以上偶数次的要求。

例如样例：
n=7，从2、3、4、5、6、7中选数
2=2^1；
3=3^1；
4=2^2；
5=5^1；
6=2^1*3^1；
7=7^1；
所以质因数2的幂是4,3的幂是2；
ans=2^4*3^2==144；
代码见下：

#include<iostream>
#include<cstdio>
#include<cstdlib>
using namespace std;
const int maxn=5000005;
long long prime[maxn];
bool not_prime[maxn];
long long num[maxn];
long long mod=100000007;
long long ksm(long long a,long long b)
{
    if(b==0) return 1ll;
    if(b==1) return a%mod;
    long long tmp=ksm(a,b/2)%mod;
    if(b%2==0)
        return ((tmp%mod)*(tmp%mod))%mod;
    else
        return ((((tmp%mod)*tmp)%mod)*(a%mod))%mod;
}
int main()
{
    freopen("hao.in","r",stdin);
    freopen("hao.out","w",stdout);
    long long n;
    int cnt=0;
    long long ans=1;
    scanf("%I64d",&n);
    not_prime[1]=true;
    for(long long i=2;i<=n;i++)
    {
        if(!not_prime[i])
            prime[++cnt]=i;
        for(int j=1;j<=cnt;j++)
        {
            if(prime[j]*i>n) break;
            not_prime[prime[j]*i]=true;
            if(i%prime[j]==0) break;
        }
    }
    for(int i=1;i<=cnt;i++)
    {
        long long aa=n;
        while(aa!=0)
        {
            num[i]+=aa/prime[i];
            aa/=prime[i];
        }
    }
    for(int i=1;i<=cnt;i++)
    {
        if(num[i]%2==0)
        ans=(ans*ksm(prime[i],num[i]))%mod;
        else
        ans=(ans*ksm(prime[i],num[i]-1))%mod;       
    }
    printf("%I64d",ans);
    return 0;
}