比赛 noip 评测结果 AAAAAAAAAAAAAATTTTTT
题目名称 __完全平方数 最终得分 70
用户昵称 liuliuliu 运行时间 12.880 s
代码语言 C++ 内存使用 34.58 MiB
提交时间 2016-11-04 19:57:14
显示代码纯文本
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <climits>
#include <cstdlib>
#include <ctime>
#include <sstream>
#include <queue>
#include <set>
#include <map>
using  namespace std;

const int maxn=5000000+10;
const int maxm=1000+10;
const int mod=100000007;
int n,m;
int f[maxn],cnt[maxn];;
long long ans=1;
bool a[maxn];
int len;

void work(int x)
{
	int p=0;
	while(x!=1)
	{
		if(x%f[p]==0)
		{
			cnt[ f[p] ]++;
			x/=f[p];
			ans=(ans*f[p])%mod;
		}
		else
			p++;
	}
}
int main()
{
	freopen("xnumber.in","r",stdin);
	freopen("xnumber.out","w",stdout); 
	cin>>n;
	for (int i=2; i<=n; i++)  //欧拉筛 
	{
		if (!a[i])
		{
			f[len++]=i;
		}
		for (int j=0; j<len && i*f[j]<=n; j++)
		{
			a[i*f[j]]=true;
			if (!(i%f[j])) break;
		}
	}
	for(int i=n; i>=2;i--)
	{
		if( a[i] )
			work(i);
		else if(!a[i])
		{
			if( cnt[ i ]%2 ==1)
				ans=(ans*i)%mod;
			}	
	}
	cout<<ans<<endl;
	return 0;
}