题目链接
题意:从区间[1,B]选择一个x,从区间[1,D]选择一个数y,使得gcd(x,y)==k 的方案数。
莫比乌斯反演学习资料
1,文库
2,百科
莫比乌斯反演学习
纯属个人理解,错误之处还望指正。莫比乌斯反演就是两个函数之间的关系,感觉类似函数中的反函数,莫比乌斯反演是数论中一个比较重要的公式,用于简化加速计算。
公式是:
F(n)=∑d|nf(d)⇔f(n)=∑d|nu(d)F(nd)
或者另外一个形式
F(n)=∑n|df(d)⇔f(n)=∑n|du(dn)F(d)
其中函数μ(d),是与n无关的,定义如下:
(1)若d=1,那么u(d)=1
(2)若
d=p1p2p3….pk
,
pi
均为互异素数,那么
μ(d)=(−1)k
(3)其它情况下
μ(d)=0
μ(d)函数还有其他的性质,资料介绍的比较多,我就不写了。
d|n:表示n%d==0
然后基本的就是这样了,下面就是怎么使用了,使用的思想类似于数学上“正难则反”的逆向思维,比如我们要计算一个函数f(n)但是这个函数很难计算或者是比较耗费时间,那么我们可以尝试找一个比较好计算的函数Ϝ(n),满足
F(n)=∑d|nf(d)
么我们就可以转化为
f(n)=∑d|nu(d)F(nd)
这样我们就能够计算出f(n)了,因为u(n)函数我们可以使用欧拉函数O(n)时间预处理出来。
题目是计算[1,B],[1,D]前一个区间选择一个x,后一个区间选择一个y,满足gcd(x,y)==k,这样的(x,y)有几个,(1,2)(2,1)算作是一对。
gcd(x,y)==k也就等价于计算
不妨设f(k)表示gcd(x,y)==k的个数,也就是我们要计算的。
在设F(k):表示gcd(x,y)==(k的倍数)的个数。其中(1<=x<=B,1<=y<=D)
下面的B,D都是B=B/k;D=D/K;
于是 F(k)=∑k|df(d) ,同时由F函数的定义可得
代码:
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
#define cl(a,b) memset(a,b,sizeof(a))
#define ll long long
#define pb push_back
#define gcd __gcd
const double EPS = 1e-8;
const int maxn = 1e5+1000;
const int inf = 0x3f3f3f3f;
const double PI = acos(-1.0);
bool check[maxn];
int prime[maxn];
int mu[maxn];
void Moblus(){
cl(check,false);
mu[1]=1;int tot=0;
for(int i=2;i<maxn;i++){
if(!check[i]){
prime[tot++]=i;
mu[i]=-1;
}
for(int j=0;j<tot;j++){
if(i*prime[j]>maxn)break;
check[i*prime[j]]=true;
if(i%prime[j]==0){
mu[i*prime[j]]=0;break;
}
else mu[i*prime[j]]=-mu[i];
}
}
}
int main(){
int T;scanf("%d",&T);
Moblus();
int cas=1;
while(T--){
int a,b,c,d,k;
scanf("%d%d%d%d%d",&a,&b,&c,&d,&k);
if(k==0){
printf("Case %d: 0\n",cas++);continue;
}
ll ans1=0;
b/=k;
d/=k;
if(b>d)swap(b,d);
for(int i=1;i<=b;i++)ans1+=(ll)mu[i]*(b/i)*(d/i);
ll ans2=0;
for(int i=1;i<=b;i++)ans2+=(ll)mu[i]*(b/i)*(b/i);
ans2>>=1;
printf("Case %d: %lld\n",cas++,ans1-ans2);
}
return 0;
}