[BZOJ 3994][SDOI 2015]约数个数 数学+反演
来源:互联网 发布:淘宝返还网 编辑:程序博客网 时间:2024/04/29 22:12
主要就是这个公式吧,其他的都是很好写的莫比乌斯反演,一股脑的推就好了.对于这个式子如果想要看严格的证明的话就去看po姐的博客吧,很好懂的,但是后来我想出来一个感性认识的,不过我觉得肯定会被打脸,没事,我已经做好准备了
感性证明:
计算一个整数x的约数有多少个,我们可以这样来算,把x因数分解:
x=
这样x的约数个数就是(1+a1)(1+a2)(1+a3)……因为对于x的每一个约数都可以用分解出来的素数和不同次方相乘得到,而每一个数都可以有ai+1种选择.例如12=2^2 * 3,那么12的约数就有
1*1=1
1*2=2
1*4=4
1*3=3
2*3=6
4*3=12
联系刚才的性质可以很显然的发现12的每一个约数都可以由两个互质的数相成得到(假设1与所有数互质),因为对于x因数分解出来的每一个因数都是两两互质的,所以我们可以枚举互质的数,计算这个数可以在哪些数中间造成贡献,即是上面的公式.
严格的证明的话还是%%po姐吧,我就YY一下
#include<cstdio>#include<cstring>#include<iostream>#define maxn 50021#define LL long longusing namespace std;int n,m,mu[maxn],vis[maxn],pri[maxn],cnt;int f[maxn],g[maxn],T;void pre(){ mu[1]=1,f[1]=1; for(int i=2;i<maxn;i++){ if(!vis[i]){pri[++cnt]=i;mu[i]=-1,f[i]=2,g[i]=1;} for(int j=1;j<=cnt&&pri[j]*i<maxn;j++){ vis[pri[j]*i]=1; mu[pri[j]*i]=-mu[i]; g[pri[j]*i]=1; f[pri[j]*i]=f[i]*2; if(i%pri[j]==0){ mu[i*pri[j]]=0; g[pri[j]*i]=g[i]+1; f[pri[j]*i]=f[i]/(g[i]+1)*(g[i]+2); break; } } } for(int i=2;i<maxn;i++)f[i]+=f[i-1],mu[i]+=mu[i-1];}void solve(){ LL ans=0; for(int j,i=1;i<=n;i=j+1){ j=min(n/(n/i),m/(m/i)); ans+=(LL)(mu[j]-mu[i-1])*f[n/i]*f[m/i]; } printf("%lld\n",ans);}int main(){ pre(); scanf("%d",&T); while(T--){ scanf("%d%d",&n,&m); if(n>m)swap(n,m); solve(); } return 0;}
0 0
- [BZOJ 3994][SDOI 2015]约数个数 数学+反演
- 【SDOI 2015】【BZOJ 3994】约数个数和
- bzoj 3994: [SDOI2015]约数个数和(反演)
- BZOJ 3994 约数个数和(莫比乌斯反演)
- [莫比乌斯反演 约数个数] BZOJ 3994 [SDOI2015]约数个数和
- BZOJ 3994 Sdoi2015 约数个数和 莫比乌斯反演
- BZOJ 3994 [SDOI2015]约数个数和 (神定理+莫比乌斯反演)
- BZOJ 3994: [SDOI2015]约数个数和|莫比乌斯反演
- BZOJ 3994 [SDOI2015]约数个数和 (莫比乌斯反演)
- BZOJ 3994: [SDOI2015]约数个数和 莫比乌斯反演
- BZOJ 3994 [SDOI2015]约数个数和 莫比乌斯反演
- 3994: [SDOI2015]约数个数和 莫比乌斯反演
- [bzoj 3994] sdoi2015 约数个数和
- BZOJ 3994 [SDOI2015]约数个数和
- bzoj 3994 [SDOI2015]约数个数和
- BZOJ 3994 [SDOI2015]约数个数和
- bzoj 3994: [SDOI2015]约数个数和
- 【bzoj3994】【SDOI2015】约数个数和【数论】【反演】
- Fragment中设置滚动条随页面滑动的方法
- Jquery选择以什么开头、结尾,包含什么的对象
- jQuery选择器之基本筛选选择器
- Windows下计算文件摘要信息MD5、SHA-1、SHA-256
- 骑马修栅栏
- [BZOJ 3994][SDOI 2015]约数个数 数学+反演
- js keyup事件 在安卓和苹果手机不同使用 input propertychange
- Install MDX on Mac
- 字符编码对程序的影响及分析
- 一个简陋的自定义用线程Thread实现的计时器Timer,可以增加定时时间
- php生成pdf文件,html转pdf文件的最佳方法(转载)
- 设计模式 命令模式 之 管理智能家电
- Python发送邮件(常见四种邮件内容)
- UE4中切换关卡