记录编号 338064 评测结果 AAAAAAAAAAAAAAAAAAAA
题目名称 __完全平方数 最终得分 100
用户昵称 Gravatarliuliuliu 是否通过 通过
代码语言 C++ 运行时间 0.603 s
提交时间 2016-11-05 07:32:58 内存使用 7.56 MiB
显示代码纯文本
#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=1000+10;
const int maxm=1000+10;
#define debug(x) cerr<<#x<<"="<<x<<endl
#define INF 0x7f7f7f7f
using namespace std;
typedef pair<int,int> pii;
typedef long long ll;
inline int init()
{
	int now=0,ju=1;char c;bool flag=false;
	while(1)
	{
		c=getchar();
		if(c=='-')ju=-1;
		else if(c>='0'&&c<='9')
		{
			now=now*10+c-'0';
			flag=true;
		}
		else if(flag)return now*ju;
	}
}
inline long long llinit()
{
	long long now=0,ju=1;char c;bool flag=false;
	while(1)
	{
		c=getchar();
		if(c=='-')ju=-1;
		else if(c>='0'&&c<='9')
		{
			now=now*10+c-'0';
			flag=true;
		}
		else if(flag)return now*ju;
	}
}
int prime[1000001];
bool isprime[5000001];
int cnt=0;
int n;
const int mod=100000007;
ll quick_pow(ll base,ll power)
{
	ll ret=1;
	for(;power;power>>=1)
	{
		if(power&1)
		{
			ret=((ret%mod)*(base%mod))%mod;
		}
		base=((base%mod)*(base%mod))%mod;
	}
	return ret;
}
ll solve(int x)
{
	ll y=n,z=0;
	while(y)
	{
		y/=x;
		z+=y;
	}
	if(z&1)
	{
		return quick_pow(x,z-1);
	}
	else return quick_pow(x,z);
}
ll ans=1;
void euler()
{
	memset(isprime,1,sizeof(isprime));
	for(int i=2;i<=n;i++)
	{
		if(isprime[i])
		{
			prime[++cnt]=i;
			ans=((ans%mod)*(solve(i)%mod))%mod;
		}
		for(int j=1;j<=cnt;j++)
		{
			if(prime[j]*i>n)break;
			isprime[prime[j]*i]=false;
			if(i%prime[j]==0)break;
		}
	}
}
int main()
{
	freopen("xnumber.in","r",stdin);
	freopen("xnumber.out","w",stdout);
	int p;
	n=init();
	euler();
	printf("%d\n",ans);
	return 0;
}