做题日志(2019 1.6~)

Posted 2021-02-01 zh-comld

1.6

1.

https://www.cnblogs.com/ZH-comld/p/10227240.html

2、

https://www.luogu.org/problemnew/show/P1501

zzLCT，而我更zz，维护标记的方式和模板线段树2的方式一样

#include<iostream>
#include<cstdio>
#define N 100002
#define ls ch[x][0]
#define rs ch[x][1]
using namespace std;
typedef long long ll;
const int mod=51061;
int ch[N][2],fa[N],n,q; 
bool rev[N];
ll sum[N],a[N],cheng[N],jia[N],size[N];
char s[2];
inline ll rd(){
    ll x=0;char c=getchar();bool f=0;
    while(!isdigit(c)){if(c==‘-‘)f=1;c=getchar();}
    while(isdigit(c)){x=(x<<1)+(x<<3)+(c^48);c=getchar();}
    return f?-x:x;
}
inline bool ge(int x){return ch[fa[x]][1]==x;}
inline bool isroot(int x){return ch[fa[x]][0]!=x&&ch[fa[x]][1]!=x;}
inline void pushup(int x){sum[x]=(sum[ls]+sum[rs]+a[x])%mod;size[x]=size[ls]+size[rs]+1;}
inline void pushdown(int x){
    if(rev[x]){rev[x]^=1;rev[ls]^=1;rev[rs]^=1;swap(ls,rs);}
    if(cheng[x]!=1){
        (sum[ls]*=cheng[x])%=mod;(sum[rs]*=cheng[x])%=mod;
        (cheng[ls]*=cheng[x])%=mod;(cheng[rs]*=cheng[x])%=mod;
        (jia[ls]*=cheng[x])%=mod;(jia[rs]*=cheng[x])%=mod;
        (a[ls]*=cheng[x])%=mod;(a[rs]*=cheng[x])%=mod;cheng[x]=1;
    }
    if(jia[x]!=0){
        (sum[ls]+=jia[x]*size[ls])%=mod;(sum[rs]+=jia[x]*size[rs])%=mod;
        (jia[ls]+=jia[x])%=mod;(jia[rs]+=jia[x])%=mod;
        (a[ls]+=jia[x])%=mod;(a[rs]+=jia[x])%=mod;jia[x]=0;
    }
}
inline void _pushdown(int x){if(!isroot(x))_pushdown(fa[x]);pushdown(x);}
inline void rotate(int x){
    int y=fa[x],o=ge(x);
    ch[y][o]=ch[x][o^1];fa[ch[y][o]]=y;
    if(!isroot(y))ch[fa[y]][ge(y)]=x;fa[x]=fa[y];
    fa[y]=x;ch[x][o^1]=y;pushup(y);pushup(x);
}
inline void splay(int x){
    _pushdown(x);
    while(!isroot(x)){
        int y=fa[x];
        if(isroot(y))rotate(x);
        else rotate(ge(y)==ge(x)?y:x),rotate(x);
    }
}
inline void access(int x){for(int y=0;x;y=x,x=fa[x])splay(x),ch[x][1]=y,pushup(x);}
inline void makeroot(int x){access(x);splay(x);rev[x]^=1;}
inline void link(int x,int y){makeroot(x);fa[x]=y;}
inline void split(int x,int y){makeroot(x);access(y);splay(y);}
inline void cut(int x,int y){split(x,y);fa[x]=ch[y][0]=0;pushup(y);}
void dfs(int x){if(ls)dfs(ls);cout<<x<<" ";if(rs)dfs(rs);}
int main(){
    n=rd();q=rd();int x,y,xx,yy,c;
    for(int i=1;i<=n;++i)a[i]=sum[i]=1,cheng[i]=1,size[i]=1;
    for(int i=1;i<n;++i){x=rd();y=rd();link(x,y);}
    while(q--){
        scanf("%s",s);
        if(s[0]==‘-‘){x=rd();y=rd();xx=rd();yy=rd();cut(x,y);link(xx,yy);}
        else if(s[0]==‘+‘){
            x=rd();y=rd();c=rd();
            split(x,y);(sum[y]+=c*size[y])%=mod;(jia[y]+=c)%=mod;(a[y]+=c)%=mod;
        }
        else if(s[0]==‘*‘){
            x=rd();y=rd();c=rd();
            split(x,y);(sum[y]*=c)%=mod;(cheng[y]*=c)%=mod;(jia[y]*=c)%=mod;(a[y]*=c)%=mod;
        }
        else{
            x=rd();y=rd();split(x,y);printf("%lld
",sum[y]);
        }
    }
    return 0;
}

View Code

3、

https://www.cnblogs.com/ZH-comld/p/10229143.html

4、

https://www.lydsy.com/JudgeOnline/upload/201804/day2(3).pdf

什么垃圾题目，交上去发现和我年轻时写的暴力一个分。。

需要对每个时间点存一个版本，二分的时候就直接用，因为C比较小，所以复制版本还是比较快的。

注意清空一切可能会出锅的地方。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#include<vector>
#define N 209
using namespace std;
vector<int>vec[N][N];
int T,C,n,m,b[N],s[N],win[N],ans[N],t;
queue<int>q;
inline int rd(){
    int x=0;char c=getchar();bool f=0;
    while(!isdigit(c)){if(c==‘-‘)f=1;c=getchar();}
    while(isdigit(c)){x=(x<<1)+(x<<3)+(c^48);c=getchar();}
    return f?-x:x;
}
struct DINIC{
  struct edge{int n,to,l;}e[N*40];
  int tot,head[N<<1],deep[N<<1],cur[N<<1];
  DINIC(){tot=1;}
  inline void clear(){tot=1;memset(head,0,sizeof(head));}
  inline void add(int u,int v,int l){
    e[++tot].n=head[u];e[tot].to=v;head[u]=tot;e[tot].l=l;
    e[++tot].n=head[v];e[tot].to=u;head[v]=tot;e[tot].l=0;
  }  
  inline bool bfs(int s,int t){
    memset(deep,0,sizeof(deep));
    memcpy(cur,head,sizeof(head));
    q.push(s);deep[s]=1;
    while(!q.empty()){
        int u=q.front();q.pop();
        for(int i=head[u];i;i=e[i].n){
            int v=e[i].to;
            if(!deep[v]&&e[i].l){deep[v]=deep[u]+1;q.push(v);}
        }
    }
    return deep[t];
  }
  int dfs(int u,int t,int l){
    if(u==t||!l)return l;
    int f,flow=0;
    for(int &i=cur[u];i;i=e[i].n){
        int v=e[i].to;
        if(deep[v]==deep[u]+1&&(f=dfs(v,t,min(l,e[i].l)))){
            e[i].l-=f;e[i^1].l+=f;l-=f;flow+=f;
            if(!l)break; 
        }
    }
    return flow;
  }
  inline int maxflow(int s,int t){
       if(bfs(s,t))return dfs(s,t,2e9);
    return 0;
  }
}dinic[N],efs;
inline int check(int mid,int id){
    efs=dinic[mid-1];
    for(int i=1;i<=s[id];++i){
      for(int j=0;j<vec[id][i].size();++j)efs.add(id,vec[id][i][j]+n,1);  
    }
    if(efs.bfs(0,t))return 1;else return 0;
}
int main(){
    T=rd();C=rd();
    while(T--){
        n=rd();m=rd();int x;t=n+m+1;
        for(int i=1;i<=m;++i)b[i]=rd(),dinic[0].add(i+n,t,b[i]);
        for(int i=1;i<=n;++i)dinic[0].add(0,i,1); 
        for(int i=1;i<=n;++i)
          for(int j=1;j<=m;++j){
            x=rd();if(x)vec[i][x].push_back(j);
        }
        for(int i=1;i<=n;++i)s[i]=rd();
        for(int i=1;i<=n;++i){
            win[i]=m+1;ans[i]=i;
            for(int j=1;j<=m;++j){
                dinic[i]=dinic[i-1];
                for(int k=0;k<vec[i][j].size();++k)dinic[i].add(i,vec[i][j][k]+n,1);
                if(dinic[i].maxflow(0,t)){win[i]=j;break;}
            }
            int l=1,r=i;
            while(l<=r){
                int mid=(l+r)>>1;
                if(check(mid,i))ans[i]=i-mid,l=mid+1;else r=mid-1;
            }
        }
        for(int i=1;i<=n;++i)printf("%d ",win[i]);puts("");
        for(int i=1;i<=n;++i)printf("%d ",ans[i]);puts("");
        for(int i=1;i<=n;++i)for(int j=0;j<=m;++j)vec[i][j].clear();
        dinic[0].clear();
    }
    return 0;
}

View Code