bzoj 4173: 数学 (数论+sqrt求phi)
2017-02-16 21:59
344 查看
4173: 数学
Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 456 Solved: 224
[Submit][Status][Discuss]
Description
Input
输入文件的第一行输入两个正整数 。Output
如题Sample Input
5 6Sample Output
240HINT
N,M<=10^15Source
[Submit][Status][Discuss]
题解:化式子(其实是可以找规律的嘛)+sqrt求phi
这道题找规律是个好方法,最后的答案就是phi(n)*phi(m)*m*n%p
但是还是需要知道应该知道推导。。。
n=x*k+r1,m=y*k+r2
n%k=r1 n/k=x
m%k=r2 m/k=y
m+n=(x+y)*k+(r1+r2)
由r1+r2>=k推得(n+m)%k=r1+r2-k,(n+m)/k=q1+q2+1
所以floor((n+m)/k)-floor(n/k)-floor(m/k)=1,不满足条件的k对应的floor((n+m)/k)-floor(n/k)-floor(m/k)=0
然后我们在看一下式子
sigma(k∈S)phi(k)=sigma(k=1..n+m) phi(k)*(floor((n+m)/k)-floor(n/k)-floor(m/k))
=simga(k=1..n+m)phi(k)*floor((n+m)/k)-sigma(k=1..n+m)phi(k)*floor(n/k)-sigma(k=1..n+m)phi(k)*floor(m/k)
如果k>m,floor(m)=0,n也是同理
进一步化简:simga(k=1..n+m)phi(k)*floor((n+m)/k)-sigma(k=1..n)phi(k)*floor(n/k)-sigma(k=1..m)phi(k)*floor(m/k)
接下来怎么搞呢?
sigma(k=1..n)i =sigma(k=1..n) sigma(i|k)phi(i)=sigma(i=1..n) phi(i)*floor(n/i)
然后发现,这个式子就是上面相似的三部分耶
所以式子最终化成:
sigma(k∈S)phi(k)=sigma(k=1..n+m)k-sigma(k=1..n)k-sigma(k=1..m)k
=n*m
phi(n)*phi(m)*sigma(k∈S)phi(k)=phi(n)*phi(m)*n*m%p 得证。
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define LL long long
#define p 998244353
using namespace std;
LL n,m;
LL phi(LL x)
{
LL k=x; LL ans=x;
for (LL i=2;i*i<=k;i++)
if (k%i==0) {
ans=ans/i*(i-1);
while (k%i==0) k/=i;
}
if (k>1) ans=ans/k*(k-1);
return ans%p;
}
int main()
{
freopen("a.in","r",stdin);
freopen("my.out","w",stdout);
scanf("%I64d%I64d",&n,&m);
LL ans=phi(n)*phi(m)%p; n%=p; m%=p;
printf("%I64d\n",ans%p*n%p*m%p);
}
相关文章推荐
- 【bzoj4173】【数学】【数论】
- [BZOJ4173]数学(数论)
- 【BZOJ 4173】 数学 数论
- BZOJ 4173 数学 数论
- BZOJ4173 数学
- 【bzoj 4173】数学
- BZOJ3643 Phi的反函数(数论+搜索)
- BZOJ 4173: 数学
- [BZOJ3643]Phi的反函数(数论+dfs)
- bzoj4173 数学
- 【BZOJ2742】【HEOI2012】Akai的数学作业 [数论]
- 【bzoj4173】数学
- BZOJ 2142 礼物 组合数学+数论
- [莫比乌斯反演 高斯消元 数学技巧] BZOJ 3601 一个人的数论
- BZOJ-4173 数学 欧拉函数+关于余数的变换
- 【BZOJ4173】数学 欧拉函数神题
- BZOJ 4173 数论
- BZOJ3505 & 洛谷P3166 [Cqoi2014]数三角形 【数学、数论】
- BZOJ 2956 模积和 (数学推导+数论分块)
- BZOJ5323 [Jxoi2018]游戏 【数论/数学】