GCD HDU - 1695

原题链接
考察:容斥原理
fw本f,二刷还是没想到.
思路
因为x,y与y,x是重复的,所以只计数x,y(y>=x)的个数.暴力枚举的话是O(n²),需要优化.一般优化第二重循环.也就是对于每一个x,找到与它gcd(x,y) = k的y.
gcd(x,y) = k —>gcd(x/k,y/k) = 1.这里比较容易想到欧拉函数,但是欧拉函数求不出来,原因在与y/k互质的数不一定会与i/k互质.那么如何找到与x/k互质的数？将x/k分解质因数,如果x/k~y/k之间有质因数的倍数,那么一定不互质,求个数的时候就是用容斥原理

#include <iostream>
#include <vector>
using namespace std;
typedef long long LL;
const int N = 100010;
bool st[N];
vector<int> v[N];
void init()
{
	for(int i=2;i<=N-10;i++)
	{
		if(st[i]) continue;
		for(int j=i;j<=N-10;j+=i)
		{
			v[j].push_back(i);
			st[j] = 1;
		}
	}
}
int main()
{
	int T,kcase = 0;
	scanf("%d",&T);
	init();
	while(T--)
	{
		int a,b,c,d,k;
		scanf("%d%d%d%d%d",&a,&b,&c,&d,&k);
		if(!k||k>b||k>d) {printf("Case %d: %lld\n",++kcase,0);continue;}
		if(b>d) swap(b,d);
		b/=k,d/=k;//找i~d范围内与i互质的数. 
		LL ans = d;//因为b>=1,所以当i=1时,互质的数为d. 
		for(int i=2;i<=b;i++)
		{
			LL temp = d-i; 
			for(int j=1;j<1<<v[i].size();j++)
			{
				LL res = 1; int cnt = 0;
				for(int p=0;p<v[i].size();p++)
			       if(j>>p&1)
			       {
			       	  if(res*v[i][p]>d) {res=-1;break;} 
			       	  res*=v[i][p]; cnt++;
				   }
				if(res==-1) continue;
				if(cnt&1) temp-=d/res-i/res;
				else temp+=d/res-i/res;//i~d之间互质的数 
			}
			ans+=temp;
		}
		printf("Case %d: %lld\n",++kcase,ans);
	}
	return 0;
}