【欧拉】HDU 3501 Calculation 2

= =本来我们应该是讲搜索专题的结果不知道为啥出来了道这个。。。容斥原理不会，所以我用欧拉函数做的

题目地址：http://acm.hdu.edu.cn/showproblem.php?pid=3501

欧拉函数的话咱其实也说不清。。。背板来的，代码如下

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#define LL long long

using namespace std;

LL prime[100000];

void init()
{
memset(prime,0,sizeof prime);
int cont=0;
for(int i=2;i<=100000;i++)
{
if(prime[i]!=0) continue;
for(int j=i;j<=100000;j+=i)
{
prime[j]=1;
}
prime[cont++]=i;
}
}

LL n;

int main()
{
init();        //筛法取素数
//for(int i=0;i<10;i++) cout<<prime[i]<<endl;

while(~scanf("%lld",&n) && n)
{
LL m=n,ans=n*n;
for(int i=0;prime[i]*prime[i]<=m;i++)
{
if(m%prime[i]==0)
{
while(m%prime[i]==0)
{
m/=prime[i];
}
ans=ans/prime[i]*(prime[i]-1);
}
}
if(m!=1) ans=ans/m*(m-1);
ans=((n-1)*n/2-ans/2)%1000000007;

printf("%lld\n",ans);
}
return 0;
}