洛谷 P1993 小K的农场(最短路/差分约束)

题目

小K在MC里面建立很多很多的农场，总共n个，以至于他自己都忘记了每个农场中种植作物的具体数量了，他只记得一些含糊的信息（共m个），以下列三种形式描述：

• 农场a比农场b至少多种植了c个单位的作物，
• 农场a比农场b至多多种植了c个单位的作物，
• 农场a与农场b种植的作物数一样多。

但是，由于小K的记忆有些偏差，所以他想要知道存不存在一种情况，使得农场的种植作物数量与他记忆中的所有信息吻合。

1<=n,m<=1e4

思路来源

https://www.luogu.org/problemnew/solution/P1993

题解

在差分约束系统中，

最长路无正环才有解，最短路无负环才有解，联系不等式性质显然

安利一波 国家集训队论文2009姜碧野《SPFA算法的优化及应用》

知道有个dfs优化的SPFA就好了，比常规SPFA表现好，虽然也可以被恶意卡

 

最长路，递归搜，判有无正环，无正环才有解

建一个0的虚点，和所有点相连，从而保证整张图连通

洛谷的题目，时限卡得太严了QAQ

代码

//x>=y+z y->x连z的边 跑最长路 这里是最长路 
//x<=y+z y->x连z的边 跑最短路  
//如果不能保证图连通 建一个虚点连向所有点从而保证连通 
//只建<=或>=的图 剩下的变号变过来
//x==y的 建x-y>=0和y-x>=0 
#include <iostream>
#include <algorithm> 
#include <cstring>
#include <cstdio>
using namespace std;
typedef long long ll;
const int INF=0x3f3f3f3f;
const int N=1e4+10; 
int rd()
{
    int a=0;char x=getchar();
    while(x<'0'||x>'9')x=getchar();
    while(x>='0'&&x<='9')a=(a<<3)+(a<<1)+x-48,x=getchar();
    return a;
}
int n,m,dis[N];
int head[N],cnt;
int op,a,b,c;
bool vis[N];
struct edge
{
 	int to,nex;
 	ll w;
}e[N*3];
void init()
{
	cnt=0;
	memset(head,-1,sizeof head);
}
//实际1-n的点 0是虚点保证图连通 
bool spfa(int u)
{
	vis[u]=1;
	for(int i=head[u];~i;i=e[i].nex)
	{
		int v=e[i].to;
		if(dis[v]<dis[u]+e[i].w)
		{
			dis[v]=dis[u]+e[i].w;
			if(vis[v])return 0;
			if(!spfa(v))return 0;
		}
	}
	vis[u]=0;
	return 1;
}
void add(int u,int v,ll w)
{
	e[cnt].to=v;
	e[cnt].w=w;
	e[cnt].nex=head[u];
	head[u]=cnt++;
}
int main()
{ 
   init();
   scanf("%d%d",&n,&m);
   while(m--)
   {
		op=rd(),a=rd(),b=rd();
		if(op==1)c=rd(),add(b,a,c);
		else if(op==2)c=rd(),add(a,b,-c);
		else if(op==3)add(a,b,0),add(b,a,0); 
   }
   for(int i=1;i<=n;++i)
   {
   	   add(0,i,0);
	   dis[i]=-INF;
   }
   if(!spfa(0))puts("No");
   else puts("Yes");
   return 0;
}