BZOJ 4173 数论
2024-10-01 00:55:12
思路:
$(m%k+n%k>=k) *phi(k)$
$我们不妨设n=q_1k+r_1 m=q_2k+r$2
$n+m=(q_1+q_2)k+r1+r2$
${\lfloor}\frac{n+m}{k}{\rfloor}-{\lfloor}\frac{m}{k}{\rfloor}-{\lfloor}\frac{n}{k}{\rfloor}=(m%k+n%k>=k)$
$原式=phi(k)*({\lfloor}\frac{n+m}{k}{\rfloor}-{\lfloor}\frac{m}{k}{\rfloor}-{\lfloor}\frac{n}{k}{\rfloor})$
$id=phi|1$
$n=\Sigma_{d|n}phi(d)$
$原式=\Sigma_{i=1}^{n+m}i-\Sigma_{i=1}^mi-\Sigma_{i=1}^ni$
$ =(n+m)*(n+m-1)/2+m*(m-1)/2+n*(n-2)/2$
$ =n*m$
//By SiriusRen
#include <cstdio>
using namespace std;
typedef long long ll;
ll n,m,mod=;
ll phi(ll x){
ll res=;
for(int i=;1LL*i*i<=x;i++){
if(x%i==){
while(x%i==)x/=i,res=res*i;
res=res/i*(i-);
}
}if(x!=)res=res*(x-);
return res;
}
int main(){
scanf("%lld%lld",&n,&m);
printf("%lld\n",((((phi(n)%mod)*(phi(m)%mod))%mod*(n%mod))%mod*(m%mod))%mod);
}
最新文章
- MQ通道配置
- python错误类型
- Hibernate - list()和iterate()的区别
- [转]怎样在cmd(命令提示符)下进行复制粘贴操作
- C# unsafe code
- MUI 列表页面绑定接口数据
- WebKit介绍及总结(一)
- nohup及/dev/null使用
- 浅谈MVC MVP MVVM
- Spring学习之二
- PO订单审批通过API
- 隐马尔可夫模型(HMM)总结
- data.table包使用应该注意的一些细节
- docker基本概念
- Python自学:第三章 使用函数sort( )对列表进行临时排序
- WPF多屏最大化
- shopnc 手机网站配置
- ASP.NET MVC中的Session设置
- linux远程方式,以及基础命令
- 聊聊javascript的null和undefined
热门文章
- Python单例模式的实现方式
- Spring AOP 学习(五)
- Apache 流框架 Flink,Spark Streaming,Storm对比分析(2)
- https://blog.csdn.net/u011495642/article/details/79958444
- HDU 1234 简单模拟题
- noip模拟赛 天天和不可描述
- nyoj_528_找球号(三)_201404152050
- 线程间的通信----wait/notify机制
- [bzoj1563][NOI2009]诗人小G(决策单调性优化)
- Spring MVC中@RequestMapping注解使用技巧(转)