BZOJ_P3561 DZY Loves Math VI(数论+莫比乌斯反演)
来源:互联网 发布:java图形界面时间控件 编辑:程序博客网 时间:2024/04/30 14:42
BZOJ传送门
Time Limit: 10 Sec Memory Limit: 256 MB
Submit: 90 Solved: 66
[Submit][Status][Discuss]
Description
给定正整数n,m。求
Input
一行两个整数n,m。
Output
一个整数,为答案模1000000007后的值。
Sample Input
5 4
Sample Output
424
HINT
数据规模:
1<=n,m<=500000,共有3组数据。
Source
By Jcvb
#include<cstdio>#include<iostream>using namespace std;#define N 500005#define P 1000000007typedef long long ll;bool b[N];int mu[N],p[N];int n,m,cnt,ans;int sum[N],a[N];inline int Pow(int a,int b,int res=1){while(b){if(b&1) res=(ll)res*a%P;a=(ll)a*a%P,b>>=1;}return res;}void GetMu(){ mu[1]=1; for(int i=2;i<=n;i++){ if(!b[i]) mu[i]=-1,p[cnt++]=i; for(int j=0;j<cnt&&i*p[j]<=n;j++){ b[i*p[j]]=1; if(i%p[j]) mu[i*p[j]]=-mu[i]; else break; } }}int main(){ scanf("%d%d",&n,&m);if(n<m) swap(n,m);GetMu(); for(int i=1;i<=n;i++) a[i]=1; for(int d=1;d<=m;d++){ int x=Pow(d,d),y=0,limm=m/d,limn=n/d; for(int j=1;j<=limn;j++) a[j]=(ll)a[j]*j%P,sum[j]=(sum[j-1]+a[j])%P; for(int k=1;k<=limm;k++)if(mu[k]) y=((ll)a[k]*a[k]%P*sum[limm/k]%P*sum[limn/k]%P*mu[k]+y+P)%P; ans=(ans+(ll)x*y%P)%P; } printf("%d",ans); return 0;}
1 0
- BZOJ_P3561 DZY Loves Math VI(数论+莫比乌斯反演)
- [数论][莫比乌斯反演][杜教筛] BZOJ 3512: DZY Loves Math IV
- BZOJ 3561 DZY Loves Math VI(莫比乌斯反演)
- BZOJ 3309 DZY Loves Math 莫比乌斯反演
- 3309: DZY Loves Math 莫比乌斯反演
- [BZOJ3309] DZY Loves Math - 莫比乌斯反演
- bzoj 3309: DZY Loves Math 莫比乌斯反演
- bzoj3561 DZY Loves Math VI 莫比乌斯函数
- BZOJ 3309 DZY Loves Math (莫比乌斯反演的应用 好题)
- bzoj3739 DZY loves math VIII 莫比乌斯函数
- bzoj3309 DZY Loves Math 数论
- BZOJ3561: DZY Loves Math VI
- 【bzoj3561】DZY Loves Math VI
- [bzoj3561] DZY Loves Math VI
- BZOJ3561: DZY Loves Math VI
- bzoj3561 DZY Loves Math VI
- bzoj 3309: DZY Loves Math (反演)
- BZOJ 3560 DZY Loves Math V 数论
- py解析xml
- Java基础面试点(1)
- 4198: [Noi2015]荷马史诗 贪心+优先队列
- zhuan:串口成帧协议
- iOS 代码修改xib中的约束
- BZOJ_P3561 DZY Loves Math VI(数论+莫比乌斯反演)
- SlidingMenu的应用实例
- 文件的读写
- 2016年司考新手复习指南
- Adjacent Node Sum (图的邻接表表示 C语言)
- LeetCode : Merge Two Sorted Lists [java]
- 网络编程遇到的一些问题
- 分布式 session
- js程度