poj 3169 BellmanFord—差分约束
2016-07-16 21:50
246 查看
点击打开链接
#include <iostream>
#include <algorithm>
#include <vector>
#include <cstring>
using namespace std;
const int M = 1100;
const int inf=1000000000;
int n,ml,md;
int d[M];
struct edge{
int to,cost;
edge(int a,int b):to(a),cost(b)
{
}
};
vector <edge> g[21000];
bool check()
{
bool flag=false;
for(int i=1;i<=n;i++)
{
for(int j=0;j<g[i].size();j++)
{
edge p=g[i][j];
if(d[p.to]>d[i]+p.cost)
{
d[p.to]=d[i]+p.cost;
flag=true;
return true;
}
}
}
return flag;
}
void solve(int s)
{
fill(d+1,d+n+1,inf);
d[s]=0;
// d(k)[u] 源点s到u最多经过k条边的最短路
// d(1)[u]=edge[s][u]
//d(k)[u]=min(d(k-1)[u],min(d(k-1)[j]+edge[j][u]))
for(int k=1;k<=n;k++) // 无负圈时 最短路最多只含 n-1条边
{
if(k==n) //如果k==n更新说明有负圈
{
if(check())
{
cout<<-1<<endl; // 有负圈即 i+1能排到i前面 不符合题意
return;
}
else
break;
}
for(int i=1;i<=n;i++)
{
for(int j=0;j<g[i].size();j++)
{
edge p=g[i][j];
if(d[p.to]>d[i]+p.cost)
{
d[p.to]=d[i]+p.cost;
}
}
}
}
if(d
>=inf) //n和1 无约束 可以无限大
{
cout<<-2<<endl;
return;
}
cout<<d
<<endl;
}
int main()
{
cin>>n>>ml>>md;
while(ml--)
{
int x,y,w;
cin>>x>>y>>w;
g[x].push_back(edge(y,w)); // d[x]+w>=d[y]
}
while(md--)
{
int x,y,w;
cin>>x>>y>>w;
g[y].push_back(edge(x,-w)); // d[x]+w<=d[y] -> d[y]+(-w)>=d[x]
}
for(int i=1;i<n;i++)
{
g[i+1].push_back(edge(i,0)); // d[i]<=d[i+1] -> d[i+1]+0>=d[i]
}
// u->v有多条边 取最小边 即求最短路
solve(1);
return 0;
}
#include <iostream>
#include <algorithm>
#include <vector>
#include <cstring>
using namespace std;
const int M = 1100;
const int inf=1000000000;
int n,ml,md;
int d[M];
struct edge{
int to,cost;
edge(int a,int b):to(a),cost(b)
{
}
};
vector <edge> g[21000];
bool check()
{
bool flag=false;
for(int i=1;i<=n;i++)
{
for(int j=0;j<g[i].size();j++)
{
edge p=g[i][j];
if(d[p.to]>d[i]+p.cost)
{
d[p.to]=d[i]+p.cost;
flag=true;
return true;
}
}
}
return flag;
}
void solve(int s)
{
fill(d+1,d+n+1,inf);
d[s]=0;
// d(k)[u] 源点s到u最多经过k条边的最短路
// d(1)[u]=edge[s][u]
//d(k)[u]=min(d(k-1)[u],min(d(k-1)[j]+edge[j][u]))
for(int k=1;k<=n;k++) // 无负圈时 最短路最多只含 n-1条边
{
if(k==n) //如果k==n更新说明有负圈
{
if(check())
{
cout<<-1<<endl; // 有负圈即 i+1能排到i前面 不符合题意
return;
}
else
break;
}
for(int i=1;i<=n;i++)
{
for(int j=0;j<g[i].size();j++)
{
edge p=g[i][j];
if(d[p.to]>d[i]+p.cost)
{
d[p.to]=d[i]+p.cost;
}
}
}
}
if(d
>=inf) //n和1 无约束 可以无限大
{
cout<<-2<<endl;
return;
}
cout<<d
<<endl;
}
int main()
{
cin>>n>>ml>>md;
while(ml--)
{
int x,y,w;
cin>>x>>y>>w;
g[x].push_back(edge(y,w)); // d[x]+w>=d[y]
}
while(md--)
{
int x,y,w;
cin>>x>>y>>w;
g[y].push_back(edge(x,-w)); // d[x]+w<=d[y] -> d[y]+(-w)>=d[x]
}
for(int i=1;i<n;i++)
{
g[i+1].push_back(edge(i,0)); // d[i]<=d[i+1] -> d[i+1]+0>=d[i]
}
// u->v有多条边 取最小边 即求最短路
solve(1);
return 0;
}
相关文章推荐
- 宏定义的使用
- [Usaco2003 Open]Lost Cows(线段树)
- SICP 习题 (2.31)解题总结 : 通用的Square-Tree
- Java 接口interface的详解
- Vmware虚拟机安装caffe,Ubuntu安装caffe,CPU only版本
- QT PySide 连连看小游戏
- C++之构造函数
- 【记录】DMA
- Objective-c分类
- 在eclipse中使用Tomcat8.0时出现Could not publish server ...错误
- C++实现R语言向量化运算(向量类:c 矩阵类:matrix)2015.9.11
- 企业级开源邮件系统搭建的全过程
- 竞赛排名 (Standard IO)
- 51Nod 1007 正整数分组 01背包
- 防止idea每次都默认以1.5版本编译的maven设置
- 我的Unity3D学习日记-06(自己动手制作FlappyBird)
- Maven+Spring+Spring MVC+MyBatis+MySQL,搭建SSM框架环境
- VS助手安装失败---VS番茄助手安装成功重启VS却没有效果!!!
- [c/c++]指针(4)
- verilog语法之memory存储器