题目地址:51Nod 1256
题意:给出2个数M和N(M < N),且M与N互质,找出一个数K满足0 < K < N且K * M % N = 1,如果有多个满足条件的,输出最小的。
思路:K*M%N=1可以写成K*M-Y*N=1,这样公式就变成了扩展欧几里德求K值。因为是要求最小的,所以求出特解K以后,要变成(K%N+N)%N。#include
#include <math...
分类:
其他好文 时间:
2015-08-30 19:39:25
阅读次数:
122
求 φ(m!) * n! / m! .拆开得 m! *∏ (p-1)/p * n! / m! =∏ (p-1)/p * n! (p|m!)我们可以预处理出1~m的∏ (p-1)/p 和 n!. 注意到取模的数为质数, 那么就可以用乘法逆元和线性筛预处理前者.时间复杂度O(MAXN)-O(1)#inc...
分类:
其他好文 时间:
2015-08-28 10:48:11
阅读次数:
232
题目链接:
http://acm.hdu.edu.cn/showproblem.php?pid=5407
解题思路:
官方题解:
The problem is just to calculate g(N) =\ LCM(C(N,0), C(N,1), ..., C(N, N))g(N) = LCM(C(N,0),C(N,1),...,C(N,N)).
Introdu...
分类:
其他好文 时间:
2015-08-25 23:48:14
阅读次数:
276
Problem DescriptionConsider a positive integer X,and let S be the sum of all positive integer divisors of 2004^X. Your job is to determine S modulo 29...
分类:
移动开发 时间:
2015-08-21 01:40:15
阅读次数:
198
给出a和b,如果一个数每一位都是a或b,那么我们称这个数为good,在good的基础上,如果这个数的每一位之和也是good,那么这个数是excellent。求长度为n的excellent数的个数mod(1e9+7)。ps:1e9+7是一个质数。...
分类:
其他好文 时间:
2015-08-18 01:19:02
阅读次数:
242
求C(n,m)%mod的方法总结1.当n,m都很小的时候可以利用杨辉三角直接求。
C(n,m)=C(n-1,m)+C(n-1,m-1);2.利用乘法逆元。
乘法逆元:(a/b)%mod=a*(b^(mod-2)) mod为素数。
逆元可以利用扩展欧几里德或欧拉函数求得: 1).扩展欧几里德:b*x+p*y=1 有解,x就是所求 2).费马小定理:b^(p-1)=1(mod p),故b*b^...
分类:
其他好文 时间:
2015-08-16 23:08:39
阅读次数:
403
题目大意:
给一个整数N,求1~N中与N互质的数的4次方的和。
解题思路:
题目简单,过程有点复杂。理清思路就简单了。
利用公式1^4 + 2^4 + … + n^4 = n*(n+1)*(2*n+1)*(3*n^2+3*n-1)/30,可以求出
1^4 + 2^4 + … + n^4,除以30可以先求出30模M的逆元,然后将上式中除以30改为
乘以30的逆元。
再来求与n不互质的数的4次方和。将n质因数分解为 n = p1^a1* p2^a2 * … * pn^an
(其中p1、p2为质数)。
这样不互...
分类:
其他好文 时间:
2015-08-12 13:22:47
阅读次数:
354
一道非常有价值的题。
【解析1】欧几里德算法求乘法逆元,前缀和
[Analysis]O(T n log n)。
[Sum]
①int运算,如果会超出界,第一个数前要加上(LL)即类型转换。
②gcd不变的欧几里德定理:可以是加,也可以是减。
[Code]/**********************************************************...
分类:
其他好文 时间:
2015-07-24 13:02:22
阅读次数:
106
【解析】Burnside引理+背包dp+乘法逆元
[Analysis]
这道题卡了好久,就是没想懂置换跟着色是不一样的。
根据burnside引理,在一个置换群作用下不等价类的个数为每个置换作用下不动点个数的平均数。
在这道题中:
置换的对象 ——
每个状态,标号为1—N(这里的N不是题目的N,而是状态的个数)。
不动点 ——
前后染色状态完全相同的状态的个数。
...
分类:
其他好文 时间:
2015-07-22 22:48:48
阅读次数:
199
inv[x] = ( Mod - Mod /x ) * inv[Mod%x] % Mod设Mod=px+q. inv[x]=rpx+q = 0 (mod Mod)pxr+qr = 0 (mod Mod)p + qr = 0 (mod Mod)r=-p/q= -p*inv[q]= -(Mod/x)*i...
分类:
其他好文 时间:
2015-07-20 18:19:52
阅读次数:
90