BZOJ3994: [SDOI2015]约数个数和
2016-03-04 16:09
1441 查看
Description
设d(x)为x的约数个数,给定N、M,求#include<cstdio> #include<cctype> #include<queue> #include<cmath> #include<cstring> #include<algorithm> #define rep(i,s,t) for(int i=s;i<=t;i++) #define dwn(i,s,t) for(int i=s;i>=t;i--) #define ren for(int i=first[x];i;i=next[i]) using namespace std; const int BufferSize=1<<16; char buffer[BufferSize],*head,*tail; inline char Getchar() { if(head==tail) { int l=fread(buffer,1,BufferSize,stdin); tail=(head=buffer)+l; } return *head++; } inline int read() { int x=0,f=1;char c=getchar(); for(;!isdigit(c);c=getchar()) if(c=='-') f=-1; for(;isdigit(c);c=getchar()) x=x*10+c-'0'; return x*f; } typedef long long ll; const int maxn=50010; int mu[maxn],f[maxn],vis[maxn],pri[maxn],cnt; ll d[maxn]; void gen(int n) { mu[1]=1;d[1]=1; rep(i,2,n) { if(!vis[i]) pri[++cnt]=i,mu[i]=-1,d[i]=2,f[i]=1; rep(j,1,cnt) { if(i*pri[j]>n) break; vis[i*pri[j]]=1; if(i%pri[j]==0) { mu[i*pri[j]]=0; f[i*pri[j]]=f[i]+1; d[i*pri[j]]=d[i]/(f[i]+1)*(f[i]+2); break; } mu[i*pri[j]]=-mu[i]; d[i*pri[j]]=d[i]*d[pri[j]]; f[i*pri[j]]=1; } } rep(i,2,n) d[i]+=d[i-1],mu[i]+=mu[i-1]; } ll solve(int n,int m) { if(n>m) swap(n,m);ll ans=0; rep(i,1,n) { int last=min(n/(n/i),m/(m/i)); ans+=(mu[last]-mu[i-1])*d[n/i]*d[m/i]; i=last; } return ans; } int main() { gen(50000); dwn(T,read(),1) printf("%lld\n",solve(read(),read())); return 0; }
View Code
相关文章推荐
- STL---容器
- 比较难的日语单词
- NSMutableAttributedString字体颜色等属性
- 在工作表多行多列的区域中查找目标值的位置
- web项目引用WEB-INF下面的图片例子
- 【反编译】的下载与使用方法 。
- 使用CXF 做webservice简单例子
- matlab中如何获取当前日期时间函数
- Tomcat在eclipse里面能正常启动,而在浏览器中访问http://localhost:8080/不能访问,且报404错
- VMware三种联网方式
- iOS之轻松上手block
- 练习1-6
- oracle中导出表的结构和数据
- 修改input-placeholder 默认色
- LINK:fatal error LNK1123: 转换到 COFF 期间失败: 文件无效或损坏
- c# uploadimage - httphelper
- OC_协议实现代理委托模式
- one-hot enconding,类别型特征,非数值型特征。
- LoRaWAN 规范 1.0(2~4章)
- IntelliJ IDEA使用教程