hdu 1695莫比乌斯反演

本文介绍了一种使用莫比乌斯反演解决特定数学问题的方法,即计算指定区间内满足最大公约数等于给定值的整数对数量。通过转换视角,将问题简化为计算相对质数对的数量,并利用μ函数实现快速计算。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

题意:
求区间[1,a],区间[1,b]中有多少对数满足gcd(i,j)=k,(2,2)只算一对 求区间[1,a] ,区间[1,b]中有多少对数满足gcd(i,j)=k,(2,2)只算一对 [1,a][1,b]gcd(i,j)=k,(2,2)

题解:莫比乌斯反演基础

f(x)f(x)f(x)gcd(i,j)=xgcd(i,j)=xgcd(i,j)=x的对数
F(x)F(x)F(x)gcd(i,j)=xgcd(i,j)=xgcd(i,j)=x的倍数的对数
b=ixb = \frac{i}{x}b=xi , d=jxd = \frac{j}{x}d=xj
那么我们的答案就是求f(x)f(x)f(x),可以转化为求gcd(b,d)=1gcd(b,d)=1gcd(b,d)=1的对数

由莫比乌斯反演公式
F(n)=∑d∣nf(d)F(n) = \displaystyle \sum_{d \mid n}f(d)F(n)=dnf(d)

f(n)=∑d∣nμ(d)F(nd)f(n) = \displaystyle \sum_{d \mid n} μ(d)F(\frac{n}{d})f(n)=dnμ(d)F(dn)

F(x)F(x)F(x)其实就是(bi)∗(di)(\frac{b}{i})*(\frac{d}{i})(ib)(id)

那么f(x)=∑i=1min(b,d)μ(i)(bi)∗(di))f(x)= \displaystyle \sum^{min(b,d)}_{i=1} μ(i)(\frac{b}{i})*(\frac{d}{i}))f(x)=i=1min(b,d)μ(i)(ib)(id))

μ函数可以在线性的时间内筛出来,最后答案累加就好了,注意这样计算,如果两个数相同那么就会被计算两次,所以要减去重复的\mu函数可以在线性的时间内筛出来,最后答案累加就好了,注意这样计算,如果两个数相同那么就会被计算两次,所以要减去重复的μ线

重复的其实只有m=min(b,d)m = min(b,d)m=min(b,d)中的情况,画个图就知道了

/* ***********************************************
Author        :pall_scall
Created Time  :2019年07月17日 星期三 16时58分45秒
File Name     :acm.cpp
************************************************ */

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <queue>
#include <vector>
#include <map>
#include <set>
#include <stack>
#include <string>
using namespace std;

typedef long long ll;
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define mem(a,b) memset(a,b,sizeof(a))
#define lowbit(x) x&-x
const int MAXN = 1e6 + 5;
const int mod = 1e9 + 7;
const double eps = 1e-6;
const double pi = acos(-1.0);

bool check[MAXN+10];
ll prime[MAXN+10];
int mu[MAXN+10];
void Moblus(){
    memset(check,false,sizeof(check));
    mu[1] = 1;
    ll tot = 0;
    for(ll i = 2; i <= MAXN; i++){
        if( !check[i] ){
            prime[tot++] = i;
            mu[i] = -1;
        }
        for(ll 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 a,b,c,d,k;
int main(){
	//freopen("in.txt","r",stdin);
	//freopen("out.txt","w",stdout);
	Moblus();
	int t;
	cin>>t;
	int ca = 0;	
	while(t--){
		cin>>a>>b>>c>>d>>k;
		if(k == 0){
			printf("Case %d: 0\n",++ca);
			continue;
		}
		b = b/k, d = d/k;
		if(b > d) swap(b,d);
		ll ans = 0,res = 0;
		for(int i = 1; i <= b; i++){
			ans += (ll)mu[i]*(b/i)*(b/i);
			res += (ll)mu[i]*(b/i)*(d/i);
		}
		printf("Case %d: %lld\n",++ca,res-ans/2);
	}
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值