快读+快写+__int128
//不用__int128的话框架是一样的
#include<stdio.h>
#include<iostream>
using namespace std;
void read(__int128 &x)
{
x=0;
int f=1;char ch=getchar();
while(ch<'0'||ch>'9')
{
if(ch=='-') f=-1;
ch=getchar();
}
while('0'<=ch&&ch<='9')
{
x=x*10+ch-'0';//x=(x<<3)+(x<<1)+ch-'0';
ch=getchar();
}
}
void write(__int128 x)
{
if(x<0) x=-x,putchar('-');
if(x>9) write(x/10);
putchar(x%10+'0');
}
int main()
{
__int128 a,b;
read(a);
read(b);
write(a*b);
}
前缀和
#include<stdio.h>
#include<iostream>
using namespace std;
typedef long long lll;
int n,m;
int s[5010][5010];
int main()
{
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)
{
int x,y,v;
scanf("%d%d%d",&x,&y,&v);
x+=1;y+=1;
s[x][y]=v;
/*s[x+m-1+1][y+m-1+1]+=v;
s[x+m-1+1][y]-=v;
s[x][y+m-1+1]-=v;*/
}
for(int i=2;i<=5001;i++)
{
s[1][i]+=s[1][i-1];
s[i][1]+=s[i-1][i];
}
for(int i=2;i<=5005;i++)
{
for(int j=2;j<=5005;j++)
{
s[i][j]=s[i-1][j]+s[i][j-1]-s[i-1][j-1]+s[i][j];
}
}
int ans=-1;
for(int i=m;i<=5005;i++)
{
for(int j=m;j<=5005;j++)
{
int t=s[i][j]-s[i-m][j]-s[i][j-m]+s[i-m][j-m];
ans=max(ans,t);
}
}
printf("%d",ans);
}
线段树
线段树1(单点修改)
#include<stdio.h>
#include<iostream>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
int n,m;
int a[100005];
struct p{
ll l,r,sum;
int f;
}tree[100005*4];
void build(int k,int l,int r)
{
tree[k].f=0;
tree[k].l=l;tree[k].r=r;
if(l==r)
{
tree[k].sum=a[l];
return;
}
int mid=(l+r)>>1;
build(k<<1,l,mid);
build(k<<1|1,mid+1,r);
tree[k].sum=tree[k<<1].sum+tree[k<<1|1].sum;
}
void push(int k,int l,int r)
{
if(tree[k].f!=0)
{
int mid=(l+r)>>1;
int s=tree[k].f;
tree[k<<1].f+=s;tree[k<<1|1].f+=s;
tree[k<<1].sum+=(mid-l+1)*s;tree[k<<1|1].sum+=(r-mid)*s;
tree[k].f=0;
}
}
void add(int k,int l,int r,int x)
{
if(l<=tree[k].l&&tree[k].r<=r)
{
tree[k].sum+=(tree[k].r-tree[k].l+1)*x;
tree[k].f+=x;
return;
}
if(tree[k].l>r||tree[k].r<l) return;
push(k,tree[k].l,tree[k].r);
int mid=(tree[k].l+tree[k].r)>>1;
if(l<=mid) add(k<<1,l,r,x);
if(r>mid) add(k<<1|1,l,r,x);
//add(k<<1,l,r,x);add(k<<1|1,l,r,x);
tree[k].sum=tree[k<<1].sum+tree[k<<1|1].sum;
}
ll ask(int k,int l,int r)
{
if(l<=tree[k].l&&tree[k].r<=r)
{
return tree[k].sum;
}
if(tree[k].l>r||tree[k].r<l) return 0;
push(k,tree[k].l,tree[k].r);
int mid=(tree[k].l+tree[k].r)>>1;ll ans=0;
if(l<=mid) ans+=ask(k<<1,l,r);
if(r>mid) ans+=ask(k<<1|1,l,r);
return ans;
}
void print()
{
printf("\n");
for(int i=1;i<=n;i++)
{
printf("%d ",ask(1,i,i));
}
printf("\n\n");
}
int main()
{
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++) scanf("%d",&a[i]);
build(1,1,n);
while(m--)
{
int t,x,y,k;
scanf("%d",&t);
if(t==1)
{
scanf("%d%d%d",&x,&y,&k);
add(1,x,y,k);
//print();
}
else
{
scanf("%d%d",&x,&y);
printf("%lld\n",ask(1,x,y));
//print();
}
}
}
线段树2(区间修改)
#include<stdio.h>
#include<iostream>
#include<string.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
struct p{
int l,r;ull sum;
ull add;ull mul;
}tree[100005*4];
ull a[100005];
int mod;int n,m,p;
void build(int k,int l,int r)
{
tree[k].add=0;tree[k].mul=1;
tree[k].l=l;tree[k].r=r;
if(l==r)
{
tree[k].sum=a[l];
return ;
}
else
{
int mid=(l+r)>>1;
build(k<<1,l,mid);
build(k<<1|1,mid+1,r);
tree[k].sum=tree[k<<1].sum+tree[k<<1|1].sum;
}
tree[k].sum%=mod;
}
void push(int k,int l,int r)
{
int mid=(l+r)>>1;
tree[k<<1].sum=(tree[k<<1].sum*tree[k].mul+tree[k].add*(mid-l+1))%mod;
tree[k<<1|1].sum=(tree[k<<1|1].sum*tree[k].mul+tree[k].add*(r-mid))%mod;
tree[k<<1].mul=(tree[k<<1].mul*tree[k].mul)%mod;
tree[k<<1|1].mul=(tree[k<<1|1].mul*tree[k].mul)%mod;
tree[k<<1].add=(tree[k<<1].add*tree[k].mul+tree[k].add)%mod;
tree[k<<1|1].add=(tree[k<<1|1].add*tree[k].mul+tree[k].add)%mod;
tree[k].mul=1;tree[k].add=0;
return;
}
void add(int k,int l,int r,int x)
{
if(l<=tree[k].l&&tree[k].r<=r)
{
tree[k].add=(tree[k].add+x)%mod;
tree[k].sum=(tree[k].sum+x*(tree[k].r-tree[k].l+1))%mod;
return;
}
if(r<tree[k].l||tree[k].r<l) return;
push(k,tree[k].l,tree[k].r);
int mid=(tree[k].l+tree[k].r)>>1;
if(l<=mid) add(k<<1,l,r,x);
if(r>mid) add(k<<1|1,l,r,x);
tree[k].sum=(tree[k<<1].sum+tree[k<<1|1].sum)%mod;
}
void mul(int k,int l,int r,int x)
{
if(l<=tree[k].l&&tree[k].r<=r)
{
tree[k].sum=(tree[k].sum*x)%mod;
tree[k].mul=(tree[k].mul*x)%mod;
tree[k].add=(tree[k].add*x)%mod;
return;
}
if(tree[k].l>r||tree[k].r<l) return;
push(k,tree[k].l,tree[k].r);
int mid=(tree[k].l+tree[k].r)>>1;
if(l<=mid) mul(k<<1,l,r,x);
if(r>mid) mul(k<<1|1,l,r,x);
tree[k].sum=(tree[k<<1].sum+tree[k<<1|1].sum)%mod;
}
ull query(int k,int l,int r)
{
if(l<=tree[k].l&&tree[k].r<=r)
{
return tree[k].sum%mod;
}
if(tree[k].l>r||tree[k].r<l) return 0;
push(k,tree[k].l,tree[k].r);
ull ans=0;
int mid=(tree[k].l+tree[k].r)>>1;
if(l<=mid) ans=(ans+query(k<<1,l,r))%mod;
if(r>mid) ans=(ans+query(k<<1|1,l,r))%mod;
return ans%mod;
}
void check()
{
cout<<endl;
for(int i=1;i<=n;i++)
{
cout<<query(1,i,i)<<" ";
}
cout<<endl<<endl;
}
int main()
{
//scanf("%d%d%d",&n,&m,&p);
ios::sync_with_stdio(false);
cin>>n>>m>>p;mod=p;
for(int i=1;i<=n;i++)cin>>a[i];
build(1,1,n);
for(int i=1;i<=m;i++)
{
int t,x,y,k;
cin>>t;
if(t==1)
{
cin>>x>>y>>k;
mul(1,x,y,k);
//check();
}
else if(t==2)
{
cin>>x>>y>>k;
add(1,x,y,k);
//check();
}
else
{
cin>>x>>y;
cout<<query(1,x,y)%mod<<endl;
//check();
}
}
}
可持续化线段树(懒人版)
#include<stdio.h>
#include<iostream>
#include<ext/rope>
using namespace std;
using namespace __gnu_cxx;
const int N=1e6+6;
rope<int> *s[N];
int n,m;
int a[N];
int main()
{
scanf("%d%d",&n,&m);
s[0]=new rope<int>();
s[0]->append(0);
for(int i=1;i<=n;i++)
{
scanf("%d",&a[i]);
}
s[0]->insert(1,a+1,n);
for(int i=1;i<=m;i++)
{
int v,k,a,b;
scanf("%d%d",&v,&k);
s[i]= new rope<int>(*s[v]);
if(k==1)
{
scanf("%d%d",&a,&b);
s[i]->replace(a,b);
}
else
{
scanf("%d",&a);
printf("%d\n",s[i]->at(a));
}
}
}
LCA
LCA(倍增)
#include<stdio.h>
#include<iostream>
#include<algorithm>
using namespace std;
const int N=5e5+5;
struct p{
int to,fr;
}all[N*4];
int g[N];int total=1;
int n,m,s;
int vis[N];
int st[N][32];
int dep[N];
void dfs(int x)
{
for(int i=g[x];i;i=all[i].fr)
{
int to=all[i].to;
if(vis[to]==0)
{
vis[to]=1;
dep[to]=dep[x]+1;
st[to][0]=x;
dfs(to);
}
}
}
int lca(int x,int y)
{
if(dep[x]<dep[y]) swap(x,y);
int maxdep=0;
while((1<<maxdep)<=dep[x]) maxdep++;
for(int i=maxdep;i>=0;i--)
{
if(dep[x]-(1<<i)>=dep[y])
{
x=st[x][i];
}
}
if(x==y) return x;
for(int i=maxdep;i>=0;i--)
{
if(st[x][i]!=st[y][i])
{
x=st[x][i];
y=st[y][i];
}
}
return st[x][0];
}
void print()
{
printf("\n");
for(int i=1;i<=n;i++)
{
for(int j=0;j<=2;j++)
{
printf("%d ",st[i][j]);
}
printf("\n");
}
printf("\n");
}
int main()
{
scanf("%d%d%d",&n,&m,&s);
for(int i=1;i<n;i++)
{
int x,y;
scanf("%d%d",&x,&y);
all[total].to=y;all[total].fr=g[x];
g[x]=total;total++;
all[total].to=x;all[total].fr=g[y];
g[y]=total;total++;
}
vis[s]=1;
dfs(s);
//print();
for(int j=1;(1<<j)<=n;j++)
{
for(int i=1;i<=n;i++)
{
st[i][j]=st[st[i][j-1]][j-1];
}
}
//print();
for(int i=1;i<=m;i++)
{
int a,b;
scanf("%d%d",&a,&b);
printf("%d\n",lca(a,b));
}
}
LCA(树剖)
#include<stdio.h>
#include<iostream>
#include<string.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const int N=5e5+5;
struct p{
int to,fr;
};
p tree[N*4];
int g[N*4];int total=1;
int n,m,s;
int size[N];
int fa[N];
int dep[N];
int wson[N];
int top[N];
void dfs1(int p)
{
size[p]=1;
int maxx=-99999;
for(int i=g[p];i;i=tree[i].fr)
{
int to=tree[i].to;
if(to!=fa[p])
{
fa[to]=p;
dep[to]=dep[p]+1;
dfs1(to);
size[p]+=size[to];
if(size[to]>maxx)
{
maxx=size[to];
wson[p]=to;
}
}
}
}
void dfs2(int p)
{
if(!top[p]) top[p]=p;
if(wson[p]!=0)
{
top[wson[p]]=top[p];
dfs2(wson[p]);
}
for(int i=g[p];i;i=tree[i].fr)
{
int to=tree[i].to;
if(to!=wson[p]&&top[to]==0)
{
dfs2(to);
}
}
}
int lca(int x,int y)
{
while(top[x]!=top[y])
{
if(dep[top[x]]<=dep[top[y]]) swap(x,y);
x=fa[top[x]];
}
if(dep[x]>dep[y])
{
swap(x,y);
}
return x;
}
int main()
{
scanf("%d%d%d",&n,&m,&s);
for(int i=1;i<=n-1;i++)
{
int x,y;
scanf("%d%d",&x,&y);
tree[total].to=y;tree[total].fr=g[x];g[x]=total;total++;
tree[total].to=x;tree[total].fr=g[y];g[y]=total;total++;
}
dfs1(s);
dfs2(s);
for(int i=1;i<=m;i++)
{
int x,y;
scanf("%d%d",&x,&y);
printf("%d\n",lca(x,y));
}
}
最短路
#include<stdio.h>
#include<iostream>
#include<string.h>
#include<vector>
#include<queue>
#include<algorithm>
using namespace std;
struct p{
int x;int value;
bool operator <(const p &a) const
{
return value>a.value;
}
};
priority_queue<p> qq;
struct q{
int from;int to;int value;
}all[200005];int dis[200005];
int g[200005];int total=1;
int vis[200005];
int n,m,s;
int main()
{
scanf("%d%d%d",&n,&m,&s);
for(int i=1;i<=m;i++)
{
int u,v,w;
scanf("%d%d%d",&u,&v,&w);
all[total].value=w;all[total].to=v;all[total].from=g[u];g[u]=total;total++;
}
memset(dis,0x3f,sizeof(dis));
dis[s]=0;//vis[s]=1;
qq.push((p){s,0});
while(qq.size())
{
p t=qq.top();qq.pop();
if(vis[t.x]==1) continue;
vis[t.x]=1;
for(int i=g[t.x];i;i=all[i].from)
{
int to=all[i].to;int v=all[i].value;
if(!vis[to]&&dis[to]>dis[t.x]+v)
{
dis[to]=dis[t.x]+v;
qq.push((p){to,dis[to]});
}
}
}
for(int i=1;i<=n;i++)
{
printf("%d ",dis[i]);
}
}
树的重心
#include<stdio.h>
#include<iostream>
#include<string.h>
#include<vector>
using namespace std;
const int N=5e4+5;
int size[N];
int dp[N];
int vis[N];
vector<int> a[N];
int n;
void dfs1(int x)
{
size[x]=1;vis[x]=1;
for(int i=0;i<a[x].size();i++)
{
int to=a[x][i];
if(vis[to]==0)
{
dfs1(to);
size[x]+=size[to];
dp[x]+=dp[to]+size[to];
}
}
}
int r;int minn=0x3f3f3f3f;
void dfs2(int x)
{
if(dp[x]<minn)
{
r=x;minn=dp[x];
}
else if(dp[x]==minn) r=min(r,x);
vis[x]=1;
for(int i=0;i<a[x].size();i++)
{
int to=a[x][i];
if(vis[to]==1) continue;
dp[to]=dp[x]-size[to]+n-size[to];
dfs2(to);
}
}
int main()
{
cin>>n;
for(int i=1;i<n;i++)
{
int x,y;
cin>>x>>y;
a[x].push_back(y);
a[y].push_back(x);
}
dfs1(1);
memset(vis,0,sizeof(vis));
dfs2(1);
printf("%d %d",r,minn);
}
持续更新ing……