时光真疯狂, 我一路执迷于匆忙.

自然数幂和模板（第一类斯特林数）ll get(ll n,ll k){ fo(i,0,k) su[i][0]=0,ss[i][i]=1; fo(i,1,k) fo(j,1,i-1) ss[i][j]=(ss[i-1][j-1]-(i-1)*ss[i-1][j]%mo)%mo; f[0]=(n+1)%mo;f[1]=n*(n+1)%mo;

int extend(int x,int p) { int np=++tot;sam[np].len=sam[p].len+1; while (p&&!sam[p].son[x]) sam[p].son[x]=np,p=sam[p].pr; if (!p) sam[np].pr=1;else { int q=sam[p].son[x]; if

accessvoid access(int x) { int y=0; while (x) { splay(x,0); f[t[x][1]]=0;p[t[x][1]]=x; t[x][1]=y;f[y]=x;p[y]=0; updata(x); y=x;x=p[x]; }}makerootvoi

带预处理的，听说精度会炸=w=void DFT(complex *a,int flag) { for(int i=0;i<len;i++) { int p=0; for(int j=i,k=0;k<lg;j/=2,k++) p=(p<<1)+(j&1); t[p]=a[i]; } for(int m=2;m<=len;m=

学习UOJ上dalao们的初始化黑科技

发现我之前写的Pollard_rho模板是假的。。。 学(抄)了一发真的，但跑的贼慢ll mult(ll x,ll y,ll Mo) { x%=Mo;y%=Mo; ll tmp=(ll)((long double)x*y/Mo+1e-8)*Mo; return (x*y-tmp+Mo)%Mo;}ll pwr(ll x,ll y,ll Mo) { ll z=1;x%

还是用0号点比较好，用数据结构维护的时候多注意一下。

Codevec记录所有在虚树中的点，方便删除 add1只需要多开一个lst就可以实现O(k)删除 为了方便默认1号点为根void build() { sort(a+1,a+tot+1,cmp); st[top=1]=1;vec.push_back(1); fo(i,1,tot) { int lca=get_lca(a[i],st[top]);

想写一个多项式全家桶就开了这个坑 注意所有运算均是在模x^2n的域中进行的Part A：多项式ln，exp，求幂 多项式ln：假设我们要求ln(f(x))ln⁡(f(x))\ln(f(x)) 我们不妨求导之后再积分，那么就是∫f′(x)f(x)∫f′(x)f(x)\int {f'(x)\over f(x)} 只需要求逆多项式exp： 先介绍牛顿迭代 假设我们知道了一个函数g(x...

不能读负数，不能各种读神奇的格式。 那你摆出来干什么，赚访问量？int get() { char ch;while (!isdigit(ch=getchar())); int o=ch-48;while (isdigit(ch=getchar())) o=o*10+ch-48; return o;}哎呀，被法线啦(>_<)

转自WerKeyTom_FTD的blog 1、整体二分及cdq分治实现时，每个区间不需要单独开队列。可以把操作弄到一个数组了，然后多两个参u、v表示这个区间的操作在u~v，做完后对每个操作打标记表示是否往右区间传即可。 2、两个long long相乘取模怎么办？ 第一种是快速乘，类似快速幂，复杂度多了个log 第二种是找一个比较小的数，然后搞一波，详见代码ll qsc(ll x,ll y){

splay旋转模板void rotate(int x) { int y=f[x],z=son(x);f[x]=f[y]; if (f[x]) t[f[x]][son(y)]=x; t[y][z]=t[x][1-z]; if (t[x][1-z]) f[t[x][1-z]]=y; f[y]=x;t[x][1-z]=y;}void splay(int x,in

忘了的时候来抄模板（我记性差~）scanf("%d%d",&n,&m);if (n>m) swap(n,m);mu[1]=1;fo(i,2,n) { if (!f[i]) p[++tot]=i,mu[i]=-1; fo(j,1,tot) { int k=p[j]*i;if (k>n) break;f[k]=1; if (!(i%p[j])) break;m

组合数取模模板（Pascal（到时候改C++））function fact(x,y:longint):note;var k:note;begin if x<a[y] then begin fact.e:=sum[x,y];fact.f:=0;exit; end; fact.f:=x div a[y];k:=fact(x div a[y],y

高斯消元模板void gauss() { fo(i,1,n-1) { fo(j,i+1,n) if (abs(matrix[j].a[i])>abs(matrix[i].a[i])) swap(matrix[i],matrix[j]); fo(j,i+1,n) { db f=matrix[j].a[i]/matrix[

sa模板（pascal倍增）for i:=1 to n do begin s[i]:=t[i];x[i]:=ord(s[i]); if x[i]>max then max:=x[i]; inc(ws[x[i]]); end; for i:=1 to max do inc(ws[i],ws[i-1]); for i:=n

离散化模板（看不懂%%%）void prepare(int *x) { fo(i,1,n) data[i]=x[i]; sort(data+1,data+n+1); int m=unique(data+1,data+n+1)-data-1; fo(i,1,n) x[i]=lower_bound(data+1,data+m+1,x[i])-data;}

树链剖分（两遍dfs）void dfs(int x,int y) { size[x]=1;fa[x]=y;d[x]=d[y]+1;int k=0; rep(i,x) if (t[i]!=y) { dfs(t[i],x);size[x]+=size[t[i]]; if (size[t[i]]>k) k=size[t[i]],son[x]=t[i];

ac自动机模板（makefail）void makefail() { for(int i=0,j=1;i<j;) { int x=d[++i]; fo(k,0,25) if (trie[x].son[k]) { int y=trie[x].son[k]; if (x) {

sa模板（c++倍增）void tsort() { memset(ws,0,sizeof(ws));int mx=0; fo(i,1,n) ws[x[y[i]]]++,mx=max(mx,x[y[i]]); fo(i,1,mx) ws[i]+=ws[i-1]; fd(i,n,1) sa[ws[x[y[i]]]]=y[i],ws[x[y[i]]]--;}void g

我之前可能学了假的BMnamespace BM{ vector<int> h[N]; int cnt,fail[N],d[N],mx; vector<int> work(int n,int *a) { h[cnt=mx=0].clear(); fo(i,1,n) { int now=-a[i]; for(unsigned j=0;j<h[c...