【差分约束系统】【强连通分量缩点】【拓扑排序】【DAG最短路】CDOJ1638 红藕香残玉簟秋,轻解罗裳,独上兰舟。

时间:2023-03-10 02:32:30
【差分约束系统】【强连通分量缩点】【拓扑排序】【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;
}