洛谷 P1993 小K的农场 【差分约束+spaf】

题目链接：https://www.luogu.org/problem/P1993；

---

首先差分约束系统主要是建图方面，举个例子：

                                 ①A - B >= 3   
                                 ②B - C >= 2     
                                 ③B - D >= 5   
                                 ④C - D >= 2                   求A - D的最小值？

第一种解决方案：①+②+④ 为 7

第二种解决方案：①+③为 8；

可以得到答案，A - D 最小值为8；可以联想到，A - B >=3 相当于从B向A连一条权值为3的有向边，每个不等式都是这样，就可以得到一张图，答案就是D到A的最长路；

---

关于差分约束不等式的一些转换：

A - B <= X 相当于 B - A >= - X；

A = B 相当于A - B >= 0;

A >= B 相当于 A - B >=0;

A > B相当于 A - B >= 1;

AC代码：

/*
最短路差分约束模板题，用spaf判断存不存在负环
*/
#include<bits/stdc++.h>
#define max(a,b) a>b?a:b
#define min(a,b) a>b?b:a
using namespace std;
typedef long long ll;
const int inf=0x3f3f3f3f;
const int mod=998244353;
const int N=1e4+10;

struct node
{
    int v,ne,w;
}q[N*5];
int f[N],vis[N],dis[N];
int e,n,m,cnt[N];
///cnt数组记录进行spaf松弛操作的次数，如果大于等于n就表示存在负环

void add(int a,int b,int c)
{
    q[e].v=b;
    q[e].w=c;
    q[e].ne=f[a];
    f[a]=e++;
}

bool spaf()
{
    queue<int>que;
    que.push(1);
    for(int i=1;i<=n;i++)
    {
        dis[i]=inf;
        vis[i]=0;
    }
    for(int i=1;i<=n;i++)///图可能不是连通的，所以遍历一遍
    {
        if(vis[i]==0)
        {
            vis[i]=1;
            dis[i]=0;
            que.push(i);
            while(!que.empty())
            {
                int u=que.front();
                que.pop();
                vis[u]=0;
                for(int j=f[u];j!=-1;j=q[j].ne)
                {
                    int v=q[j].v,w=q[j].w;
                    if(dis[v]>dis[u]+w)
                    {
                        dis[v]=dis[u]+w;
                        if(vis[v]==0)
                        {
                            vis[v]=1;
                            cnt[v]++;
                            que.push(v);
                        }
                        if(cnt[v]==n) return false;
                    }
                }
            }
        }
    }
    return true;
}

int main()
{
    memset(f,-1,sizeof(f));
    e=0;
    scanf("%d %d",&n,&m);
    int a,b,c,op;
    for(int i=1;i<=m;i++)
    {
        scanf("%d",&op);
        if(op==1)
        {
            scanf("%d %d %d",&a,&b,&c);
            add(b,a,c);
        }
        else if(op==2)
        {
            scanf("%d %d %d",&a,&b,&c);
            add(a,b,-c);
        }
        else
        {
            scanf("%d %d",&a,&b);
            add(b,a,0);
        }
    }
    if(spaf()) printf("Yes\n");
    else printf("No\n");
    return 0;
}