摘要:
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),输出最短的序列长度和元素。方法:米勒拉宾测试素数或者线性筛+试除法,然后枚... 阅读全文