bzoj 4173: 数学 (数论+sqrt求phi)

4173: 数学

Time Limit: 10 Sec   Memory Limit: 256 MB
Submit: 456   Solved: 224
[ Submit][ Status][ Discuss]

Description

 

Input

 输入文件的第一行输入两个正整数 。 

Output

 如题

Sample Input

5 6

Sample Output

240

HINT

 N,M<=10^15

Source

[ 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);
}




  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值