洛谷1037 NOIP2009 最優貿易

好吧 我太zz了 居然花了3個小時 還是找了一組大數據測的 wa無數次

有一個很簡單的思路 就是正反跑2遍spfa存最大值和最小值 每個點做差取max 這個想法好像比較(特別)好想 也很好實現

可是我這個蒟蒻先寫了個tarjan 然後就一直順着這個思路 陷入深坑 不過跑完灰常快(O(∩_∩)O哈哈~) 我的思路是先縮點 然後建反向邊跑一邊 看那些點是可以到的 然後正向建邊dfs 當然得剪枝啦 我這裏t了好幾次

具體的看代碼吧 考前切切noip 順帶把變量名設成各位dalao膜一膜

#include<bits/stdc++.h>
using namespace std;
const int N=1e5+5;
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int n,m,ans,sum,cnt,top,tot,dfn[N],low[N],head[N],bel[N],mn[N],mx[N],in[N],v[N],a[5*N],b[5*N],c[5*N],amaz[N];
bool vis[N],to[N];
struct Edge{
    int v,nxt;
}e[N*2]; 
void add(int u,int v){
    e[++tot].v=v,e[tot].nxt=head[u],head[u]=tot; 
}
void fc(int x){
    to[x]=1;vis[x]=1;
    for(int i=head[x];i;i=e[i].nxt){
        int j=e[i].v;
        if(!vis[j]) fc(j);  
    }
}
void dfs(int x,int need,int zjq){
//  bug(need);
    if(x==bel[1]) ans=max(ans,zjq);//,cout<<x<<" "<<need<<" "<<zjq<<endl;
    else
        for(int i=head[x];i;i=e[i].nxt){
            int j=e[i].v;
            int wly=max(need,mx[j]);
            if(amaz[j]<wly){
                amaz[j]=wly;        
                dfs(j,wly,max(zjq,wly-mn[j])); 
            }
            else if(to[j]){
                //cout<<x<<" "<<need<<" "<<zjq<<endl;
                ans=max(ans,zjq);
            }
        }
}
void tarjan(int x){
    dfn[x]=low[x]=++cnt;
    vis[x]=1;in[++top]=x;
    for(int i=head[x];i;i=e[i].nxt){
        int j=e[i].v;
        if(!dfn[j]){
            tarjan(j); 
            low[x]=min(low[j],low[x]);
        }
        else if(vis[j]){
            low[x]=min(low[x],dfn[j]);
        }
    }
    if(dfn[x]==low[x]){
//      cout<<x<<" : ";
        mn[++sum]=INT_MAX;
        do{
            vis[in[top]]=0;
            bel[in[top]]=sum;
//          cout<<in[top]<<" ";
            mx[sum]=max(mx[sum],v[in[top]]);
            mn[sum]=min(mn[sum],v[in[top]]);
        }while(in[top--]!=x);
//      cout<<endl;
    }
}
int main(){
#ifdef Devil_Gary
    freopen("in.txt","r",stdin);
    freopen("out.txt","w",stdout);
#endif
    n=read(),m=read();
    for(int i=1;i<=n;i++) v[i]=read();
    for(int i=1;i<=m;i++){
        a[i]=read(),b[i]=read(),c[i]=read();
        if(c[i]==1) add(a[i],b[i]);
        else add(a[i],b[i]),add(b[i],a[i]);
    }
    tarjan(1);
    for(int i=1;i<=tot;i++) e[i].v=e[i].nxt=0;tot=1;
    for(int i=1;i<=n;i++) head[i]=vis[i]=0;tot=1;
    for(int i=1;i<=m;i++) if(c[i]==1) if(bel[a[i]]!=bel[b[i]]) add(bel[a[i]],bel[b[i]]);
    fc(bel[1]);
    for(int i=1;i<=tot;i++) e[i].v=e[i].nxt=0;tot=1;
    for(int i=1;i<=n;i++) head[i]=0;tot=1;
    for(int i=1;i<=m;i++) if(c[i]==1) if(bel[a[i]]!=bel[b[i]]) add(bel[b[i]],bel[a[i]]);
    dfs(bel[n],mx[bel[n]],ans);
    printf("%d",ans);
}
發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章