【2015广东工业大学新生赛D】【数论-欧拉函数】求互质对数 1~n中互质数对数

最新推荐文章于 2021-07-31 11:00:09 发布

snowy_smile

最新推荐文章于 2021-07-31 11:00:09 发布

阅读量1.6k

点赞数

分类专栏：数论-欧拉函数文章标签：数论-欧拉函数

本文链接：https://blog.csdn.net/snowy_smile/article/details/50209391

版权

数论-欧拉函数专栏收录该内容

3 篇文章 0 订阅

订阅专栏

#include<stdio.h>
#include<iostream>
#include<string.h>
#include<string>
#include<ctype.h>
#include<math.h>
#include<set>
#include<map>
#include<vector>
#include<queue>
#include<bitset>
#include<algorithm>
#include<time.h>
using namespace std;
void fre(){freopen("c://test//input.in","r",stdin);freopen("c://test//output.out","w",stdout);}
#define MS(x,y) memset(x,y,sizeof(x))
#define MC(x,y) memcpy(x,y,sizeof(x))
#define MP(x,y) make_pair(x,y)
#define ls o<<1
#define rs o<<1|1
typedef long long LL;
typedef unsigned long long UL;
typedef unsigned int UI;
template <class T1,class T2>inline void gmax(T1 &a,T2 b){if(b>a)a=b;}
template <class T1,class T2>inline void gmin(T1 &a,T2 b){if(b<a)a=b;}
const int N=0,M=0,Z=1e9+7,ms63=1061109567;
int casenum,casei;
int phi[1010];
void EULERinit()
{
	int top=1000;
	for(int i=2;i<=top;++i)phi[i]=i;
	for(int i=2;i<=top;++i)
	{
		if(phi[i]==i)for(int j=i;j<=top;j+=i)phi[j]-=phi[j]/i;
		phi[i]+=phi[i-1];
	}
}
int main()
{
	EULERinit();
	scanf("%d",&casenum);
	for(casei=1;casei<=casenum;++casei)
	{
		int n;scanf("%d",&n);
		printf("%d\n",phi[n]);
	}
	return 0;
}
/*
【题意】
问你，有多少对（x，y），1<=x<y<=n，满足gcd（x，y）==1

【类型】
欧拉函数

【分析】
我们发现，这个就是忽略（1，1）后，欧拉函数的前n项和。
于是暴力打个表就可以AC啦！

【时间复杂度&&优化】
O（nlogn）

*/