UVA 11426 GCD - Extreme (II)(欧几里得定理+欧拉函数)
2017-05-03 11:07
387 查看
解这道题,需要以下几步:
1.建立递推关系,s(n)=s(n-1)+gcd(1,n)+gcd(2,n)+……+gcd(n-1,n);
2.设f(n)=gcd(1,n)+gcd(2,n)+……+gcd(n-1,n)。
gcd(x,n)=i是n的约数(x<n),按照这个约数进行分类。设满足gcd(x,n)=i的有g(n,i)个,则有f(n)=sum(i*g(n,i))。
而gcd(x,n)=i等价于gcd(x/i,n/i)=1,因此g(n,i)等价于phi(n/i)
phi(x)为欧拉函数。
3.降低时间复杂度。用筛法预处理phi[x]表
用筛法预处理f(x)->枚举因数,更新其所有倍数求解。
代码如下
1.建立递推关系,s(n)=s(n-1)+gcd(1,n)+gcd(2,n)+……+gcd(n-1,n);
2.设f(n)=gcd(1,n)+gcd(2,n)+……+gcd(n-1,n)。
gcd(x,n)=i是n的约数(x<n),按照这个约数进行分类。设满足gcd(x,n)=i的有g(n,i)个,则有f(n)=sum(i*g(n,i))。
而gcd(x,n)=i等价于gcd(x/i,n/i)=1,因此g(n,i)等价于phi(n/i)
phi(x)为欧拉函数。
3.降低时间复杂度。用筛法预处理phi[x]表
用筛法预处理f(x)->枚举因数,更新其所有倍数求解。
代码如下
#include<cstdio> #include<cstring> #include<cstdlib> #include<cmath> #include<ctime> #include<iostream> #include<algorithm> #include<vector> #include<string> #include<queue> using namespace std; typedef long long ll; const int MAXN = 4000010; bool check[MAXN+10]; int phi[MAXN+10]; int prime[MAXN+10]; int tot; void phi_and_prime_table(int N) { memset(check,false,sizeof(check)); phi[1] = 1; tot = 0; for(int i = 2; i <= N; i++) { if( !check[i] ) { prime[tot++] = i; phi[i] = i-1; } for(int j = 0; j < tot; j++) { if(i * prime[j] > N)break; check[i * prime[j]] = true; if( i % prime[j] == 0) { phi[i * prime[j]] = phi[i] * prime[j]; break; } else { phi[i * prime[j]] = phi[i] * (prime[j] - 1); } } } } ll f[MAXN],s[MAXN]; int main() { phi_and_prime_table(MAXN); memset(f,0,sizeof(f)); for(int i=1; i<=MAXN; i++) for(int j=i+i; j<=MAXN; j+=i) f[j]+=i*phi[j/i]; memset(s,0,sizeof(s)); s[1]=0; for(int i=2; i<=MAXN; i++) s[i]=s[i-1]+f[i]; int n; while(~scanf("%d",&n)&&n) { printf("%lld\n",s ); } return 0; }
相关文章推荐
- UVa 11426 GCD - Extreme (II) / 素数筛选 + 欧拉函数
- 【欧拉函数 && 求gcd(n, x) = y,x的个数】UVA - 11426 GCD - Extreme (II)
- UVa 11426 - GCD - Extreme (II) 转化+筛法生成欧拉函数表
- UVA - 11426 GCD - Extreme (II) (欧拉函数)
- GCD - Extreme (II) UVA - 11426 (欧拉函数)
- UVA 11426 GCD - Extreme (II) (欧拉函数)
- UVA11426:GCD - Extreme (II) (欧拉函数)
- 【欧拉函数】uva 11426 GCD - Extreme (II)
- UVA - 11426 GCD - Extreme (II) (欧拉函数)
- UVa 11426 (欧拉函数 GCD之和) GCD - Extreme (II)
- UVA 11426 GCD - Extreme (II) (欧拉函数)
- GCD - Extreme (II) UVA - 11426 (欧拉函数)
- UVa 11426 GCD - Extreme (II) (欧拉函数应用·O(N*logN))
- UVA 11426 GCD - Extreme (II) (数论|欧拉函数)
- UVA11426 GCD - Extreme (II)---欧拉函数的运用
- UVA 11426 GCD Extreme (II) 欧拉函数
- UVA 11426 GCD - Extreme (II) (欧拉函数+筛法)
- UVa 11426 - GCD - Extreme (II) (数学 欧拉函数)
- UVA 11426 - GCD - Extreme (II) 欧拉函数-数学
- UVA 11426 GCD-Extreme(II) ★ (欧拉函数)