YYHSOI模拟赛题解（T4完全平方数）

　　

题目描述

一个数如果是另一个整数的完全平方，那么我们就称这个数为完全平方数（Pefect Sqaure），也称平方数。

小A认为所有的平方数都是很perfect的~

于是他给了小B一个任务：用任意个不大于n的不同的正整数相乘得到完全平方数，并且小A希望这个平方数越大越好。

请你帮助小B告诉小A满足题意的最大的完全平方数。

输入

输入文件名为number.in 

 

输入仅 1行，一个数n。

 

输出

输出文件名为number.out

输出仅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。

这一道题目，我一看就觉得是质因数分解，但是由于数据范围是500W，所以对于一般的筛法就不可行了，所以这一道题目需要用到欧拉线性筛来筛出所有的质数。但是筛完后对于分解质因数仍然没有帮助。然后我就想到了初中数学所学过的，如果对于一个范围内的数质因数分解后存在几个P（P为质数），我们有一个求和式，我就把这个用上去了。然后这道题目的时间复杂度就变成了约为2000W的一个算法

#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>

using namespace std;

int N;
long long Mod=100000007;
bool flag[5000005];
int p[5000005];
long long s[5000005];
int num;
int top=0;
int q[5000005];

long long mul(int a,long long b,long long Mod)
{
    long long tmp=1;
    long long A=(long long) a;
    while (b)
    {
        if (b & 1) tmp=((tmp % Mod) * (A % Mod)) % Mod;
        A=(A % Mod) * (A % Mod) % Mod;
        b=b >> 1;
    }
    return tmp;
}

int main()
{
    scanf("%d",&N);
    for (int i=2; i<=N; i++)
    {
        if (! flag[i])
        {
            p[++num]=i;
        }
        for (int j=1; p[j] * i <= N && j<=num; j++)
        {
            flag[p[j]*i]=1;
            if (i % p[j] == 0) break;
        }
    }
    /*for (int i=2; i<=N; i++)
    {
        if (! flag[i])
        {
            continue;
        }
        int X=i;
        for (int j=1; j<=num; j++)
        {
            while (X % p[j] == 0)
            {
                s[p[j]]++;
                X=X / p[j];
            }
            if (X == 0) break;
        }
    }*/
    for (int i=1; i<=num; i++)
    {
        for (long long j=p[i]; j<= N; j*=p[i])
        {
            s[p[i]]+=(long long) N / j;
        }
    }
    long long ans=1;
    for (int j=1; j<=num; j++)
    {
        //printf("%d %d\n",p[j],s[p[j]]);
        int i=p[j];
        long long b=s[i];
        if (b & 1) b--;
        if (b == 0) continue;
        long long tmp=mul(i,b,Mod);
        ans=(ans % Mod) * (tmp % Mod) % Mod;
    }
    printf("%lld\n",ans);
    return 0;
}