码迷,mamicode.com
首页 > 其他好文 > 详细

bzoj4173数学(欧拉函数)

时间:2018-08-25 00:40:09      阅读:193      评论:0      收藏:0      [点我收藏+]

标签:ssi   次数   src   floor   define   details   scan   targe   const   

技术分享图片

解答:https://blog.csdn.net/popoqqq/article/details/46820313

   https://blog.csdn.net/zhhx2001/article/details/52300924

技术分享图片

倒数第二行的拓展(借用某群提问的一个问题):

技术分享图片

对于d*phi(d),当d为k时,在左式,他出现在 i 为k,2k,3k...(n/k)*k  共出现n/k次,在右式,他出现在 n/i 为n/1,n/2...n/(n/k) 共出现n/k次,画个表格,每行就是每个 i 对应1,2,3...floor(n/i),i 从1~n画出后就是一个逐渐缩减的三角形,每列对应一个d的值,每个d出现在 i 为floor(5/i)>=3的地方,终止在floor(5/i)=3的地方,所以每个d出现次数是 floor(n/i)

技术分享图片左式在图中的体现是把每个纵列d的每个安放位置调整一下,例如②,第一个2在i==1,第二个2调整到i==4,总出现次数不变(即看n中有几个大小为2的块

 

这里d用两种不同的方式出现在这些个 i 中,分析第一种有n/k个值能d|i 时,第二种就是把这个方式解释一下 ~~~

而等号右边换成技术分享图片同样成立,证明更简单了,就是参数为k时出现n/k次!!!

 

第一个等式用莫比乌斯反演列出f(1),f(2)...f(n)那些个等式,就很容易发现F(n)=n

 

小总结:模‘%’  整除‘|’ 这些抽象符号尽量化简成+-*/ 的形式!!!

ac代码:

#include<bits/stdc++.h>
#define per(i,a,b) for(int i=a;i<=b;i++)
#define mod 998244353//注意mod最好define,不要定义类型,否则需要ll的时候会错,不过呢,ll的时候int也会变成ll把
using namespace std;
typedef long long ll;
const int inf =0x3f3f3f3f;
const double eps=1e-8;
//#define siz 1000005
ll read()
{
    ll res=0,f=1;
    char ch=getchar();
    while(ch<1||ch>9){f=ch==-?-1:1;ch=getchar();}
    while(ch>1&&ch<9){res=res*10+ch-0;ch=getchar();}
    return res*f;
}
ll n,m;
ll phi(ll n)
{
    ll i,re=n;
    for(i=2;i*i<=n;i++){
        if(n%i==0){
            re=re/i*(i-1);//注意先出后乘,如果先乘可能会超出数据范围
            while(n%i==0){n/=i;}//n除完素因子后为1
        }
    }
    if(n!=1){re=re/n*(n-1);}
    return re;
}
int main()
{
    scanf("%lld%lld",&n,&m);
    printf("%lld\n",phi(n)%mod*(phi(m)%mod)%mod*(n%mod)%mod*(m%mod)%mod);
    return 0;
}

 

bzoj4173数学(欧拉函数)

标签:ssi   次数   src   floor   define   details   scan   targe   const   

原文地址:https://www.cnblogs.com/WindFreedom/p/9532418.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!