luogu1073 最优贸易 (tarjan+dp)

tarjan缩点，然后按照拓扑序，做1号点能到达的点的答案
具体做法是对每个点记一个min[i],max[i],vis[i]和ans[i]
做拓扑序的时候，假设在从u点开始做，有边u到v，如果vis[u]=1，则则
vis[v]=1（初始时vis[bel[1]]=1）；
更新在v点及以前买进的最小进价：min[v]=min{min[v],min[u]}
统计在v点卖出的答案：ans[v]=max{ans[u],ans[v],max[v]-min[v]}
则最后答案就是ans[bel[N]]

貌似还有做两遍spfa的做法...不管了

 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 #include<queue>
 #include<vector>
 #define LL long long int
 #define inf 0x3f3f3f3f
 using namespace std;
 const int maxn=,maxm=;
 
 LL rd(){
     LL x=;char c=getchar();int neg=;
     while(c<''||c>''){if(c=='-') neg=-;c=getchar();}
     while(c>=''&&c<='') x=x*+c-'',c=getchar();
     return x*neg;
 }
 
 int N,M;
 int eg[maxm][],egh[maxn],ect;
 int eg2[maxm][],egh2[maxn],ect2;
 int val[maxn],in[maxn],bel[maxn];
 int dfn[maxn],low[maxn],stk[maxn],sct,tot,pct;
 int mi[maxn],ma[maxn],ans[maxn];
 vector<int> hav[maxn];
 bool instk[maxn],flag[maxn],vis[maxn];
 
 inline void adeg(int a,int b){
     eg[ect][]=b;eg[ect][]=egh[a];egh[a]=ect++;
 }
 inline void adeg2(int a,int b){
     eg2[ect2][]=b;eg2[ect2][]=egh2[a];egh2[a]=ect2++;
     in[b]++;
 }
 
 void tarjan(int x){
     dfn[x]=low[x]=++tot;instk[x]=;stk[++sct]=x;
     for(int i=egh[x];i!=-;i=eg[i][]){
         int j=eg[i][];
         if(instk[j]) low[x]=min(low[x],dfn[j]);
         else if(!dfn[j]) {
             tarjan(j);low[x]=min(low[x],low[j]);
         }
     }if(dfn[x]==low[x]){
         pct++;mi[pct]=inf;ma[pct]=-inf;
         while(sct){
             mi[pct]=min(mi[pct],val[stk[sct]]);
             ma[pct]=max(ma[pct],val[stk[sct]]);
             bel[stk[sct]]=pct;
             instk[stk[sct]]=;
             hav[pct].push_back(stk[sct]);
             if(stk[sct--]==x) break;
         }
     }
 }
 
 int main(){
     int i,j,k;
     //freopen("testdata.in","r",stdin);
     N=rd(),M=rd();
     for(i=;i<=N;i++) val[i]=rd();
     memset(egh,-,sizeof(egh));memset(egh2,-,sizeof(egh2));
     for(i=;i<=M;i++){
         int a=rd(),b=rd(),c=rd();
         adeg(a,b);if(c==) adeg(b,a);
     }for(i=;i<=N;i++) if(!dfn[i]) tarjan(i);
 
     for(i=;i<=pct;i++){
         memset(flag,,pct+);
         for(j=;j<hav[i].size();j++){
             for(k=egh[hav[i][j]];k!=-;k=eg[k][]){
                 if(i==bel[eg[k][]]) continue;
                 if(!flag[bel[eg[k][]]]) adeg2(i,bel[eg[k][]]);
                 flag[bel[eg[k][]]]=;
             }
         }
     }
     queue<int> q;vis[bel[]]=;
     for(i=;i<=pct;i++){
         if(!in[i]) q.push(i);
     }
     while(!q.empty()){
         int p=q.front();q.pop();
         for(i=egh2[p];i!=-;i=eg2[i][]){
             j=eg2[i][];
             if(vis[p]){
                 vis[j]=;ans[j]=max(ans[j],ans[p]);
                 mi[j]=min(mi[j],mi[p]);
                 ans[j]=max(ans[j],ma[j]-mi[j]);
             }
             if(!(--in[j])) q.push(j);
         }
     }
     printf("%d\n",ans[bel[N]]);
 }