———数学知识———
Hacheylight
ZJ OIer
展开
-
UVA10140 Prime Distance (素数的判定) HQG_AC的博客
题意: 求l到r的素数中距离最小的是多少,最大的是多少,没有就输出 There are no adjacent primes.思路: 发现l和r很大,但是r-l的值不大。于是考虑缩短存储: prime[i]表示i+l是不是质数。 于是就只用打一波素数筛就可以了(可以不用线性筛也能过),之后模拟。注意:l,r都要long long 不然 像2147483647这样的数据就RE了。...原创 2018-05-21 16:04:37 · 180 阅读 · 0 评论 -
7.27日比赛 总结+题解
小T的GCD分数:10 预期分数:80~90(由于机子比较快,预期会高几分)题意 求GCD(al,al+1...,ar)==1GCD(al,al+1...,ar)==1GCD(a_l,a_{l+1}...,a_r)==1 的 max(r−l+1)max(r−l+1)max(r-l+1)求LCM(al,al+1...,ar)==∏rk=lakLCM(al,al+1...,ar)==∏...原创 2018-07-27 20:04:04 · 143 阅读 · 0 评论 -
数论模板总结
数论模板总结错误什么的都请大佬就出来,博主会飞速改!!!快速乘int fastmul(int a, int b, int p) { int x = 0 ; while (b) { if (b & 1) x = (x + a) % p ; a = (a + a) % p ; b >>= 1 ; } return x ;}原创 2018-10-29 21:18:19 · 205 阅读 · 4 评论 -
CF1097D Makoto and a Blackboard 题解
题意就是给一个数,每次可以把它替换成为它的一个因数,问替换k次后期望的数为多少一看就不会做,只好先从简单的开始着手如果一个数为质数,肯定非常好处理,有12k\frac{1}{2^k}2k1种可能为原数,其余都为1如果一个数的分解质因数形式为x=p1e1x=p_1^{e_1}x=p1e1怎么搞?很容易想出dpdpdpdp[i][j]dp[i][j]dp[i][j]表示做了iii次操作...原创 2019-01-05 20:44:32 · 552 阅读 · 0 评论 -
CF1062B Math 题解
本题居然没有题解,我来写一篇我们发现如果一个数xxx的分解质因数结果中有一些因数的质数不为零,我们总能够通过乘再开根使得他变为零,但不可能消失于是第一个答案就是该数的所有质因子的积yyy对于第二问,我们也就是要求从xxx到yyy的变化次数我们发现只会有一次乘法,因为只需要一次乘法+多次开根就可以搞定那我们只需要每次把yyy平方一下,如果是xxx的倍数的话就break特判正好等于xxx的...原创 2019-01-02 16:12:31 · 195 阅读 · 0 评论 -
数论初步总结
前面的简单小知识到时候再补充先写类欧原创 2019-01-27 09:33:32 · 379 阅读 · 2 评论 -
数论好题总结
。。。原创 2019-01-29 09:21:43 · 243 阅读 · 0 评论 -
0x37 ZAP (Mobius)
看着蓝书就入了 MobiusMobiusMobius 反演的坑了(但是蓝书并没有提到反演二字首先根据基本的数论常识,我们知道 ∑x=1a∑y=1b[gcd(x,y)=k]\sum\limits_{x=1}^a\sum\limits_{y=1}^b[\gcd(x,y)=k]x=1∑ay=1∑b[gcd(x,y)=k] 等价于 ∑x=1ak∑y=1bk[gcd(x,y)=1]\sum\lim...原创 2019-02-05 13:32:19 · 190 阅读 · 2 评论 -
Codeforces 1110C Meaningless Operations
没有摘要原创 2019-02-08 21:32:46 · 257 阅读 · 0 评论 -
数论ACM总结
莫名其妙发现有几道题目不是数论。。。原创 2019-01-30 12:58:08 · 314 阅读 · 0 评论 -
7.25 题目总结 +题解
考题1传送门小T的矩阵 得分:303030 预期:100100100题意:sumi=xor∑nj=1aijsumi=xor∑j=1naijsum_i=xor\sum_{j=1}^na_{ij} 第i列ans=xor∑ni=1sumians=xor∑i=1nsumians=xor \sum_{i=1}^nsum_i 学了数论,看这道题目就特别的像数论,但并不是数论。。。...原创 2018-07-27 16:37:57 · 173 阅读 · 0 评论 -
【模板】乘法逆元
逆元有4种求法:费马小定理exgcd递推式(本方法)快速幂// luogu-judger-enable-o2#include <bits/stdc++.h>using namespace std ;typedef long long ll ;const int N = 3000010 ;ll n,p ;ll a[N] ;int main(){ ...原创 2018-07-26 21:07:30 · 648 阅读 · 0 评论 -
【NOI2018】屠龙勇士
我写的是85分的算法。至于我怎么过的,你一看就知道了。请不要看本蒟蒻的Code了。#include &lt;bits/stdc++.h&gt;using namespace std ;typedef long long ll ;const int N = 100010 ;multiset &lt;ll&gt; s ;multiset &lt;ll&gt; ::iterato...原创 2018-07-25 08:19:57 · 334 阅读 · 0 评论 -
Map Generator POJ - 3557 DP+概率
给你n个点,生成这条边的概率位P,判断联通的概率是多少?设连通块的大小为k,成功的概率不好枚举,但可以枚举不成功的概率设F[n]=1−G[n],G[n]F[n]=1−G[n],G[n]F[n]=1-G[n],G[n]表示成功的概率G[i]+=c[i−1][j−1]∗F[j]∗pow(1−p,double(j∗(i−j)))G[i]+=c[i−1][j−1]∗F[j]∗pow(1−p,d...原创 2018-07-21 16:15:34 · 183 阅读 · 0 评论 -
青蛙的约会 (同余,线性同余方程) HQG_AC的博客
模拟样例易知 如果青蛙们走k步,则必有km+x≡kn+y(modkm+x≡kn+y(modkm+x ≡kn+y (mod L)L)L)将该式化简,−&gt;(m−n)∗k≡y−x(mod−&gt;(m−n)∗k≡y−x(mod->(m-n)*k ≡y-x (mod L)L)L)−&gt;ak≡b(modL)−&gt;ak≡b(modL)->ak ≡b (mod L)之后参见NO...原创 2018-07-26 19:28:29 · 141 阅读 · 0 评论 -
上帝与集合的正确用法 (欧拉函数,快速幂)HQG_AC
这个欧拉函数裸题先用筛法(埃氏筛还是欧拉筛都可以,反正我用了欧拉筛,因为快)然后运用欧拉函数积性的性质,solve就可以了,至于求解,快速幂帮您搞定#include <bits/stdc++.h>using namespace std ;const int N = 10000010 ;int phi[N],p[N] ;bool flag[N] ;int T...原创 2018-07-26 19:49:13 · 167 阅读 · 0 评论 -
Biorhythms POJ 1006 uva 756 (同余) HQG_AC
实际上这是一个同余问题。可以应用中国剩余定理。我是直接暴力求解。注意格式#include<bits/stdc++.h>using namespace std;int main(){ int p,e,i,d,count,x,n; scanf("%d",&n); while(n--){ count=0 ; wh...原创 2018-07-26 19:54:45 · 233 阅读 · 0 评论 -
余数之和 (暴力,思维,前缀和) HQG_AC
#include <bits/stdc++.h>using namespace std;typedef long long ll ;ll ans,k,n,check;ll f(ll l,ll r){ l--; return (r+1)*r/2-(l+1)*l/2;}int main(){ scanf("%lld%lld",&n,&...原创 2018-07-26 19:58:25 · 348 阅读 · 0 评论 -
计算器(BSGS,快速幂,exgcd)
这个题目坑啊,查来查去,排查了好久,发现自己快速幂写错了。。。这个题目有三个询问1、yzyzy^z modmodmod ppp2、xyxyxy ≡ z(modz(modz(mod p)p)p) (求x)3、yx≡z(modyx≡z(mody^x ≡z(mod p)p)p)(求x)第一个用快速幂求。第二个就是裸的ecgcd,见 同余方程前两个是35分第三个就是裸的BS...原创 2018-07-26 20:10:15 · 618 阅读 · 0 评论 -
仪仗队 (欧拉函数) HQG_aC
我们可以将图沿对角线切成两个直角三角形观察可以发现,从第三行开始,第i行中当前直角三角形范围内能观察到的点数等于(i-1)的欧拉函数值// luogu-judger-enable-o2#include<bits/stdc++.h>using namespace std;int a[40001],n;int fact(int n){ if(a[n]>=0...原创 2018-07-26 20:15:10 · 319 阅读 · 0 评论 -
沙拉公主的困惑 (乘法逆元+素数判定)
加了各种优化,BZOJ还是没过,洛谷倒是过了主要预处理不是O(n),所以挂了。乘法逆元+素数#pragma GCC diagnostic error "-std=c++11"#pragma GCC optimize("-fdelete-null-pointer-checks,inline-functions-called-once,-funsafe-loop-optimizatio...原创 2018-07-26 20:18:47 · 224 阅读 · 0 评论 -
同余方程
直接exgcd#include <bits/stdc++.h>using namespace std ;typedef long long ll ;int a,b,x,y ;void exgcd(int a,int b,int &x,int &y){ if (b==0){ x=1 ; y=0 ; }...原创 2018-07-26 20:25:07 · 133 阅读 · 0 评论 -
考试总结
膜你赛考的好差啊原创 2019-02-19 10:59:10 · 319 阅读 · 0 评论