4 wspl654321

尚未进行身份认证

我要认证

。。。。。。

等级
TA的排名 2w+

vim使用

http://www.xuebuyuan.com/1116162.html https://blog.csdn.net/luojj26/article/details/50935921 https://blog.csdn.net/lhy2932226314/article/details/69668891# set nu “nu=number set cul “行亮 cuc 列 s...

2019-06-12 21:29:41

NTT

#include<cstdio>#include<iostream>using namespace std;const int P=998244353,M=510000;int ans=0,g[M],inv[M],n,m,w[2][M],jc[M];int abs(int x){return x>0?x:-x;}int pow(int a,int k){i...

2019-06-12 21:29:15

FFT

#include<cstdio>#include<iostream>#include<cmath>using namespace std;const double PI=acos(-1);const int M=110000; struct complex{//虚数a+bi double a,b;complex(){}complex(doub...

2019-06-12 21:28:53

鬼东西

fft 点分治 李超线段树 cdq分治 后缀数组 分块大发 多项式除法,求逆 ,提答构造 后缀自动机 https://blog.csdn.net/clover_hxy/article/details/68059043硬币 http://www.cnblogs.com/CQzhangyu/p/7054998.html 有趣 https://blog.sengxia...

2019-06-12 21:28:34

莫比乌斯反演与杜教筛

莫比乌斯反演我们知道积性函数可以迪雷克卷积 我们知道F(x)=Σd|xf(d) 要求f(x) 两面卷上一个μ 得 F∗μ=f∗1∗μF∗μ=f∗1∗μF*μ=f*1*μ 得 F∗μ=f∗(1∗μ)F∗μ=f∗(1∗μ)F*μ=f*(1*μ) 得 F∗μ=f∗eF∗μ=f∗eF*μ=f*e 得F∗μ=fF∗μ=fF*μ=f 也就是f=Σd|xμ(x/d)f(d)f=Σd|xμ(x...

2019-06-12 21:27:59

P4022 [CTSC2012]熟悉的文章

这真是个神dp 我肯定想不出来,太牛逼了题解假设我们知道了len[i]为i之前的最长匹配熟悉子串 我们可以二分答案,得到一个L 那么就可以得到 dp[i]表示匹配到i的最长熟悉长度 dp[i]=max(dp[j]+i−j)(i−len[i]<=j<=i−L)dp[i]=max(dp[j]+i−j)(i−len[i]<=j<=i−L)dp[i]=max(d...

2019-06-12 21:26:44

杂题

HH去散步 矩阵乘法优化dp 对边dp dp[i]=西格玛dp[j]; 最后找是连接节点的边,计数// luogu-judger-enable-o2#include<cstdio>#include<cstring>using namespace std;const int M=121,p=45989; int u[M],v[M],n,m,t,S,T,st...

2019-06-12 21:26:06

扩展欧拉定理

ab≡x(modm)求xab≡x(modm)求x a^b≡x(mod m) 求xgcd(a,m)=1,欧拉定理:ab≡a的(bmodφ(m))次方(modm)gcd(a,m)=1,欧拉定理:ab≡a的(bmodφ(m))次方(modm)gcd(a,m)=1,欧拉定理:a ^ b ≡a 的( b mod φ ( m ) )次方 (mod m)gcd(a,m)>1,且b&gt...

2019-06-12 21:25:50

范德蒙矩阵

https://wenku.baidu.com/view/70fe14f8f705cc1755270933.html 范德蒙矩阵的行列式值 =所有可能差值的积 | 1 , 1 , 1 , 1 , 1 ||x1,x2,x3,x4,x5||x1^2,x2^2,x3^2,x4^2,x5^2||x1^3,x2^3,x3^3,x4^3,x5^3||x1^4,x2^...

2019-06-12 21:24:31

更牛逼的高斯消元(辗转相除法)

没有精度问题,但多了个log,不明白的话手动模拟一下for(int i=1;i<n;i++) { for(int j=i+1;j<n;j++) while(a[j][i]) { int t=a[i][i]/a[j][i]; for(int k=i;k<n...

2019-06-12 21:23:36

最小顶标和与最大权值匹配

结论:不知为什么。 满足任意xi+xj>=cost[i][j] xi到xj连cost[i][j]的边 求出最大权值匹配就是最小顶标和 最大权职匹配可以用最大费用可行流

2019-06-12 21:22:17

KD_TREE总结

用途KD_TREE可以求K维空间最短/长两点的距离 或者前T长的距离,蒟蒻大概就知道这一点。。。思想把空间切割,完成缩短时间的半目的实现(建树,插入)怎样切割??? 用二叉树来实现,每一层按层数%K这一维切割 记录每一个子空间的状态就可以了,例如处理出每个空间每维的上限下限 然后合并子树的状态,插入类似复杂度就可以做到随机数据logN的复杂度,但毒瘤出题...

2019-06-12 21:21:24

P4211 [LNOI2014]LCA

又一鬼题 对于每次查询 把l-r到跟的路上+1 然后深度就是这个位置的全职 然后查分一下,查询r-l 这样就可以离线记下查询位置 只添加,不用删除 就可以利用之前的信息了// luogu-judger-enable-o2#include<cstdio>#include<iostream>#include<cstring>#inclu...

2018-04-01 09:21:45

线性基

P3292 [SCOI2016]幸运数字 亦或和 最大值就是二维线性基了 倍增维护线性基 本想写树剖维护,想了想还是倍增简单。。。// luogu-judger-enable-o2#include<cstdio>#include<cstring>#include<iostream>#define ll long longusing name...

2018-03-31 18:59:10

灾难

搞出灭绝树。 一群动物的祖先就是他们的灭绝祖先 倍增搞lca// luogu-judger-enable-o2#include<cstdio>#include<iostream>#include<cstring>#include<queue>using namespace std;const int M=110000;int n...

2018-03-31 18:48:06

栗栗的书架

本以为是线段树套主席书 不过看数据做题,大数据只有一条线。。 二分一个厚度直接大数据主席书 小数据暴力dp// luogu-judger-enable-o2// luogu-judger-enable-o2// luogu-judger-enable-o2#include<cstdio>#include<iostream>#define mid (l+r...

2018-03-31 18:44:54

费用刘

二分一个费用 限制上限,跑最大流// luogu-judger-enable-o2#include<cstdio>#define ll long long#include<cstring> #include<queue>using namespace std;const int M=10000,inf=0x7fffffff;const l...

2018-03-31 18:41:03

杂题

P1955 [NOI2015]程序自动分析 离散化+并查集判约束 可以先弄0的再弄1 避开了一些情况 一开始写了个反集,不知为啥错了。。。// luogu-judger-enable-o2#include<cstdio>#include<iostream>#include<algorithm>#define ll long long u...

2018-03-31 18:39:52

三分法

维护凸函数 可以用导数二分水做。 但某些不知通项的函数只能拉格朗日插值做了。。 不过我不会。。。。 但不是多项式的也无法作,还是三分吧。。。 宅男计划 猜出是一个凸函数 三分+ 神TM贪心,不懂就这么写了。。。// luogu-judger-enable-o2// luogu-judger-enable-o2#include<cstdio>#include...

2018-03-31 17:41:09

小Z的袜子

简单莫队 对袜子分块,对操作排序 概率是 某颜色个数/总个数 然后暴力修改就行了// luogu-judger-enable-o2#include<cstdio>#include<iostream>#include<cmath>#include<algorithm> #define ll long longusing names...

2018-03-31 17:25:16

查看更多

勋章 我的勋章
  • 持之以恒
    持之以恒
    授予每个自然月内发布4篇或4篇以上原创或翻译IT博文的用户。不积跬步无以至千里,不积小流无以成江海,程序人生的精彩需要坚持不懈地积累!
  • 勤写标兵Lv3
    勤写标兵Lv3
    授予每个自然周发布7篇到8篇原创IT博文的用户。本勋章将于次周周三上午根据用户上周的博文发布情况由系统自动颁发。