UVA12995

U V A 12995 UVA12995

  • 题目
    UVA12995
  • 题意
  • 分析
    g c d ( a , b ) = = 1 gcd(a,b) == 1 考虑欧拉函数
    由题意得:
    i = 1 n j = 1 i 1 g c d ( i , j ) = 1 \sum_{i=1}^{n}\sum_{j=1}^{i-1}gcd(i,j) = 1
    \Rightarrow i = 1 n φ ( i ) \sum_{i=1}^{n}\varphi(i)
    欧拉函数线性筛预处理就可以了。
  • 代码
/*
独立思考
*/
#include<bits/stdc++.h>
using namespace std;
template <typename T>
void read(T &x)
{
	x = 0;
	char c = getchar();
	int sgn = 1;
	while (c < '0' || c > '9') {if (c == '-')sgn = -1; c = getchar();}
	while (c >= '0' && c <= '9')x = x * 10 + c - '0', c = getchar();
	x *= sgn;
}
template <typename T>
void out(T x)
{
	if (x < 0) {putchar('-'); x = -x;}
	if (x >= 10)out(x / 10);
	putchar(x % 10 + '0');
}
typedef long long ll;
typedef unsigned long long ull;
const int N = 1e6 + 5;
int v[N], prime[N], phi[N];
int m = 0;
void euler(int n)
{
	memset(v, 0, sizeof(v));
	for (int i = 2; i <= n; i++)
	{
		if (v[i] == 0)
		{
			v[i] = i;
			prime[++m] = i;
			phi[i] = i - 1;
		}
		for (int j = 1; j <= m; j++)
		{
			if (prime[j] > v[i] || prime[j] > n / i) break;
			v[i * prime[j]] = prime[j];
			phi[i * prime[j]] = phi[i] * (i % prime[j] ? prime[j] - 1 : prime[j]);
		}
	}
}
ll sum[N];
int main ()
{

	int n;
	euler(N);
	phi[1] = 0;
	for (int i = 1; i <= N; i++) sum[i] = sum[i - 1] + phi[i];
	while (scanf("%d", &n) != EOF && n)
	{
		printf("%lld\n", sum[n]);
	}
	return 0 ;
}
  • 方法
    欧拉函数线性筛
  • 总结
    把题目描述抽象处理,变成已知模型能力要加强。

猜你喜欢

转载自blog.csdn.net/strategist_614/article/details/88615005