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;
}