• [LUOGU] P3275 [SCOI2011]糖果


    显然的差分约束,除了题目给出的约束以外,边界是(A_x-Ageq1)
    判环有两种方法,要不就是SPFA记录最长路经过次数,要不就是tarjan缩了SCC,看SCC内有没有正权边

    #include<algorithm>
    #include<iostream>
    #include<cstdio>
    #include<queue>
    #include<vector>
    #define int long long
    using namespace std;
    
    inline int rd(){
      int ret=0,f=1;char c;
      while(c=getchar(),!isdigit(c))f=c=='-'?-1:1;
      while(isdigit(c))ret=ret*10+c-'0',c=getchar();
      return ret*f;
    }
    #define space() putchar(' ')
    #define nextline() putchar('
    ')
    void pot(int x){if(!x)return;pot(x/10);putchar('0'+x%10);}
    void out(int x){if(!x)putchar('0');if(x<0)putchar('-'),x=-x;pot(x);}
    
    const int MAXN=4000005;
    
    int n,m;
    
    int nex[MAXN],to[MAXN],wi[MAXN];
    int ecnt,head[MAXN];
    inline void add(int x,int y,int w){
    //	cout<<x<<" "<<y<<" "<<w<<endl;
        nex[++ecnt] = head[x];
        to[ecnt] = y;
        wi[ecnt]=w;
        head[x] = ecnt;
    }
    
    int dis[MAXN],inq[MAXN],cnt[MAXN];
    queue<int> Q;
    bool spfa(int st){
    //	memset(dis,0x3f,sizeof(dis));
        dis[st]=0;Q.push(st);
        while(!Q.empty()){
            int top=Q.front();Q.pop();inq[top]=0;
     		
     		cnt[top]++;if(cnt[top]>n) return 1;
            for(int i=head[top];i;i=nex[i]){
                int v=to[i];
                if(dis[top]+wi[i]>dis[v]){
                    dis[v]=dis[top]+wi[i];
                    if(!inq[v]) Q.push(v),inq[v]=1;
                }
            }
        }
        return 0;
    }
    
    int dfn[MAXN],low[MAXN],tim;
    int sta[MAXN],ins[MAXN],top; 
    int bl[MAXN],scc;
    vector<int> vec[MAXN];
    void tarjan(int x){
        dfn[x]=low[x]=++tim;
        ins[sta[++top]=x]=1;
        for(int i=head[x];i;i=nex[i]){
            int v=to[i];
            if(!dfn[v]){
                tarjan(v);
                low[x]=min(low[x],low[v]);
            }else if(ins[v]){
                low[x]=min(low[x],dfn[v]);	
            }
        }
        if(dfn[x]!=low[x]) return;
        int elm=-1;++scc;
        do{
            elm=sta[top--];
            ins[elm]=0;
            bl[elm]=scc;
            vec[scc].push_back(elm);
        }while(elm!=x);
    }
    
    
    signed main(){
        n=rd();m=rd();
        int t,x,y;
        for(int i=1;i<=m;i++){
            t=rd();x=rd();y=rd();
            if(t==1) add(x,y,0),add(y,x,0);
            if(t==2) add(x,y,1);
            if(t==3) add(y,x,0);
            if(t==4) add(y,x,1);
            if(t==5) add(x,y,0);
        }
        for(int i=n;i>=1;i--) add(n+1,i,1);
        for(int i=1;i<=n+1;i++){
            if(!dfn[i]) tarjan(i);	
        }
        for(int i=1;i<=scc;i++){
            int sz=vec[i].size();
            for(int j=0;j<sz;j++){
                for(int k=head[vec[i][j]];k;k=nex[k]){
                    int v=to[k];
                    if(bl[v]!=i) continue;
                    if(wi[k]!=0) return puts("-1"),0;	
                }
            }
        }
        if(spfa(n+1)) return puts("-1"),0;
        int ans=0;
        for(int i=1;i<=n;i++) ans+=dis[i];
        out(ans);
    }
    
  • 相关阅读:
    数据结构之数组
    数据结构之链表
    MongoDB使用笔记
    数据结构之ArrayList
    java设计模式之--装饰者模式
    JAVA设计模式之--模板方法模式
    HashTable、HashMap、ConcurrentHashMap源码分析
    Docker使用笔记
    First-blog:解决mybatis 用mysql进行模糊搜索时,查不了中文问题
    css cursor: url() 使用火狐浏览器问题,鼠标没有效果
  • 原文地址:https://www.cnblogs.com/ghostcai/p/9836695.html
Copyright © 2020-2023  润新知