HDU 2582 f(n) 数学

题意:求f(n)=GCD(3)+GCD(4)+….GCD(n)
GCD(n)=gcd(C(1,n) , C(2,n) , … C(n-1,n) )
C(1,n)为组合数里的从n个东西中挑出1个东西的方案数
思路:
每个数都能表示为

如果该数是素数很容易得出GCD(n)=n;
否则的话有两种情况
1. n=p1^e1
2. n=p1^e1 * p2^e2 * …pk ^ek
如果是情况二
例如是12=2^2 * 3^1
我们可以通过因为C(1,12)==12，所以我们的答案只能是12的因子
我们可以通过C（4，12）来得到一个数x1，该数一定是gcd(2^2, x1)==1
我们可以通过C（3，12）来得到一个数x2，该数一定是gcd(3, x2)==1
PS:可以试着手推一下
所以gcd(C(1,12),C(4,12),C(3,12))==1
以此推广为方案 2 GCD(n)==1
情况一的情况下
例如是8=2^3
因为C（1，8）==2，所以只能是8的因子
我们可以通过C（4，8）来得到一个数x1，x1一定=2^1 * k1| gcd(k1,2)==1
以此推广方案一 一定是GCD(n)=p1(p1为n的素因子)

#include<stdio.h>
#include<string.h>
#include<iostream>
#include<algorithm>
#include<math.h>
#include<queue>
#include<stack>
#include<string>
#include<vector>
#include<map>
#include<set>
using namespace std;
#define mem(a,b) memset(a,b,sizeof(a))
#define lowbit(x) (x&(-x))
typedef long long LL;
#define maxn 1000005
const int inf=(1<<28)-1;
bool notprime[maxn];
int primes[100005];
int ans[maxn];
LL sum[maxn];
void get_prime()
{
    notprime[1]=true;
    for(int i=2;i<maxn;++i)
    if(!notprime[i])
    {
        primes[++primes[0]]=i;
        for(LL j=(LL)i*i;j<maxn;j+=i)
        notprime[j]=true;
    }
}
int main()
{
    get_prime();
    for(int i=3;i<maxn;++i)
    ans[i]=1;
    for(LL i=1;i<=primes[0];++i)
    {
        ans[primes[i]]=primes[i];
        for(LL j=primes[i];j<maxn;j*=primes[i])
        ans[j]=primes[i];
    }
    sum[3]=ans[3];
    for(int i=4;i<maxn;++i)
    sum[i]=sum[i-1]+ans[i];
    int n;
    while(~scanf("%d",&n))
    {
        printf("%lld\n",sum[n]);
    }
    return 0;
}