hdu_2582 f(n) (找规律)

http://acm.hdu.edu.cn/showproblem.php?pid=2582

题意：
f(n)= Gcd(3)+Gcd(4)+…+Gcd(i)+…+Gcd(n).
Gcd(n)=gcd(C[n][1],C[n][2],……,C[n][n-1])
C[n][k]从N个里面取k个
思路：
开始推公式，没推出，遂YY了下，打了个表，发现了规律：
Gcd(n)可以直接算出：
1,当n为质数时，Gcd(n)=n;
2,当n为某个质数k的q次方时（q>=2） Gcd(n)=k;
3,其它情况都为1。
根据以上规律，我先把所有的数+1,再把质数或质数的次方加上质数本身-1 

我的代码：

/*
program hdu_2582
author:BlackAndWhite
*/ 
#include<stdio.h>
#include<math.h>
#include<string.h>
__int64 n;
__int64 i,ans;
bool is[1000005];
int prm[1000005];
int getprm(int n)
{
    int i, j, k = 0;

    int s, e = (int)(sqrt(0.0 + n) + 1);
    memset(is, 1, sizeof(is));
    prm[k++] = 2; is[0] = is[1] = 0;
    for (i = 4; i < n; i += 2) is[i] = 0;
    for (i = 3; i < e; i += 2) if (is[i])
    {
        prm[k++] = i;
        for (s = i * 2, j = i * i; j < n; j += s)
            is[j] = 0;
    }
    for ( ; i < n; i += 2) if (is[i]) prm[k++] = i;
    return k;
}
__int64 fun(int x)
{
    __int64 i,z=x;
    for(i=0;;i++)
        if(z<=n) z*=x;
        else break;
    return (i)*(x-1);
}
int main()
{
    getprm(1000001);
    while(~scanf("%I64d",&n))
    {
        ans=n-2;
        for(i=0;prm[i];i++)//从2加起
        {
            if(prm[i]>n) break;
            ans+=fun(prm[i]);
        }
        printf("%I64d\n",ans-1);
    }
    return 0;
}