模板合集(未完

dijkstra

priority_queue<pii,vector<pii>,greater<pii> >q;
void dijkstra(int s){
	memset(dis,inf,sizeof(dis));
    memset(vis,0,sizeof(vis));
    q.push(make_pair(dis[s]=0,s));
    while(!q.empty()){
		int u=q.top().second;q.pop();
        if(vis[u]) continue;vis[u]=1;
        for(int i=head[u],v;i;i=e[i].nxt)
            if(dis[v=e[i].v]>dis[u]+e[i].w) dis[v]=dis[u]+e[i].w,q.push(make_pair(dis[v],v));
    }
}

prim

priority_queue<pii,vector<pii>,greater<pii> >q;
void prim(){
	memset(dis,inf,sizeof(dis));
    memset(vis,0,sizeof(vis));
    q.push(make_pair(dis[1]=0,1));
    while(!q.empty()){
		int u=q.top().second;q.pop();
        if(vis[u]) continue;
        ans+=dis[u],vis[u]=1;
        for(int i=head[u],v,w;i;i=e[i].nxt)
            if(!vis[v=e[i].v]&&dis[v]>(w=e[i].w)) q.push(make_pair(dis[v]=w,v));
    }
}

kruskal

void kruskal(){
    sort(e+1,e+m+1);
	for(int i=1;i<=n;++i) f[i]=i;
    for(int i=1,u,v,cnt=0;i<=m&&cnt<n-1;++i)
        if(find(u=e[i].u)!=find(v=e[i].v)) f[f[u]]=f[v],ans+=e[i].w,++cnt;
}

tarjan有向图缩点

void tarjan(int u){
	dfn[u]=low[u]=++idx,s.push(u),inst[u]=1;
    for(int i=head[u],v;i;i=e[i].nxt)
		if(!dfn[v=e[i].v]) tarjan(v),low[u]=min(low[u],low[v]);
        else if(inst[v]&&dfn[v]<low[u]) low[u]=dfn[v];
    if(dfn[u]==low[u]){
        int v;++Bcnt;
        do{
            v=s.top(),s.pop();
            bl[v]=Bcnt,++sz[Bcnt],inst[v]=0;
        }while(u!=v);
	}
}

2-SAT

线性筛

void prime(int N){
	memset(prime,0,sizeof(prime));
	memset(v,0,sizeof(v));
	for(int i=2;i<=n;++i){
		if(!v[i]) v[i]=i,prime[++num_prime]=i;
		for(int j=1;j<=num_prime&&i*prime[j]<=n;++j){
			v[i*prime[j]]=prime[j];
			if(!(i%prime[j])) break;
		}
	}
}

扩欧

求关于x的同余方程(axequiv 1(mod b)) 的最小正整数解

void exgcd(ll a,ll b,ll &x,ll &y){
    if(!b){x=1,y=0;return;}
    exgcd(b,a%b,x,y);
    ll t=x;x=y,y=t-(a/b)*y;
}

逆元

线性筛逆元

void Qinv(int P){
    inv[1]=1;
    for(int i=2;i<=n;++i) inv[i]=(ll)inv[P%i]*(P-P/i);
}

费马小定理求逆元

中国剩余定理

void exgcd(ll a,ll b,ll &x,ll &y){
    if(!b) {x=1,y=0;return;}
    exgcd(b,a%b,x,y;
    ll t=x;x=y,y=t-a/b*y;
}
ll china(int n,int *a,int *m){
    ll ans=0,M=1ll,x,y;
    for(int i=1;i<=n;++i) M*=m[i];
    for(int i=1;i<=n;++i){
        ll mi=M/m[i];
        exgcd(mi,m[i],x,y);
        x=(x%m[i]+m[i])%m[i];
        ans=(ans+a[i]*mi*x)%M;
	}
    return (ans+M)%M;
}

扩展中国剩余定理

BSGS

求一个(x)使得(y^xequiv z(mod p))

int qpow(int a,int b){
	int res=1;
    while(b){
        id(b&1) res=(ll)res*a%p;
        a=(ll)a*a%p,b>>=1;
    }
    return res;
}
map<int,int>hash;
int BSGS(){
    hash.clear(),y%=p,z%=p;
    int t=(int)sqrt(p)+1;
    for(int i=0,val;i<t;++i)
    	val=(ll)z*qpow(y,i)%p,hash[val]=i;
    y=qpow(y,t);
    if(!y) return !z?1:-1;
    for(int i=0,val,j;i<=t;++i){
        val=qpow(y,i);
        j=hash.find(val)==hash.end()?-1:hash[val];
        if(j>=0&&i*t-j>=0) return i*t-j;
    }
    return -1;
}

线段树

void pushup(int o){
	sum[o]=(sum[ls]+sum[rs])%p;
}
void updatenode1(int o,int l,int r,int k){
	sum[o]=1ll*sum[o]*k%p;
	add[o]=1ll*add[o]*k%p;
	mul[o]=1ll*mul[o]*k%p;
}
void updatenode2(int o,int l,int r,int k){
	sum[o]=(sum[o]+1ll*(r-l+1)*k%p)%p;
	add[o]=(add[o]+k)%p;
}

void pushdown(int o,int l,int r){
	int mid=(l+r)>>1;
	if(mul[o]!=1){
		updatenode1(ls,l,mid,mul[o]);
		updatenode1(rs,mid+1,r,mul[o]);
		mul[o]=1;
	}
	if(add[o]){
		updatenode2(ls,l,mid,add[o]);
		updatenode2(rs,mid+1,r,add[o]);
		add[o]=0; 
	}
}

void update1(int o,int l,int r,int x,int y,int k){
	if(r<x||y<l) return;
	if(x<=l&&r<=y) {updatenode1(o,l,r,k);return;}
	pushdown(o,l,r);
	int mid=(l+r)>>1;
	update1(ls,l,mid,x,y,k);
	update1(rs,mid+1,r,x,y,k);
	pushup(o);
	return;
}
void update2(int o,int l,int r,int x,int y,int k){
	if(r<x||y<l) return;
	if(x<=l&&r<=y) {updatenode2(o,l,r,k);return;}
	pushdown(o,l,r);
	int mid=(l+r)>>1;
	update2(ls,l,mid,x,y,k);
	update2(rs,mid+1,r,x,y,k);
	pushup(o);
	return;
}

void buildtree(int o,int l,int r){
	add[o]=0,mul[o]=1;
	if(l==r) {sum[o]=a[l];return;}
	int mid=(l+r)>>1;
	buildtree(ls,l,mid);
	buildtree(rs,mid+1,r);
	pushup(o);
}
int query(int o,int l,int r,int x,int y){
	int ans=0;
	if(r<x||y<l) return 0;
	if(x<=l&&r<=y) return sum[o];
	pushdown(o,l,r);
	int mid=(l+r)>>1;
	ans=(ans+query(ls,l,mid,x,y))%p;
	ans=(ans+query(rs,mid+1,r,x,y))%p;
	return ans;
}

int main(){
	rd(n),rd(m),rd(p);
	for(rg int i=1;i<=n;++i) rd(a[i]);
	buildtree(1,1,n);
	for(rg int i=1,op,x,y,k;i<=m;++i){
		rd(op),rd(x),rd(y);
		if(op!=3) rd(k);
		if(op==1) update1(1,1,n,x,y,k);
		else if(op==2) update2(1,1,n,x,y,k);
		else printf("%d
",query(1,1,n,x,y));
	}
	return 0;
}