HDU 1695 GCD 容斥原理/莫比乌斯反演

2024-06-12 11:38

本文主要是介绍HDU 1695 GCD 容斥原理/莫比乌斯反演,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

题意:

给你两个区间[a,b],[c,d],还有一个k。让你从区间[a,b]中找出x,[c,d]中找出y,问共有多少组(x,y)使得gcd(x,y)=k。

(x,y)和(y,x)算一组。

思路:

参考:http://blog.csdn.net/yang_7_46/article/details/9072533

容斥。


普通容斥:

*如果gcd(x,y)=k,则gcd(x/k,y/k)=1。那么对于两个区间来说,我们都默认进行了b = b/k,d = d/k操作。(a,c固定为1)。(这样做并不会使得组数遗漏)

因为(x,y)和(y,x)只算一组,因此在容斥过程中,我们要保证x < y,这样才不会重复计算。

剩下就是容斥的过程,可以戳这里。

其它一些细节要自己处理一下,例如b = 1或者d = 1的时候。


莫比乌斯反演:

莫比乌斯资料:http://blog.csdn.net/acdreamers/article/details/8542292

定义f(n):gcd(x, y)为n的方案数。

定义F(n):gcd(x, y)是n的倍数的方案数。

则我们要求的就是f(1)。默认b < d,并且都已经除以k。

F(n) = (b/n)*(d/n);

套用莫比乌斯公式即可。


code(普通容斥):

#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <iostream>
#include <algorithm>
#include <vector>
using namespace std;const int MAXN = 1e5+5;
typedef long long LL;int a, b, c, d, k;
bool isp[MAXN+5];
vector <int> vec[MAXN+5];void prime()
{memset(isp, false, sizeof(isp));for(int i = 2;i <= MAXN; i++){if(!isp[i]){vec[i].push_back(i);for(int j = i*2; j <= MAXN; j += i){isp[j] = true;vec[j].push_back(i);}}}
}
LL calc(int t, int p)
{LL ret, v = 1;int cnt = 0;for(int i = 0;i < vec[t].size(); i++){if((1<<i)&p){cnt++;v *= vec[t][i];}}if(v == 0) return 0;ret = (LL)(d-t)/v;if(cnt%2 == 0) ret = -ret;return ret;
}void solve()
{//tepan 1LL res = d;for(int i = 2;i <= b; i++){res += d-i;for(int j = 1;j < (1<<vec[i].size()); j++)res -= calc(i, j);}printf("%I64d\n", res);
}int main()
{prime();int T, cas = 0;scanf("%d", &T);while(T--){scanf("%d%d%d%d%d", &a, &b, &c, &d, &k);printf("Case %d: ", ++cas);if(k == 0){puts("0");continue;}b /= k, d /= k;if(b > d) swap(b, d);if(b == 0 || d == 0)puts("0");else solve();}return 0;
}

 code(莫比乌斯反演):

#include <bits/stdc++.h>
using namespace std;const int N = 1e5+5;
typedef long long LL;int b, d, k;
int prime[N], cnt;
int mu[N];
bool vis[N];
void Mobius() {mu[1] = 1;cnt = 0;for(int i = 2;i < N; i++) {if(!vis[i]) {mu[i] = -1;prime[cnt++] = i;}for(int j = 0;j < cnt; j++) {if(i*prime[j] >= N) break;vis[i*prime[j]] = true;if(i%prime[j] != 0)mu[i*prime[j]] = -mu[i];else {mu[i*prime[j]] = 0;break;}}}
}int main() {Mobius();int T, cas = 0;scanf("%d", &T);while(T--) {scanf("%*d%d%*d%d%d", &b, &d, &k);printf("Case %d: ", ++cas);if(k == 0) {puts("0");continue;}b /= k, d /= k;if(b > d) swap(b, d);LL ans = 0;for(int i = 1;i <= b; i++)ans += 1ll*mu[i]*(b/i)*(d/i);LL tmp = 0;for(int i = 1;i <= b; i++)tmp += 1ll*mu[i]*(b/i)*(b/i);printf("%I64d\n", ans-tmp/2);}return 0;
}



这篇关于HDU 1695 GCD 容斥原理/莫比乌斯反演的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1054115

相关文章

Java编译生成多个.class文件的原理和作用

《Java编译生成多个.class文件的原理和作用》作为一名经验丰富的开发者,在Java项目中执行编译后,可能会发现一个.java源文件有时会产生多个.class文件,从技术实现层面详细剖析这一现象... 目录一、内部类机制与.class文件生成成员内部类(常规内部类)局部内部类(方法内部类)匿名内部类二、

Python中随机休眠技术原理与应用详解

《Python中随机休眠技术原理与应用详解》在编程中,让程序暂停执行特定时间是常见需求,当需要引入不确定性时,随机休眠就成为关键技巧,下面我们就来看看Python中随机休眠技术的具体实现与应用吧... 目录引言一、实现原理与基础方法1.1 核心函数解析1.2 基础实现模板1.3 整数版实现二、典型应用场景2

Java的IO模型、Netty原理解析

《Java的IO模型、Netty原理解析》Java的I/O是以流的方式进行数据输入输出的,Java的类库涉及很多领域的IO内容:标准的输入输出,文件的操作、网络上的数据传输流、字符串流、对象流等,这篇... 目录1.什么是IO2.同步与异步、阻塞与非阻塞3.三种IO模型BIO(blocking I/O)NI

JAVA封装多线程实现的方式及原理

《JAVA封装多线程实现的方式及原理》:本文主要介绍Java中封装多线程的原理和常见方式,通过封装可以简化多线程的使用,提高安全性,并增强代码的可维护性和可扩展性,需要的朋友可以参考下... 目录前言一、封装的目标二、常见的封装方式及原理总结前言在 Java 中,封装多线程的原理主要围绕着将多线程相关的操

kotlin中的模块化结构组件及工作原理

《kotlin中的模块化结构组件及工作原理》本文介绍了Kotlin中模块化结构组件,包括ViewModel、LiveData、Room和Navigation的工作原理和基础使用,本文通过实例代码给大家... 目录ViewModel 工作原理LiveData 工作原理Room 工作原理Navigation 工

Java的volatile和sychronized底层实现原理解析

《Java的volatile和sychronized底层实现原理解析》文章详细介绍了Java中的synchronized和volatile关键字的底层实现原理,包括字节码层面、JVM层面的实现细节,以... 目录1. 概览2. Synchronized2.1 字节码层面2.2 JVM层面2.2.1 ente

MySQL的隐式锁(Implicit Lock)原理实现

《MySQL的隐式锁(ImplicitLock)原理实现》MySQL的InnoDB存储引擎中隐式锁是一种自动管理的锁,用于保证事务在行级别操作时的数据一致性和安全性,本文主要介绍了MySQL的隐式锁... 目录1. 背景:什么是隐式锁?2. 隐式锁的工作原理3. 隐式锁的类型4. 隐式锁的实现与源代码分析4

MySQL中Next-Key Lock底层原理实现

《MySQL中Next-KeyLock底层原理实现》Next-KeyLock是MySQLInnoDB存储引擎中的一种锁机制,结合记录锁和间隙锁,用于高效并发控制并避免幻读,本文主要介绍了MySQL中... 目录一、Next-Key Lock 的定义与作用二、底层原理三、源代码解析四、总结Next-Key L

Spring Cloud Hystrix原理与注意事项小结

《SpringCloudHystrix原理与注意事项小结》本文介绍了Hystrix的基本概念、工作原理以及其在实际开发中的应用方式,通过对Hystrix的深入学习,开发者可以在分布式系统中实现精细... 目录一、Spring Cloud Hystrix概述和设计目标(一)Spring Cloud Hystr

MySQL中的MVCC底层原理解读

《MySQL中的MVCC底层原理解读》本文详细介绍了MySQL中的多版本并发控制(MVCC)机制,包括版本链、ReadView以及在不同事务隔离级别下MVCC的工作原理,通过一个具体的示例演示了在可重... 目录简介ReadView版本链演示过程总结简介MVCC(Multi-Version Concurr