Hao

【问题描述】
从中找一些数乘起来使得答案是一个完全平方数，求这个完全平方数最大可能是多少。
【输入格式】
第一行一个数字。
【输出格式】
一行一个整数代表答案对取模之后的答案。
【样例输入】
7
【样例输出】
144
【样例解释】
但是塔外面有东西。
0=< N< =5*1E6

---

思路：想让一个数为完全平方数，那么这个数分解质因数后每个质因数的次数一定是偶数次个。
反其道而行之，对n的阶乘分解质因数
如果一个质因数出现的次数为偶数次，直接加入答案
如果是奇数，需要去掉一个再加入答案。
如何快速得到一个质数在n!出现的次数？
统计1~N的所有因数个数,实际上可以发现1的因数个数=N/1,2的因数个数=N/2…那么则有如下式子:

F(N)=N/1+N/2+…+N/N。
那么这个问题就便于解决了,举个例子,如求1~81的3的个数是多少,可以发现ans=81/3+27/3+9/3+3/3=40。233333
工具：快速幂，分解质因数，线性筛（逃）！

#include <cstdio>
#include <iostream>
#define ll long long
using namespace std;
const int maxm=5*1e6+10;
const int mod=100000007;
ll prime[maxm],cnt;
bool vis[maxm];
ll num[maxm];
void get_p(int n)
{
    vis[1]=1;
    for(int i=2;i<=n;i++)
    {
        if(!vis[i])
         prime[++cnt]=i;

        for(int j=1;j<=cnt&&1ll*i*prime[j]<=n;j++)
        {
            vis[i*prime[j]]=1;
            if(!(i%prime[j]))
             break;
        }
    }
}
void rts(int x,int n)
{
    while(n)
     num[x]+=n/prime[x],n/=prime[x];
}
ll fast_pow(ll x,ll y)
{
    ll tmp=1;
    while(y)
    {
        if(y&1)
        tmp=tmp*x%mod;
        x=x*x%mod;
        y>>=1;
    }
    return tmp;
}
int main()
{
    freopen("hao.in","r",stdin);
    freopen("hao.out","w",stdout);
    int n;
    scanf("%d",&n);

    get_p(n);

    for(int i=1;i<=cnt;i++)
     rts(i,n);

    ll ans=1;
    for(int i=1;i<=cnt;i++)
    {
        if(num[i]%2) num[i]--;
        ans=(1ll*ans*1ll*fast_pow(prime[i],num[i]))%mod;
    }

    printf("%lld",ans); 
}