http://www.lydsy.com/JudgeOnline/problem.php?id=4407

题意:

给下N,M,K.求
思路:
 来自:http://blog.csdn.net/ws_yzy/article/details/50670213

#include<cstdio>
#include<cmath>
#include<iostream>
#include<algorithm>
#include<cstring>
#define ll long long
const ll Mod=;
ll f[],sum[],p[],s[];
bool mark[];
ll K,n,m;
ll read(){
ll t=,f=;char ch=getchar();
while (ch<''||''<ch){if (ch=='-') f=-;ch=getchar();}
while (''<=ch&&ch<=''){t=t*+ch-'';ch=getchar();}
return t*f;
}
ll Pow(ll x,ll y){
ll res=;
while (y){
if (y%) res=(res*x)%Mod;
y/=;
x=(x*x)%Mod;
}
return res;
}
void init(){
f[]=;
for (int i=;i<=;i++){
if (!mark[i]){
p[++p[]]=i;
s[p[]]=Pow(i,K);
f[i]=s[p[]]-;
}
for (int j=;j<=p[]&&p[j]*i<=;j++){
mark[p[j]*i]=;
if (i%p[j]==){
f[i*p[j]]=f[i]*s[j]%Mod;
break;
}
f[i*p[j]]=f[i]*f[p[j]]%Mod;
}
}
for (int i=;i<=;i++)
sum[i]=sum[i-]+f[i]%Mod;
}
int main(){
int T=read();K=read();
init();
while (T--){
n=read();m=read();
if (n>m) std::swap(n,m);
int j=;
ll ans=;
for (int i=;i<=n;i=j+){
j=std::min(n/(n/i),m/(m/i));
ans+=(((n/i)*(m/i)%Mod)*(sum[j]-sum[i-]))%Mod;
ans%=Mod;
}
printf("%lld\n",ans);
}
return ;
}

最新文章

  1. ABP(现代ASP.NET样板开发框架)系列之21、ABP展现层——Javascript函数库
  2. 通过js实现回到顶部功能
  3. Apache Sqoop - Overview——Sqoop 概述
  4. s3c2440液晶屏驱动 (内核自带) linux-4.1.24
  5. .Net程序员飞扬有用的85个工具
  6. 创建javascript对象的几种方式
  7. 年轻的团队Mono玩转Dalvik
  8. spring校验和文件上传
  9. sticky bit
  10. 基于Flume的美团日志收集系统(二)改进和优化
  11. 个人实验记录之EIGRP基本配置
  12. DemoExample
  13. 安装win7 32位系统出现的问题解决办法
  14. leetcode面试准备:Multiply Strings
  15. 【Linux】 诊断工具-strace
  16. Android 子线程请求ASP.NET后台
  17. Cookie中文乱码问题
  18. DevOps之平台架构
  19. Java代码风格和在idea中的一些设置
  20. LeetCode(5):最长回文子串

热门文章

  1. 夏普比率(Sharpe Ratio)
  2. powershell利用winform批量执行tsql语句
  3. SpringMVC进行文件的上传以及多文件的上传(转)
  4. Hibernate 关于load和get方法区别
  5. 遍历INI文件和删除指定域内容
  6. JavaScript MVC 框架[转载]
  7. atitit.提升研发效率的利器---重型框架与类库的差别与设计原则
  8. POJ2230 Watchcow【欧拉回路】
  9. Android入门之ActionBar实现Tab导航
  10. MM32初识(兼容STM32)