题目链接:
http://acm.hdu.edu.cn/showproblem.php?pid=2582
题目大意:
给出公式Gcd(n)=gcd(C[n][1],C[n][2],……,C[n][n-1]),让求f(n)= Gcd(3)+Gcd(4)+…+Gcd(i)+…+Gcd(n)。
解题思路:
对于公式:,小范围打表可以找出规律:
(1),如果n为素数,那么G=n;
(2),如果n有多个素因子,那么G=1;
(3),如果n只有一个素因子,那么G=该素因子。
所以可以直接利用筛素数的模板,筛出一个素数时就往后更新答案即可
1 #include<bits/stdc++.h> 2 using namespace std; 3 typedef long long ll; 4 const int maxn = 1000000+10; 5 ll ans[maxn]; 6 bool is_prime[maxn]; 7 int sieve(int n)//返回n以内素数的个数 8 { 9 int p = 0; 10 for(int i = 0; i <= n; i++)is_prime[i] = ans[i] = 1; 11 is_prime[0] = is_prime[1] = 0; 12 for(ll i = 2; i <= n; i++) 13 { 14 if(is_prime[i]) 15 { 16 ans[i] = i; 17 for(ll j = i * 2; j <= n; j += i)is_prime[j] = 0; 18 for(ll j = i * i; j <= n; j *= i)ans[j] = i; 19 } 20 } 21 for(int i = 4; i <= n; i++)ans[i] += ans[i - 1]; 22 return p; 23 } 24 int main() 25 { 26 int tot = sieve(1000000), n; 27 while(scanf("%d", &n) != EOF) 28 printf("%lld\n", ans[n]); 29 return 0; 30 }