【正文】 mi[lc],mi[rc])。i++) scanf(%d, amp。 M。const int N=3000100。 //printf(mat[%d]get\n,i)。 if (!iamp。 while (inamp。 }else j=next[j]。void getNext(int m){ int i=0,j=1。char a[N],b[N]。}KMP 字符串匹配includecstdioincludecstringusing namespace std。i++){ st0+=st[i]。 while (~scanf(%s,st)){ string st0=$。 delete(p)。} } int ma=0。 else p[i]=1。 int mx=0,id=0。const int N=233333。 ++j) temp = x[j] * a[i][j]。 i=0。im 。j=n。 for (int i=k+1。continue。 for (int i=k+1。kmamp。int m,n。 if (k=n)return nk。++j)a[i][j]^=a[k][j]。 for (int i=k+1。 continue。 //r為第col列的一個(gè)1 for (int i=k。kmamp。 }}。jcol。 } int* operator [](int x){return a[x]。struct matrix{ int a[N][N]。高斯消元，判斷有無(wú)解的includecstdioincludecmathincludecstringincludeiostreamincludevectorusing namespace std。j=n+1。} return ret。 for (。 for (int i=0。k++) tmp[i][j]=(tmp[i][j]+a[i][k]*x[k][j])%P。j=n+1。} matrix operator * (matrix x){ matrix tmp 。 int row,col。const int P = 9973。代碼倉(cāng)庫(kù)目錄：01.【數(shù)學(xué)方法】矩陣快速冪02.【數(shù)學(xué)方法】高斯消元（na239。const int N=13。 matrix():row(N),col(N){memset(a,0,sizeof(a))。 for (int i=0。j++){ tmp[i][j]=0。 } return tmp。i=n+1。x。 } void print(){ for (int i=0。j++) printf(%d ,a[i][j])。typedef long long LL。 int row,col。} void print(){ for (int i=0。j++) printf(%d ,a[i][j])。int Gauss(matrix a,int m,int n){ int x_t = 0。amp。im。} if (r!=k)for (int i=col。im。 } for (int i=k。 //返回自由元個(gè)數(shù)}高斯消元，求出一組解的include iostreaminclude algorithminclude cstdioinclude cstringinclude cmathusing namespace std。double a[N][N],x[N]。amp。im。}//列全為0 if (r!=k)for(int i=col。im。j++)a[i][j]=a[k][j]*t。++i)//無(wú)解 if (fabs(a[i][n])EPS) return 1。 i){//回帶求解 double temp = a[i][n]。 x[i] = (temp / a[i][i])。//20W//在o(n)時(shí)間內(nèi)算出以每個(gè)點(diǎn)為中心的最大回文串長(zhǎng)度int Manacher(string st){ int len=()。 for (int i=1。 while (st[i+p[i]]==st[ip[i]])p[i]++。 for(int i=1。 return ma1。 for (int i=0。 st0+=。typedef long long ll。bool mat[N]。 next[0]=1。 }}void KMP(int n,int m){ memset(mat,0,sizeof(mat))。amp。amp。 j=next[j]。struct linetree{ define lc (t1) define rc (t1^1) int mi[N],M。 memset(mi,INF,sizeof(mi))。mi[i])。 } void change(int t,int x){ for (mi[t+=M]=x,t=1。 for (l+=M1,r+=M+1。 if ( ramp。int main(){ int n,q,ord,x,y。){ (n)。){ scanf(%d%d%d,amp。 if (ord)(x,y)。const int INF=0x3f3f3f3f。 node () { } node(int x,int y,int z){l=x。inline bool cmp1(node a,node b){return 。 while(Mn)M=1。 memset(ta, 0 ,sizeof(ta))。i++)l[i]=r[i]= iM 。 } inline void down(int t){ if (tM)return 。 ta[rc]=max(ta[rc],ta[t])。 ti[lc]=min(ti[lc],ti[t])。 mi[t]=min(mi[lc],mi[rc])。 ti[t]=min(ti[t],x)。return。 maintain(t)。 query(lc)。線段樹(shù)（區(qū)間加+賦值）includecstdioincludecstringincludeiostreamincludealgorithmusing namespace std。 inline void build(int n){ M=1。 memset(tag,0,sizeof(tag))。i++){//leaf if(i=n+M)scanf(%d,amp。 } for (int t=M。 len[t]=len[lc]+len[rc]。 sum[rc]=Set[t]*len[rc]。 tag[lc]=tag[rc]=0。 tag[rc]+=tag[t]。 } inline void _set(int t,int x){ sum[t]=x*len[t]。amp。 if (L=mid)change(lc,L,R,x)。amp。 if (L=mid)set(lc,L,R,x)。amp。 if (L=mid)ans+=query(lc,L,R)。SplayTreeincludecstdioincludealgorithmusing namespace std。class SplayTree{public: void Init(){rt=NULL。//if x has left child xf =yf。 xl=y。 xf = yf。 xr=y。 }else if (x==pl){ if (p==pfl){Zig(p)。} }else {//x==pr if (p==pfr){Zag(p)。} } } rt=x。} else if (xTkey)T=Tl。 while (T){ fa=T。//two the same keys } T=(Node*)malloc(sizeof(Node))。 if (fa){ if (fakeyx)fal=T。 if (NULL==T)return 。 Splay(T)。 Splay(T)。 return (Maxnum(T))。 } Node *Join(Node *t1,Node *t2){ if (NULL==t1)return t2。 return T。 t1=Tl。 printf(%d,Tkey)。 Delete(Tl)。//root}。//bf=balance factor,h=height Node *l,*r。} int BF(Node *l,Node *r){//get balance factor if (NULL==l amp。 else if (NULL == r) return lh。 ar=bl。 abf=BF(al,ar)。 b=al。 bh=max(H(bl),H(br)) + 1。 } Node *LRrorate(Node *a){//left then right al = Lrorate(al)。 } Node *RLrorate(Node *a){//right then left ar=Rrorate(ar)。 } void Insert(int x){_Insert(rt,x)。 Tbf=0。 } if (x Tx) _Insert(Tl,x)。//maintain Tbf=BF(Tl,Tr)。 else if (Tbf 0 amp。amp。 Trbf 0)T=RLrorate(T)。 int flag。 else if (large==INF)flag=0。 ans=(ans+xsmall)%BASE。small,int amp。 return Find(Tl,x,small,large)。T,int x){ if (NULL==T)return 。 else T=Lrorate(T)。 else T=Rrorate(T)。amp。 _Delete(Tl,tx)。//bf==0 or 1 } }else {//left || right Node *t=T。T=NULL。} void InOrder(Node *T){//print l rt r if (NULL==T)return 。 } void Free