摘要：
const int MAXN = 1000000;bool check[MAXN+10];int prime[MAXN+10];int mu[MAXN+10];void Moblus(){ memset(check,false,sizeof(check)); ... 阅读全文
摘要：
题意：求满足gcd(x,n)>=m的x的个数，其中，1=m则ans+=phi(n/i)，判断s*s是否等于n即如果是完全平方数的时候，只算一遍；如果不是完全平方数，为了计算根号n后面的s，比如6%2=0，则6%3也等于0，所以如果n/s>=m且不是完全平方数则有ans... 阅读全文
摘要：
参考大神https://blog.csdn.net/a27038/article/details/77203892 阅读全文
摘要：
直接套欧拉函数模板，链接：欧拉函数#include#include#include#include#include#include#include#include#include#include#include#include#include#include#incl... 阅读全文
摘要：
欧拉函数有直接求法和打欧拉函数表法。欧拉函数的定义：对正整数n，欧拉函数是少于或等于n的数中与n互质的数的数目。例如euler(8)=4，因为1,3,5,7均和8互质。Euler函数表达通式：euler(x)=x(1-1/p1)(1-1/p2)(1-1/p3)(1-1... 阅读全文
摘要：
思路：设第x天高峰同时出现，则x≡p（mod23）,x≡e(mod28),x≡i(mod33)，根据中国剩余定理，求出乘法逆元，直接套公式求和就行了。注意：求和后要减去初始天数d，如果是小于等于0，则要加上m=23*28*33#include#include#incl... 阅读全文
摘要：
ll kgcd(ll a,ll b,ll &x,ll &y)//扩欧{ if(!b) { x=1;y=0;return a; } ll t=kgcd(g,a%b,y,x); y-=a/b*x; return t;}ll... 阅读全文
摘要：
题意：求有多少对(x,y)满足x在x1~x2,y在y1~y2里使方程ax+by+c=0方法：c取相反数，然后扩欧解得x,y，再x*=c/d,y*=c/d求得一组解，计算x、y的改变量，x=b/d,y=a/d*-1。根据 x1ly||rx>ry||lp>rp,则输出0#... 阅读全文
摘要：
ll qpow(ll a,ll b,ll M) //快速幂{ ll ans =1; while(b) { if(b&1) ans*=a;ans%=M; a*=a;a%=M;b>>=1; } ... 阅读全文
摘要：
题意：给出一个长度为n的正整数序列，一个素序列是一个长度至少为2的连续的子序列，总和是大于或等于2的一个素数。例如给出3 5 6 3 8，存在两个长度为2的素序列（5+6=11，3+8=11），输出最短的序列长度和元素。方法：米勒拉宾测试素数或者线性筛+试除法，然后枚... 阅读全文