多校赛上最近有道题目是gugufishtion

Problem Description

Today XianYu is too busy with his homework, but the boring GuGu is still disturbing him!!!!!! At the break time, an evil idea arises in XianYu’s mind. ‘Come on, you xxxxxxx little guy.’ ‘I will give you a function ϕ(x) which counts the positive integers up to x that are relatively prime to x.’ ‘And now I give you a fishtion, which named GuGu Fishtion, in memory of a great guy named XianYu and a disturbing and pitiful guy GuGu who will be cooked without solving my problem in 5 hours.’ ‘The given fishtion is defined as follow:

Gu(a,b)=ϕ(ab)ϕ(a)ϕ(b)Gu(a,b)=ϕ(ab)ϕ(a)ϕ(b)

And now you, the xxxxxxx little guy, have to solve the problem below given m,n,p.

(∑a=1m∑b=1nGu(a,b)(modp)(∑a=1m∑b=1nGu(a,b)(modp)

So SMART and KINDHEARTED you are, so could you please help GuGu to solve this problem? ‘GU GU!’ GuGu thanks.

看题解的时候看到这个通过观察,感到十分崩溃

后来发现原来这个φ 是欧拉函数

欧拉函数,就是求小于n的与n互质的数的数目

通式是:

其中p1, p2……pn为x的所有质因数,x是不为0的整数。

φ(1)=1(和1互质的数(小于等于1)就是1本身)。

注意:每种质因数只一个。 比如12=2*2*3那么φ(12)=12*(1-1/2)*(1-1/3)=4

若n是质数p的k次幂,

,因为除了p的倍数外,其他数都跟n互质。

设n为正整数,以 φ(n)表示不超过n且与n互素的正整数的个数,称为n的欧拉函数值

φ:N→N,n→φ(n)称为欧拉函数。

欧拉函数是积性函数——若m,n互质,

特殊性质:当n为奇数时,

, 证明与上述类似。

若n为质数则

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

28

29

30

/*

特性 :

1.若a为质数,phi[a]=a-1;

2.若a为质数,b mod a=0,phi[a*b]=phi[b]*a

3.若a,b互质,phi[a*b]=phi[a]*phi[b](当a为质数时,if b mod a!=0 ,phi[a*b]=phi[a]*phi[b])

*/

int m[n],phi[n],p[n],nump;

//m[i]标记i是否为素数,0为素数,1不为素数;p是存放素数的数组;nump是当前素数个数;phi[i]为欧拉函数

int make()

{

        phi[1]=1;

    for (int i=2;i<=n;i++)

    {

        if (!m[i])//i为素数

        {

            p[++nump]=i;//将i加入素数数组p中

            phi[i]=i-1;//因为i是素数,由特性得知    

        }    

        for (int j=1;j<=nump&&p[j]*i<n;j++)  //用当前已的到的素数数组p筛,筛去p[j]*i

        {

            m[p[j]*i]=1;//可以确定i*p[j]不是素数 

            if (i%p[j]==0) //看p[j]是否是i的约数,因为素数p[j],等于判断i和p[j]是否互质 

            {

                phi[p[j]*i]=phi[i]*p[j]; //特性2

                break;

            }

            else phi[p[j]*i]=phi[i]*(p[j]-1); //互质,特性3其,p[j]-1就是phi[p[j]]   

        }

    }

}

莫比乌斯(Möbius)函数

  对于每个正整数n(n ≥ 2),设它的质因数分解式为:

  

  根据这个式子定义n的莫比乌斯函数为:

  

  也就是如果n有平方因子,则为0. 否则是-1的质因数个数次方。

  举个简单的例子:6 = 2 × 3,所以;  9 = 3×3, 所以 

【命题一】

  对于正整数n有:

  

  也就是n>2时,所有n的约数对应函数值之和为0.

莫比乌斯反演

  若定义在正整数集上的两个函数,f(n)和g(n)满足对任意n有:

        (1)

  

  则可以通过f来表示g:

      (2)

  反之,亦可以由关系(2)得到(1)

莫比乌斯函数和欧拉函数的关系:

  

转载于:https://www.cnblogs.com/wyboooo/p/9643390.html

欧拉函数和莫比乌斯函数相关推荐

  1. Comet OJ - Contest #8 E神奇函数(莫比乌斯函数容斥)

    Comet OJ - Contest #8 E神奇函数(莫比乌斯函数容斥) 题目大意 定义d(x)d(x)d(x)的值为x的最小素因子,定义 f(x)={1x=1d(x)f(xd2(x))x>1 ...

  2. 筛表合集(素数筛 欧拉函数筛 莫比乌斯函数筛)

    [目录] 一.素数筛 1.素数判断 2.素数普通筛 3.素数线性筛 4.素数区间筛 二.欧拉函数筛 三.莫比乌斯函数筛 [素数筛] 1.直接判定质数 bool judgePrime( int num ...

  3. 容斥原理与欧拉函数与莫比乌斯函数,狄利克雷卷积与莫比乌斯变换,反演

    莫比乌斯函数可以看成是一种被内化了的容斥原理,许多数论上的结论定理根据容斥原理和数学归纳法可以推导出来,但是有关容斥原理的表达式的构造往往并不容易,运气不好很难找到,而莫比乌斯函数则是巧妙的把容斥原理 ...

  4. 专题·莫比乌斯函数与欧拉函数【including 整除分块,积性函数,狄利克雷卷积,欧拉函数,莫比乌斯函数,莫比乌斯反演

    初见安~又是好久没写博客了--加上CSP才炸了一波. 目录 一.整除分块 题解 二.积性函数 三.狄利克雷卷积 四.欧拉函数 五.莫比乌斯函数(mu) 六.莫比乌斯反演 一.整除分块 看个例题:洛谷P ...

  5. 线性筛与欧拉函数、莫比乌斯函数

    网上关于素数筛的资料很多,这里只是给出弱鸟整理的几个线性筛和应用. 最朴素的素数筛--埃拉托斯特尼筛法(Sieve of Eratosthenes) 复杂度 Olognlogn int primes[ ...

  6. c莫比乌斯函数_莫比乌斯函数总结

    莫比乌斯函数总结 性质:\(\sum_{d|n}\mu(d)=[n==1]\) 这个可以用组合数的性质来证,形象点的话就是杨辉三角. 因为恒等式:\(\sum_{i=0}^{n}(-1)^nC_{n} ...

  7. 25 欧拉积分: (伽马)函数、(贝塔)函数

  8. [莫比乌斯反演]莫比乌斯函数

    莫比乌斯函数定义 μ ( n ) = { 1 n = 1 ( − 1 ) k n = p 1 p 2 p 3 - p k 0 p 2 ∣ n \mu(n)=\begin{cases}1&n=1 ...

  9. python莫比乌斯_莫比乌斯函数 - osc_7eqzxl4g的个人空间 - OSCHINA - 中文开源技术交流社区...

    前导 要学习莫比乌斯函数 需要学习 到 积性函数,深度理解欧拉筛. 先说说什么是积性函数吧. 积性函数 其实积性函数非常好理解, 定义 积性函数:若gcd(a,b)=1,且满足f(ab)=f(a)f( ...

最新文章

  1. 金九银十,不要跳槽!
  2. python将二维列表内容写入和读取.csv文件
  3. 线上会议丨中国中文信息学会2020学术年会将于12月27日举行
  4. linux scp远程复制命令使用
  5. 开源杀毒引擎 ClamAV
  6. python爬虫小说代码示例-中文编程,用python编写小说网站爬虫
  7. c++求矩阵的秩_一篇文章搞定矩阵相关概念及意义通俗解释汇总
  8. OSPF邻接形成过程
  9. “敏捷版”全链路压测
  10. emc整改措施及案例_12.5EMC一点通||EMC整改措施
  11. linux目标机无法连接到网络,linux 安装了xunsearch服务,但是PHP连接显示:目标计算机积极拒绝...
  12. TCP/IP协议详解:IP、ARP、RARP、ICMP、IGMP
  13. wpsa3版面怎么变成a4、排版不变_wpsa4变a3排版打印
  14. 俄亥俄州立大学哥伦布分校计算机科学,【精选】俄亥俄州立大学哥伦布分校计算机科学与工程本科专业申请.pdf...
  15. 百度地图点击城市的出现下钻
  16. 使用Python自动化管理Exchange邮箱
  17. autoit-v3 编辑工具SciTE Script Editor 解决中文乱码+中文换行崩溃问题
  18. 视频流媒体播放器EasyPlayer-RTSP-Android 如何随意切换播放视频流?
  19. vue父组件中调用子组件中的方法 及vue父组件调用孙子组件的方法
  20. 3GPP TS 23501-g51 中英文对照 | 4.2.5a Radio Capabilities Signalling optimisation

热门文章

  1. MDCC创新应用展 应用秀场之《程序员》杂志封面人物秀(组图)
  2. Bug改不完,迭代总延期,咋办?
  3. 终止一个TCP连接要经过四次挥手内容
  4. 计算机专业英语app推荐信,计算机系教授英文推荐信_自荐信.docx
  5. Linux c 八种排序
  6. pahoMQTT移植
  7. IEEEtrans论文LaTeX排版
  8. html doctype 作用介绍
  9. Android事件传递可以这样理解
  10. javascript延迟加载方式