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;
}
资源下载链接为: https://pan.quark.cn/s/9e7ef05254f8 在苹果的生态系统中,IAP(应用内购买)是苹果应用商店(App Store)中应用开发者常采用的一种盈利模式,允许用户在应用内直接购买虚拟商品或服务。苹果为开发者提供了一份详细的人民币(CNY)IAP定价表,这份定价表具有以下特点: 价格分级:定价表由多个价格等级组成,开发者可根据虚拟商品的价值选择相应等级,等级越高,价格越高。例如,低等级可能对应基础功能解锁,高等级则对应高级服务或大量虚拟道具。 税收与分成:苹果会从应用内购买金额中抽取30%作为服务费或佣金,这是苹果生态的固定规则。不过,开发者实际到手的收入会因不同国家和地区的税收政策而有所变化,但定价表中的价格等级本身是固定的,便于开发者统一管理。 多级定价策略:通过设置不同价格等级,开发者可以根据商品或服务的类型与价值进行合理定价,以满足不同消费能力的用户需求,从而最大化应用的总收入。例如,一款游戏可以通过设置不同等级的虚拟货币包,吸引不同付费意愿的玩家。 特殊等级:除了标准等级外,定价表还包含备用等级和特殊等级(如备用等级A、备用等级B等),这些等级可能是为应对特殊情况或促销活动而设置的额外价格点,为开发者提供了更灵活的定价选择。 苹果IAP定价表是开发者设计应用内购机制的重要参考。它不仅为开发者提供了标准的收入分成模型,还允许开发者根据产品特性设定价格等级,以适应市场和满足不同用户需求。同时,开发者在使用定价表时,还需严格遵守苹果的《App Store审查指南》,包括30%的分成政策、使用苹果支付接口、提供清晰的产品描述和定价信息等。苹果对应用内交易有严格规定,以确保交易的透明性和安全性。总之,苹果IAP定价表是开发者在应用内购设计中不可或缺的工具,但开发者也需密切关注苹果政策变化,以确保应用的合规运营和收益最大化。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值