首页 > 代码库 > 【差分约束系统】【强连通分量缩点】【拓扑排序】【DAG最短路】CDOJ1638 红藕香残玉簟秋,轻解罗裳,独上兰舟。
【差分约束系统】【强连通分量缩点】【拓扑排序】【DAG最短路】CDOJ1638 红藕香残玉簟秋,轻解罗裳,独上兰舟。
题意: 给定n个点(点权未知)和m条信息:u的权值>=v的权值+w 求点权的极小解和极大解(无解则输出-1)
极小解即每个点的点权可能的最小值 极大解即每个点的点权可能的最大值
题解: 差分约束系统
对于val[u]>=val[v]+w 要得到极小解,v是没有受限制的,其最小值为0 而u受到v的限制,显然,val[u]的最小值就是val[v]+w
在多条件限制下,我们用v连向u边权为w的边表示每个限制条件val[u]>=val[v]+w 那么如果得到的是拓扑图,则按拓扑序求到每个点的最长路,就得到极小解
如果得到的不是拓扑图,即图中存在回路 那么如果存在回路边权和>0,则无解(成立的前提是边权都是大于等于零的)
总的来说,求极小解的做法就是,先对val[u]>=val[v]+w建立v连向u边权为w的边 对得到的图求强连通分量,将每个强连通分量缩成一个点 若存在边权和>0的强连通分量,则无解
否则在缩点后的拓扑图上,从入度为0的点出发按拓扑序求到每个点的最长路 该最长路就是每个点的最小值
不缩点直接判是否有回路,没有回路再拓扑,这样做会出错 因为会有边权和为0的强连通分量
对于求极大解,则将条件写成val[v]<=val[u]-w的形式 建立u连向v边权为-w的边,同样求强连通分量并缩点 在缩点后的拓扑图上做最短路,该最短路就是每个点的最大值
差分约束系统如果要求最优解而非合法解,并且整个图不能转化成从某个特定点出发的话。需要进行缩点+拓扑排序……往往只在边权都大于零或者都小于零的时候才成立。
因为对于多起点的情况而言,spfa判负环是失效的,互相制约的关系也很难通过添加虚拟结点来弥补。
#include<cstdio> #include<vector> #include<cstring> #include<algorithm> #include<queue> using namespace std; int us[1000010],vs[1000010],ws[1000010]; struct Edge{ int v,w; }; vector<Edge>G[100010]; vector<int>rG[100010],tv; bool used[100010]; int cmp[100010]; int ru[100010],f[100010],ansa[100010]; int n,m; void AddEdge(int U,int V,int W){ G[U].push_back((Edge){V,W}); rG[V].push_back(U); } void dfs(int U){ used[U]=1; for(int i=0;i<G[U].size();++i){ if(!used[G[U][i].v]){ dfs(G[U][i].v); } } tv.push_back(U); } void rdfs(int U,int K){ used[U]=1; cmp[U]=K; for(int i=0;i<rG[U].size();++i){ if(!used[rG[U][i]]){ rdfs(rG[U][i],K); } } } int scc(){ memset(used,0,sizeof(used)); tv.clear(); for(int i=1;i<=n;++i){ if(!used[i]){ dfs(i); } } memset(used,0,sizeof(used)); int K=0; for(int i=tv.size()-1;i>=0;--i){ if(!used[tv[i]]){ rdfs(tv[i],++K); } } return K; } int main(){ int x,y,z; scanf("%d%d",&n,&m); for(int i=1;i<=m;++i){ scanf("%d%d%d",&us[i],&vs[i],&ws[i]); AddEdge(vs[i],us[i],ws[i]); } int sccs=scc(); for(int i=1;i<=m;++i){ if(cmp[vs[i]]==cmp[us[i]] && ws[i]>0){ puts("-1"); return 0; } } for(int i=1;i<=n;++i){ G[i].clear(); } for(int i=1;i<=m;++i){ if(cmp[vs[i]]!=cmp[us[i]]){ ++ru[cmp[us[i]]]; G[cmp[vs[i]]].push_back((Edge){cmp[us[i]],ws[i]}); } } queue<int>q; for(int i=1;i<=sccs;++i){ if(!ru[i]){ q.push(i); } } while(!q.empty()){ int U=q.front(); q.pop(); for(int i=0;i<G[U].size();++i){ f[G[U][i].v]=max(f[G[U][i].v],f[U]+G[U][i].w); --ru[G[U][i].v]; if(!ru[G[U][i].v]){ q.push(G[U][i].v); } } } if(*max_element(f+1,f+sccs+1)>100){ puts("-1"); return 0; } for(int i=1;i<=n;++i){ ansa[i]=f[cmp[i]]; } for(int i=1;i<=n;++i){ G[i].clear(); rG[i].clear(); } for(int i=1;i<=m;++i){ AddEdge(us[i],vs[i],-ws[i]); } sccs=scc(); for(int i=1;i<=m;++i){ if(cmp[us[i]]==cmp[vs[i]] && ws[i]<0){ puts("-1"); return 0; } } for(int i=1;i<=n;++i){ G[i].clear(); } for(int i=1;i<=m;++i){ if(cmp[us[i]]!=cmp[vs[i]]){ ++ru[cmp[vs[i]]]; G[cmp[us[i]]].push_back((Edge){cmp[vs[i]],-ws[i]}); } } memset(f,0x7f,sizeof(f)); for(int i=1;i<=sccs;++i){ if(!ru[i]){ q.push(i); f[i]=100; } } while(!q.empty()){ int U=q.front(); q.pop(); for(int i=0;i<G[U].size();++i){ f[G[U][i].v]=min(f[G[U][i].v],f[U]+G[U][i].w); --ru[G[U][i].v]; if(!ru[G[U][i].v]){ q.push(G[U][i].v); } } } if(*min_element(f+1,f+sccs+1)<0){ puts("-1"); return 0; } for(int i=1;i<=n;++i){ printf("%d %d\n",ansa[i],f[cmp[i]]); } return 0; }
【差分约束系统】【强连通分量缩点】【拓扑排序】【DAG最短路】CDOJ1638 红藕香残玉簟秋,轻解罗裳,独上兰舟。