hdu1695 容斥原理 莫比乌斯反演
给定两个数b,d,问[1,b]和[1,d]区间上有多少对互质的数。(x,y)和(y,x)算一个。
对于[1,b]部分,用欧拉函数直接求。对于大于b的部分,求n在[1,b]上有多少个互质的数,用容斥原理。
主要学习容斥原理的写法,本题使用DFS。容斥原理复杂度比较高,是指数复杂度。
输出长整型不能用lld,必须用I64d。
#include<stdio.h>
#include<iostream>
#include<stdlib.h>
#include<string.h>
#include<algorithm>
using namespace std;
typedef long long ll;
const int maxn = 1e5 + 7;
int a, b, c, d, k;
bool isPrime[maxn];
ll euler[maxn];
ll sumeuler[maxn];
int p[maxn][20], ps[maxn];
void initPrimes(){
memset(isPrime, 1, sizeof(isPrime));
for (int i = 1; i < maxn; i++){
euler[i] = i;
}
for (int i = 2; i < maxn; i++){
if (isPrime[i]){
for (int j = i; j < maxn; j += i){
isPrime[j] = false;
p[j][ps[j]++] = i;
euler[j] = euler[j] * (i - 1) / i;
}
}
}
sumeuler[1] = 1;
for (int i = 2; i < maxn; i++){
sumeuler[i] = sumeuler[i - 1] + euler[i];
}
}
ll dfs(int ind, int n, int x){//容斥原理核心代码
ll s = 0;
for (int i = ind; i < ps[x]; i++){
s += n / p[x][i] - dfs(i + 1, n / p[x][i], x);
}
return s;
}
ll huzhi(int n, int x){//0~n之间,与x互质的数字个数
return n - dfs(0, n, x);
}
int main(){
freopen("in.txt", "r", stdin);
initPrimes();
int caseCount;
scanf("%d", &caseCount);
for (int caseid = 1; caseid <= caseCount; caseid++){
scanf("%d%d%d%d%d", &a, &b, &c, &d, &k);
ll ans;
if (k == 0)
ans = 0;
else{
if (b > d)swap(b, d);
b /= k, d /= k;
ans = sumeuler[b];
for (int i = b + 1; i <= d; i++){
ans += huzhi(b, i);
}
}
printf("Case %d: %I64d\n", caseid, ans);
}
return 0;
}
容斥原理的另一种写法:
int calc(int n,int m)//n < m,求1-n内和m互质的数的个数
{
getFactors(m);
int ans = 0;
for(int i = 1;i < (1<<fatCnt);i++)
{
int cnt = 0;
int tmp = 1;
for(int j = 0;j < fatCnt;j++)
if(i&(1<<j))
{
cnt++;
tmp *= factor[j][0];
}
if(cnt&1)ans += n/tmp;
else ans -= n/tmp;
}
return n - ans;
}
容斥原理项的个数为2的幂次,肯定不会太大,所以一定可以用一个int来表示所有情况。
本题还可以用莫比乌斯反演来解决。
#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>
#include <vector>
#include <queue>
#include <set>
#include <map>
#include <string>
#include <math.h>
#include <stdlib.h>
#include <time.h>
using namespace std;
const int MAXN = 1e5 + 7;
typedef long long ll;
bool isprime[MAXN];
int prime[MAXN], psize;
int mu[MAXN];
ll euler[MAXN], sumeuler[MAXN];
void init(){
memset(isprime, false, sizeof(isprime));
mu[1] = euler[1] = sumeuler[1] = 1;
psize = 0;
for (int i = 2; i <= MAXN; i++) {
if (!isprime[i]) {
prime[psize++] = i;
mu[i] = -1;
euler[i] = i - 1;
}
for (int j = 0; j < psize&&i*prime[j] < MAXN; j++) {
isprime[i * prime[j]] = true;
if (i % prime[j] == 0) {
mu[i * prime[j]] = 0;
euler[i*prime[j]] = euler[i] * prime[j];
break;
}
else {
mu[i * prime[j]] = -mu[i];
euler[i*prime[j]] = euler[i] * (prime[j] - 1);
}
}
sumeuler[i] = sumeuler[i - 1] + euler[i];
}
}
ll eu(int b){
long long ans2 = 0;
for (int i = 1; i <= b; i++)
ans2 += (long long)mu[i] * (b / i)*(b / i);
return ans2;
}
void debug(){
for (int i = 1; i < MAXN; i++){
if (eu(i) / 2 != sumeuler[i] - 1)
cout <<i<<" "<< eu(i) / 2 << " " << sumeuler[i] - 1 << endl;
}
exit(0);
}
int main(){
freopen("in.txt", "r", stdin);
int T;
int a, b, c, d, k;
init();
//debug();
scanf("%d", &T);
int iCase = 0;
while (T--)
{
iCase++;
scanf("%d%d%d%d%d", &a, &b, &c, &d, &k);
if (k == 0) {
printf("Case %d: 0\n", iCase);
continue;
}
b /= k; d /= k;
if (b > d)swap(b, d);
long long ans1 = 0;
for (int i = 1; i <= b; i++)
ans1 += (long long)mu[i] * (b / i)*(d / i);
//这里ans2表示重复的部分,这部分可以用欧拉函数直接求,完全不需要for循环,但是提交却是wa,经过本机检测,完全没有问题
long long ans2 = 0;
for (int i = 1; i <= b; i++)
ans2 += (long long)mu[i] * (b / i)*(b / i);
ans1 -= ans2 / 2;
//ans1 -= (sumeuler[b] - 1);
printf("Case %d: %I64d\n", iCase, ans1);
}
return 0;
}
hdu1695 容斥原理 莫比乌斯反演的更多相关文章
- HDU1695 GCD(莫比乌斯反演)
传送门 看了1个多小时,终于懂了一点了 题目大意:给n,m,k.求gcd(x,y) = k(1<=x<=n, 1<=y<=m)的个数 思路:令F(i)表示i|gcd(x,y)的 ...
- hdu1695(莫比乌斯反演)
传送门:GCD 题意:求[1,n],[1,m]gcd为k的对数. 分析:莫比乌斯入反演门题,gcd(x,y)==k等价于gcd(x/k,y/k)==1,求出[1,n][1,m]互质的对数,在减去[1, ...
- BZOJ 2440 中山市选2011 全然平方数 二分答案+容斥原理+莫比乌斯反演
题目大意:求第k个无平方因子数是多少(无视原题干.1也是全然平方数那岂不是一个数也送不出去了? 无平方因子数(square-free number),即质因数分解之后全部质因数的次数都为1的数 首先二 ...
- hdu1695 GCD 莫比乌斯反演做法+枚举除法的取值 (5,7),(7,5)看做同一对
/** 题目:hdu1695 GCD 链接:http://acm.hdu.edu.cn/status.php 题意:对于给出的 n 个询问,每次求有多少个数对 (x,y) , 满足 a ≤ x ≤ b ...
- BZOJ2005: [Noi2010]能量采集(容斥原理 莫比乌斯反演)
Time Limit: 10 Sec Memory Limit: 512 MBSubmit: 4727 Solved: 2877[Submit][Status][Discuss] Descript ...
- [HAOI2011][bzoj2301] Problem b [莫比乌斯反演+容斥原理+分块前缀和优化]
题面: 传送门 有洛谷就尽量放洛谷链接呗,界面友好一点 思路: 和HDU1695比较像,但是这一回有50000组数据,直接莫比乌斯反演慢慢加的话会T 先解决一个前置问题:怎么处理a,c不是1的情况? ...
- BZOJ2301: [HAOI2011]Problem b[莫比乌斯反演 容斥原理]【学习笔记】
2301: [HAOI2011]Problem b Time Limit: 50 Sec Memory Limit: 256 MBSubmit: 4032 Solved: 1817[Submit] ...
- 【BZOJ-2440】完全平方数 容斥原理 + 线性筛莫比乌斯反演函数 + 二分判定
2440: [中山市选2011]完全平方数 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 2371 Solved: 1143[Submit][Sta ...
- hdu1695 莫比乌斯反演
莫比乌斯反演:可参考论文:<POI XIV Stage.1 <Queries>解题报告By Kwc-Oliver> 求莫比乌斯函数mu[i]:(kuangbin模板) http ...
随机推荐
- RecyclerView 缓存机制详解
一 前言 RecyclerView据官方的介绍,该控件用于在有限的窗口中展示大量数据集,其实这样功能的控件我们并不陌生,例如:ListView.GridView.RecyclerView可以用来代替传 ...
- SHA加密算法工具类
1.利用Apache的工具类实现加密: maven: <dependency> <groupId>commons-codec</groupId> <artif ...
- 如何获取隔壁wifi密码,非暴力破解
目前常见的Wi-Fi加密方式有WEP.WPA2和WPS(链接为各自的破解方式),不过有网友反映以往破解WPA2的方法耗时太长,而且不适用于所有WPS启动的接入点.而今天介绍的这种方法则更加省时省力. ...
- tyvj P2018 「Nescafé26」小猫爬山 解题报告
P2018 「Nescafé26」小猫爬山 时间: 1000ms / 空间: 131072KiB / Java类名: Main 背景 Freda和rainbow饲养了N只小猫,这天,小猫们要去爬山.经 ...
- 【BZOJ】【1419】Red is good
数学期望/期望DP 还是戳<浅析竞赛中一类数学期望问题的解决方法>这篇论文…… $$ f[i][j]= \begin{cases} 0 &, &i==0 \\ f[i-1] ...
- wiggle-subsequence
// 参考了:https://discuss.leetcode.com/topic/51893/two-solutions-one-is-dp-the-other-is-greedy-8-lines ...
- smtp发送带附件的邮件(直接将string类型结果保存为附件)
该方式直接保存为HTML文件,也可以是文本文件,其它格式效果不是很好 MailMessage mmsg = new MailMessage(); mmsg.Subject = " ...
- Cognos事件工作室Event Studio开发步骤
Cognos本身是很强大的,只是很多人的思维只是局限在数据-模型-展示的层面上,下面我们就来介绍一下Cognos中稍微有些冷门的一个组件事件工作室(IBM Cognos Event Studio),概 ...
- 动态改变UIPopupList选项(NGUI)
NGUI的UIPopupList 可以通过修改items属性来动态改变菜单选项: public class popListvahnge : MonoBehaviour { public UIPopup ...
- sqoop安装部署(笔记)
sqoop是一个把关系型数据库数据抽向hadoop的工具.同时,也支持将hive.pig等查询的结果导入关系型数据库中存储.由于,笔者部署的hadoop版本是2.2.0,所以sqoop的版本是:sqo ...