题目大意 :给出两个范围a-b、c-d,以及一个数字k,求GCD(x,y) = k的组合数。
GCD()为最大公约数,x为范围a到b的整数,y为c到d的整数。
思路: 两边同时除以k,变成GCD(x/k,y/k) = 1,也就是在1-b/k,和1-d/k中找到gcd为1的数,也就是互质的数的组合。
这样可以分为两个阶段,(1, b),以及(b+1, d)中对于(1,b)互质的数。
分两段的原因是(1, b)之间会有相同组合的数,(1,2)(2,1),这一段中出了(1,1)只有一次之外都会变成两次的,所以要+1再除以2。
第二段就不会有重叠的问题。
然后就(1,b)中每一个数分解质因子,用容斥原理来求(1,b)中互质的数。
也就是减去都能被2整除的数,减去都能被3整除的数,加上能被6整除的数…
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cstdio>
using namespace std;
typedef long long ll;
const int maxn = 1e6+10;
int prime[maxn];
int len = 0;
void divi(int n)
{
len = 0;
for(int i = 2;i*i<=n;i++)
{
if(n%i == 0)
prime[len++] = i;
while(n%i == 0)
n/=i;
}
if(n >1)
prime[len++] = n;
}
int solve(int n,int b)
{
divi(n);
int ans = b;
for(int i = 1;i<(1<<len);i++)
{
int cnt = 0;
int temp = 1;
for(int j = 0;j<len;j++)
{
if(i&(1<<j))
{
cnt++;
temp *= prime[j];
}
}
if(cnt%2 == 0)
ans += b/temp;
else
ans -=b/temp;
}
return ans;
}
int main()
{
int t;
scanf("%d", &t);
int kase = 0;
while(t--)
{
int a,c,b,d, k;
scanf("%d%d%d%d%d", &a, &b, &c,&d,&k);
if(k == 0)
{
printf("Case %d: 0\n", ++kase);
continue;
}
b=b/k;
d = d/k;
ll ans = 0;
if(b>d)
swap(b,d);
for(int i = 1;i<=b;i++)
{
ans += solve(i, b);
}
ans = (ans+1)/2;
//除以2的时候多减去了一个(1,1)
for(int i = b+1;i<=d;i++)
{
ans += solve(i,b);
}
printf("Case %d: %lld\n", ++kase, ans);
}
return 0;
}