【欧拉】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;
}