Bzoj2301: [HAOI2011]Problem b:莫比乌斯反演
2016-04-19 09:10
453 查看
题目链接:[HAOI2011]Problem b
妈妈我终于会反演啦【雾
其实我就是来码个公式复习一下QwQ
前置技能:
![](https://oscdn.geek-share.com/Uploads/Images/Content/202007/03/76ad09be3790c6157af934dd69723b86)
那么按照题意,我们令
![](https://oscdn.geek-share.com/Uploads/Images/Content/202007/03/fa853a1d88b6797f54f3b1966deb34de)
首先明确我们为什么要反演,因为我们要求一个不好直接算的东西,而他和一个容易计算的东西有关系,所以我们要同过那个容易计算的东西来反向计算出我们真正要求的东西
于是我们令
![](https://oscdn.geek-share.com/Uploads/Images/Content/202007/03/d7bab0e62a48120d3d5832c10927be00)
,其中?是整数,整个式子表示k的倍数
根据前置技能得
![](https://oscdn.geek-share.com/Uploads/Images/Content/202007/03/e3d88446dc8c7c4a058d40d4ac5d7cd4)
这时候发现d/k求起来很麻烦,换个元,令d'=d/k
于是
![](https://oscdn.geek-share.com/Uploads/Images/Content/202007/03/b43f3211153857bd103b27b8d6fcdfef)
再来化简f函数得到
![](https://oscdn.geek-share.com/Uploads/Images/Content/202007/03/16a45d6a9d362e90c1faefe0f3ee8a62)
所以
![](https://oscdn.geek-share.com/Uploads/Images/Content/202007/03/3a9c8bf69b780b54c95b4433126306a7)
所以
![](https://oscdn.geek-share.com/Uploads/Images/Content/202007/03/2d9c7408a6d7a55bc0140b714e45e819)
然后分块求就好了,线性筛出莫比乌斯函数并求出前缀和,注意容斥一下
#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define LL long long
using namespace std;
const int maxn=100010;
int a,b,c,d,k;
int p[maxn],N=0,mu[maxn];
bool vis[maxn];
void get_pre(){
mu[1]=1;
for (int i=2;i<=maxn-5;++i){
if (!vis[i]){
p[++N]=i;
mu[i]=-1;
}
for (int j=1;j<=N&&p[j]*i<=maxn-5;++j){
vis[p[j]*i]=1;
if (i%p[j]==0){mu[p[j]*i]=0;break;}
mu[p[j]*i]=-mu[i];
}
}
for (int i=2;i<=maxn;++i) mu[i]+=mu[i-1];
}
LL solve(int l,int r){
int last; LL ret=0;
l/=k; r/=k;
for (int i=1;i<=l&&i<=r;i=last+1){
last=min(l/(l/i),r/(r/i));
ret+=1LL*(l/i)*(r/i)*(mu[last]-mu[i-1]);
}return ret;
}
int main(){
get_pre(); int T; scanf("%d",&T);
while (T--){
scanf("%d%d%d%d%d",&a,&b,&c,&d,&k);
printf("%lld\n",solve(b,d)-solve(a-1,d)-solve(b,c-1)+solve(a-1,c-1));
}
}
妈妈我终于会反演啦【雾
其实我就是来码个公式复习一下QwQ
前置技能:
那么按照题意,我们令
首先明确我们为什么要反演,因为我们要求一个不好直接算的东西,而他和一个容易计算的东西有关系,所以我们要同过那个容易计算的东西来反向计算出我们真正要求的东西
于是我们令
,其中?是整数,整个式子表示k的倍数
根据前置技能得
这时候发现d/k求起来很麻烦,换个元,令d'=d/k
于是
再来化简f函数得到
所以
所以
然后分块求就好了,线性筛出莫比乌斯函数并求出前缀和,注意容斥一下
#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define LL long long
using namespace std;
const int maxn=100010;
int a,b,c,d,k;
int p[maxn],N=0,mu[maxn];
bool vis[maxn];
void get_pre(){
mu[1]=1;
for (int i=2;i<=maxn-5;++i){
if (!vis[i]){
p[++N]=i;
mu[i]=-1;
}
for (int j=1;j<=N&&p[j]*i<=maxn-5;++j){
vis[p[j]*i]=1;
if (i%p[j]==0){mu[p[j]*i]=0;break;}
mu[p[j]*i]=-mu[i];
}
}
for (int i=2;i<=maxn;++i) mu[i]+=mu[i-1];
}
LL solve(int l,int r){
int last; LL ret=0;
l/=k; r/=k;
for (int i=1;i<=l&&i<=r;i=last+1){
last=min(l/(l/i),r/(r/i));
ret+=1LL*(l/i)*(r/i)*(mu[last]-mu[i-1]);
}return ret;
}
int main(){
get_pre(); int T; scanf("%d",&T);
while (T--){
scanf("%d%d%d%d%d",&a,&b,&c,&d,&k);
printf("%lld\n",solve(b,d)-solve(a-1,d)-solve(b,c-1)+solve(a-1,c-1));
}
}
相关文章推荐
- 1.m分解阶乘之和
- 2.几种递推数
- 3.欧拉函数
- 4.快速幂模m算法
- 5.扩展欧几里得&&中国剩余定理
- 6.数论_web
- Project Euler Problem 387 - Harshad Numbers - 深度优先
- 编程之美2015初赛A
- 数论题集
- BZOJ3275 Number (最小割)
- 原根
- 阶与原根学习笔记
- HDU 1299 Diophantus of Alexandria
- Leftmost Digit(HDU 1060)
- Rightmost Digit(HDU 1061)
- Python-在奇数中寻找素数
- ZOJ 2674 Strange Limit 欧拉定理
- LeetCode-Palindrome Number
- 组合数求模总结
- Sicily 1047 Super Snooker