数论
文章平均质量分 56
SSL_MYD
这个作者很懒,什么都没留下…
展开
-
YBT金牌导航元素
YBT金牌导航元素#include<iostream> #include<cstdio>#include<algorithm>using namespace std;struct jgt{ long long id,magic;}a[1010];long long d[64];bool cmp(jgt t1,jgt t2){ return t1.magic>t2.magic;}bool add(long long x){ for(lo原创 2021-06-22 18:22:58 · 178 阅读 · 0 评论 -
YBT金牌导航彩灯
YBT金牌导航彩灯#include<iostream> #include<cstdio>using namespace std;long long d[110],mods=2008;void add(long long x){ long long i; for(i=63;i>=0;--i) if(x&(1ll<<i)) if(d[i])x^=d[i]; else{d[i]=x;return;} return;}int m原创 2021-06-21 18:32:19 · 77 阅读 · 0 评论 -
YBT高效进阶 6.2.5 余数之和
YBT余数求和思路(n,k)(n,k)(n,k)=∑i=1nk mod i=\sum_{i=1}^nk\bmod i=∑i=1nkmodi=∑i=1nk−⌊ki⌋=\sum_{i=1}^nk-\left\lfloor \dfrac ki \right\rfloor=∑i=1nk−⌊ik⌋=n∗k−∑i=1n⌊ki⌋=n*k-\sum_{i=1}^n\left\lfloor \dfrac ki \right\rfloor=n∗k−∑i=1n⌊ik⌋分块i∈[x,⌊k⌊kx⌋⌋]时⌊ki原创 2021-06-10 13:59:53 · 116 阅读 · 0 评论 -
洛谷P4777【模板】扩展中国剩余定理(EXCRT)
洛谷P4777【模板】扩展中国剩余定理(EXCRT)TITLE思路中国剩余定理CRTx≡a1(modm1)x\equiv a_1\pmod{m_1}x≡a1(modm1)………………x≡an(modmn)x\equiv a_n\pmod{m_n}x≡an(modmn)gcd(mi,mj)=1,(i≠j)\gcd(m_i,m_j)=1,(i\not=j)gcd(mi,mj)=1,(i=j)扩展中国剩余定理EXCRTx≡a1(modm1)x\equiv a_1\pmod{m原创 2021-05-22 08:21:06 · 122 阅读 · 0 评论 -
洛谷P3868 [TJOI2009]猜数字
洛谷P3868 [TJOI2009]猜数字TITLE思路中国剩余定理(CRT)用龟速乘,注意a[i]=(a[i]%m+m)%mCODE#include<iostream> #include<cstdio>using namespace std;void exgcd(long long a,long long b,long long &x,long long &y){ if(!b){x=1,y=0;return;} exgcd(b,a%b,x,y原创 2021-05-21 22:07:12 · 93 阅读 · 0 评论 -
洛谷P1495 【模板】中国剩余定理(CRT)/曹冲养猪+扩展欧几里得求逆元
P1495 【模板】中国剩余定理(CRT)/曹冲养猪TITLE思路fadjgoajdfiajbjaorwdCODE原创 2021-05-21 20:54:19 · 136 阅读 · 0 评论 -
洛谷P1447 [NOI2010] 能量采集
洛谷P1447 [NOI2010] 能量采集TITLE思路∑i=1n∑j=1m2gcd(i,j)−1=−n∗m+2∑i=1n∑j=1mgcd(i,j)\sum_{i=1}^n\sum_{j=1}^m2\gcd(i,j)-1=-n*m+2\sum_{i=1}^n\sum_{j=1}^m\gcd(i,j)∑i=1n∑j=1m2gcd(i,j)−1=−n∗m+2∑i=1n∑j=1mgcd(i,j)∵∑k∣nϕ(k)=n\because \sum_{k|n}\phi(k)=n∵∑k∣nϕ(k)原创 2021-05-20 17:20:13 · 93 阅读 · 0 评论 -
P5431 【模板】乘法逆元2
P5431 【模板】乘法逆元2TITLE思路线性求逆元预处理前缀积,后缀积扩展欧几里得算出所有数的积的逆元CODE#include<iostream> #include<cstdio>using namespace std;long long a[5000010],b[5000010],c[5000010];inline long long read()//快读 { long long ans=0; char u=getchar(); while (u&原创 2021-05-20 16:06:48 · 81 阅读 · 0 评论 -
洛谷P1516 青蛙的约会(扩展欧几里得EXGCD)
洛谷P1516 青蛙的约会TITLE思路ax+by=c存在整数解当且仅当gcd(a,b)|c扩展欧几里得算法实现求解二元一次不定方程的整数解若求出ax’+by’=gcd(a,b),那么将等式两边变为原来的c/gcd(a,b)倍long long exgcd(long long a,long long b,long long &x,long long &y){ if(!b){x=1,y=0;return a;} long long ans=exgcd(b,a%b,x,y),t原创 2021-05-19 14:05:09 · 105 阅读 · 0 评论 -
洛谷P4139 上帝与集合的正确用法
洛谷P4139 上帝与集合的正确用法TITLE思路欧拉定理ab≡abmod ϕ(m)+ϕ(m)a^b \equiv a^{b\mod \phi(m)+\phi(m)}ab≡abmodϕ(m)+ϕ(m)代入公式计算即可若ab,b<=2,return 0否则return 2ϕ(x)+函数自己(ϕ(x))mod b2^{\phi(x)+函数自己(\phi(x))}\mod b2ϕ(x)+函数自己(ϕ(x))modblong long ExEulerTheorem(long long x原创 2021-05-18 21:24:11 · 85 阅读 · 0 评论 -
洛谷P5253 [JSOI2013]丢番图
洛谷P5253 [JSOI2013]丢番图TITLE思路1/x+1/y=1/n1/x+1/y=1/n1/x+1/y=1/n化简得(x−n)(y−n)=n2化简得(x-n)(y-n)=n^2化简得(x−n)(y−n)=n2ans=(n2的因数个数+1)/2ans=(n^2的因数个数+1)/2ans=(n2的因数个数+1)/2ans=(∏(n的每个质因数的指数∗2+1)+1)/2ans=(\prod{(n的每个质因数的指数*2+1)}+1)/2ans=(∏(n的每个质因数的指数∗2+1)+1)/2原创 2021-05-18 21:11:05 · 84 阅读 · 0 评论 -
洛谷P5495 Dirichlet 前缀和
洛谷P5495 Dirichlet 前缀和TITLE思路∵质数唯一分解\because 质数唯一分解∵质数唯一分解∴每个数对倍数产生贡献\therefore 每个数对倍数产生贡献∴每个数对倍数产生贡献先线性筛出质数从小到大计算 每个数对倍数产生贡献贡献累加CODE#include<iostream> #include<cstdio>using namespace std;const int maxn=20000000;unsigned int seed,a[原创 2021-05-17 17:28:55 · 67 阅读 · 0 评论 -
洛谷P3383 【模板】线性筛素数
洛谷P3383 【模板】线性筛素数TITLE思路每个合数只被筛一次for(b[1]=1,i=2;i<=n;++i){ if(!b[i])z[++tot]=i;//标记质数 for(j=1;j<=tot&&i*z[j]<=n;++j)//用每个质数*i筛合数 { b[i*z[j]]=1;//标记为合数 if(!(i%z[j]))break;//如果(第j个质数)整除i,那么后面的数(=i*x)会被(第j个质数*x*(i/第j个质数))筛掉 }}原创 2021-05-17 17:21:26 · 99 阅读 · 0 评论 -
UVA12716 GCD等于XOR GCD XOR
UVA12716 GCD等于XOR GCD XORTITLE思路∵gcd(i,j)≤i−j≤ixorj\because \gcd(i,j) \leq i-j \leq i xor j∵gcd(i,j)≤i−j≤ixorj要gcd(i,j)=ixorjgcd(i,j)=ixorjgcd(i,j)=ixorj所以gcd(i,j)=i−j=ixorjgcd(i,j)=i-j=ixorjgcd(i,j)=i−j=ixorjgcd(i,j)=gcd(i,i−j)=i−jgcd(i,j)=gcd(i,i-原创 2021-05-17 16:59:44 · 86 阅读 · 0 评论