1.题目描述:
2.题意:
问A到B的范围内,与N互质的数字有多少个。
3.思路:
唯一分解定理+容斥原理。 这题我们先来看看互质的条件:gcd(a,b)==1。
1)那么我们可以先把N分解,然后划去A到B范围内,N的质因数的倍数。 那么问题就转化为,求A到B范围内,N的质因数的倍数的多少(用num表示),然后最后的答案就是:B-A+1-num。
2)那么我们下一个问题就是,如何求倍数的数量了。由于我们直接求A到B的数量,边界比较麻烦。所以我们转换为:1到B的以及1到A-1的数量,然后相减就好啦!
3)那么对于[1,n]数量内的质因数p的倍数有多少呢,有:[n/p] (就是n整除p)个。
4)但是我们得考虑一个问题,一个数既是p1的倍数,也是p2的倍数,那么它会被计算两次,这样的话我们就会减多了,那么我们就要加回一次。那么对于一个由k个质因数组成的数字d,那么我们要加还是减呢?我们来看个例子:
假设n=12,N是30,那么质因数就是2,3,5。
2的倍数:2,4,6,8,10,12
3的倍数:3,6,9,12
5的倍数:5,10
我们会发现:6,10,12出现了两次,他们的共同特点是:他们是两个集合的交集。我们现在要求的是三个集合的并集,那么我们用容斥原理 就有:集合的并集的数量=三个集合的数量,减去他们两两之间集合的交集的数量,加上三三之间集合的交集的数量,减去四四集合之间的交集的数量… 我们发现,是偶数个集合的交就减,奇数个集合的交就加。
对于我们的倍数集合也是如此,2和3的倍数的交集合,就是2*3的倍数的集合。我们用0,1来代表选与不选,这样就可以拼凑出所有的可能选项。因为不可能一个集合都不选,所以我们从1开始,到(1<<k),就是k个质因数组成的所有状态。
4.代码:
//E
//#include<bits/stdc++.h>
//#pragma GCC optimize(3,"Ofast","inline")
#include<iostream>
#include<string>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<vector>
//#include<random>
#include<cstdlib>
#include<ctime>
#include<map>
#include<set>
#include<stack>
#include<queue>
#define FAST ios::sync_with_stdio(false)
#define DEV_RND ((int)rand()*RAND_MAX+rand())
#define RND(L,R) (DEV_RND%((R)-(L)+1)+(L))
#define abs(a) ((a)>=0?(a):-(a))
#define sz(x) ((int)(x).size())
#define all(x) (x).begin(),(x).end()
#define mem(a,b) memset(a,b,sizeof(a))
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
#define rep(i,a,n) for(int i=a;i<n;++i)
#define per(i,n,a) for(int i=n-1;i>=a;--i)
#define pb push_back
#define mp make_pair
#define fi first
#define se second
#define li inline
#define re register
using namespace std;
//typedef uniform_int_distribution<int> RNDI;
typedef pair<int,int> PII;
typedef vector<int> VI;
typedef double db;
typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
const int Hash = 131;//13331
const int maxn = 1e5+5;
const int maxm = 100000+5;
const int mod = 1e9+7;
const int inf = 0x3f3f3f3f;
const double eps = 1e-7;
const double pi = acos(-1);
//int dir[4][2]={{-1,0},{1,0},{0,-1},{0,1}};
//li int f(int x){return x==par[x]?par[x]:par[x]=f(par[x]);}
//mt19937 eng(time(0));
//li int RND(int L,int R){RNDI rnd(L,R);return rnd(eng);}
li ll lowbit(ll x){return x&(-x);}
ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
void ex_gcd(ll a,ll b,ll &d,ll &x,ll &y){if(!b){d=a,x=1,y=0;}else{ex_gcd(b,a%b,d,y,x);y-=x*(a/b);}}//x=(x%(b/d)+(b/d))%(b/d);
li ll qpow(ll a,ll b,ll MOD=mod){ll res=1;a%=MOD;while(b>0){if(b&1)res=res*a%MOD;a=a*a%MOD;b>>=1;}return res%MOD;}
li ll qmul(ll a,ll b,ll MOD=mod){return (a*b-(ll)((long double)a/MOD*b)*MOD+MOD)%MOD;}
li ll inv(ll x,ll p){return qpow(x,p-2,p);}
li ll jos(ll n,ll k,ll s=1){ll res=0;rep(i,1,n+1) res=(res+k)%i;return (res+s)%n;}
db f(db x){return x;}
li db sim(db l,db r){return (f(l)+4.*f((l+r)/2.)+f(r))*(r-l)/6.;}
db asr(db l,db r,db ans,db eps){db m=l+(r-l)/2.,L=sim(l,m),R=sim(m,r);return fabs(L+R-ans)<=15.*eps?L+R+(L+R-ans)/15.:asr(l,m,L,eps/2)+asr(m,r,R,eps/2);}
db asr(db l,db r,db eps){return asr(l,r,sim(l,r),eps);}
namespace IO
{
li ll read()
{
ll x=0,sign=1;char c=getchar();
while(c>'9'||c<'0') {if(c=='-') sign=-1;c=getchar();}
while('0'<=c&&c<='9') x=x*10+c-'0',c=getchar();
return x*sign;
}
template<typename T>
li void write(T x,char t='\n')
{
if(x<0){x=-x;putchar('-');};
static int sta[25];int top=0;
do{sta[top++]=x%10;}while(x/=10);
while(top) putchar(sta[--top]+'0');
putchar(t);
}
}
using namespace IO;
/*-------------head-------------*/
//
ll l,r,n;
VI p;
li ll solve(ll x,ll n)
{
p.clear();
for(int i=2;i*i<=n;++i)
if(!(n%i))
{
p.pb(i);
while(!(n%i)) n/=i;
}
if(n>1) p.pb(n);
ll ans=0;
rep(pr,1,(1<<sz(p)))
{
ll k=0,d=1;
rep(i,0,sz(p))
{
if(pr&(1<<i))
{
k++;
d*=p[i];
}
}
if(k&1) ans+=x/d;
else ans-=x/d;
}
return x-ans;
//puts("");
}
int main()
{
//srand(time(0));
//freopen("C:\\Users\\Administrator\\Desktop\\in.txt","r",stdin);
int QAQ=1;
for(int QwQ=read();QwQ;QwQ--)
{
l=read(),r=read(),n=read();
printf("Case #%d: %lld\n",QAQ++,solve(r,n)-solve(l-1,n));
}
//while(~scanf("%d",&n)) solve();
return 0;
}