thy的专栏

传送门：http://www.lydsy.com:808/JudgeOnline/problem.php?id=1477 思路：扩展欧几里德求解同余方程。 设x步后碰面，初始位置为x0，y0。 那么就有x0+mx=y0+nx(mod L) (m-n)x=y0-x0(mod L) //程序里要先把m-n变成正数，L已经是正数，不用变，负数模会有问题，这两步的等号是同余。 于是就可以把这个同

传送门：http://www.lydsy.com:808/JudgeOnline/problem.php?id=1407 思路：因为M 山洞数量满足要求，就是对于任何两个野人组成的同余方程最小解大于寿命短的野人的寿命，或者没有解。 代码：#include #include #include #define abs(a) (a<0?-a:a) const int maxn=100000

题目大意：求长度为n且每项均在[1,n]的不上升数列与不下降数列的个数和。 思路：总数就是不下降数列的个数*2-n(常数列的个数) 然后考虑不下降数列的个数 为了方便，把第0项设为0，把第n+1项设为n。 差分，然后不下降数列就是差分数组a[i]每一项大于等于0，且Σa[i]=n。 每项+1，就相当于在2n-1个空位（本来是2n+1，首尾不能放）放n个板子。 于是答案就是C(2n-1,

传送门：http://www.lydsy.com:808/JudgeOnline/problem.php?id=2705 思路：首先要求的是Σgcd（i,n）(1 于是可以化成ΣΣd*[gcd(i,n/d)==1]    (d|n,1 把d提出来 Σd*Σ[gcd(i,n/d)==1]  (d|n,1 右边的不就是phi（欧拉函数）吗 于是有Σd*phi(n/d)  (d|n) 因为

传送门：http://www.spoj.com/problems/GCDEX/ 思路：令g(n)=ΣΣgcd(i,j)(i 很显然ans=g(n)-(n+1)*n/2 那么g(n)=Σf(i)(i f(n)在上一题已证为积性函数，线筛求出，就可以在O(n)时间打出表来。 为了方便先设mindiv[i]为i的最小质因子，T[i]为i的最小质因子次数，d[i]=mindiv[i]^t[i]

传送门：http://www.lydsy.com:808/JudgeOnline/problem.php?id=1101 思路：设a 求的就是f(a,b,d)=ΣΣ[gcd(i,j)==d](1 又有公式Σmiu(d) =[n==1] (d|n) //miu是莫比乌斯函数 把n用gcd(i,j)代入得 ΣΣΣmiu(d)(d|gcd(i,j),1 =Σmiu(d)*(a/d)*(b/d

传送门：http://www.lydsy.com:808/JudgeOnline/problem.php?id=2820 思路： #include #include #include const int N=10000010; using namespace std; typedef long long ll; int mu[N+10],g[N+10],pri[

传送门：http://www.lydsy.com/JudgeOnline/problem.php?id=3667 思路：首先我们说说Miller_Rabin算法 我们发现了费马小定理 那它倒过来对不对呢 如果a^(p-1)=1(mod p),那么p一定是素数吗？ 很不幸，是错的 虽然出错概率很低，但是可以被卡 于是我们就给它打补丁 我们又找到了一个二次探测的方法 如果p是质数，那