图论(7)负环和差分约束

news2024/11/23 21:36:11

一、概念

给定一张有向图,如果存在一个环,环上各边权值之和是负数,则称这个环为负环。

判断方式:bellman-ford算法和spfa算法。+抽屉原理

这里只介绍spfa。设立cnt数组表示从1到x的最短路径包含的边数,如果cnt[i]大于等于n,则存在负环。

二、例题

904. 虫洞

 虫洞是负边,存在负环则true。

#include<iostream>
#include<algorithm>
#include<cstring>
#include<queue>
using namespace std;

const int N =510,M=5210;

int f,n,m1,m2;
int dist[N];
int h[N],e[M],ne[M],w[M],idx;
bool st[N];
int cnt[N];
int q[N];

void add(int a,int b,int c)
{
    e[idx]=b,w[idx]=c,ne[idx]=h[a],h[a]=idx++;
}

bool spfa()
{
    memset(dist,0,sizeof dist);
    memset(cnt, 0, sizeof cnt);
    queue<int> q;
    for(int i=1;i<=n;i++)
    {
        q.push(i);
        st[i]=true;
    }

    while(q.size())
    {
        int t=q.front();
        q.pop();
        st[t]=false;

        for(int i=h[t];~i;i=ne[i])
        {
            int j=e[i];
            if(dist[j]>dist[t]+w[i])
            {
                dist[j]=dist[t]+w[i];
                cnt[j]=cnt[t]+1;
                if(cnt[j]>=n) return true;
                if(!st[j])
                {
                    st[j]=true;
                    q.push(j);
                }
            }
        }
    }
    return false;
}


int main()
{
    cin>>f;
    while(f--)
    {
        cin>>n>>m1>>m2;
        memset(h,-1,sizeof h);
        idx=0;

        for(int i=0;i<m1;i++)
        {
            int a,b,c;
            cin>>a>>b>>c;
            add(a,b,c),add(b,a,c);
        }
        for(int i=0;i<m2;i++)
        {
            int a,b,c;
            cin>>a>>b>>c;
            add(a,b,-c);
        }
        if(spfa()) puts("YES");
        else puts("NO");
    }
    return 0;
}

作者:yankai
链接:https://www.acwing.com/activity/content/code/content/4455339/
来源:AcWing
著作权归作者所有。商业转载请联系作者获得授权,非商业转载请注明出处。

 361. 观光奶牛

 这个和01分数规划有关,涉及01分数规划,一般会用二分。

设二分答案为mid:

1.如果图中存在一个环,若\sum w(e)*mid-\sum w(v)<0,则说明

\exists S=\left \{ v,e \right \} \ \ st.\frac{\sum w(v)}{\sum w(e)}>mid

即本题所求的最大值大于mid

2如果对任意的环,都大于等于0,则说明最大值小于等于mid。

因此可以二分,并且我们发现(1)中的环,实际上可以替代为:没有点权,边权e(x,y)的权值是mid*w(e)-w(x),即本来的边权乘以mid再减去入点的权值。的这样一个环。因此可以构建一个新图。

(1)中的条件变为,这样的新图有负环。

因此二分求解答案即可。

#include<iostream>
#include<algorithm>
#include<cstring>
#include<queue>
using namespace std;

const int N =1010,M=5010;

int n,m;
int wf[N];
int h[N],e[M],ne[M],wt[M],idx;
bool st[N];
int cnt[N];
double dist[N];
int q[N];

void add(int a,int b,int c)
{
    e[idx]=b,wt[idx]=c,ne[idx]=h[a],h[a]=idx++;
}

bool check(double  mid)
{
    memset(dist,0,sizeof dist);
    memset(cnt,0,sizeof cnt);
    memset(st,0,sizeof st);
    queue<int> q;

    for(int i=1;i<=n;i++)
    {
        st[i]=true;
        q.push(i);
    }
    while(q.size())
    {
        int t=q.front();
        q.pop();
        st[t]=false;

        for(int i=h[t];i!=-1;i=ne[i])
        {
            int j=e[i];
            if(dist[j]>dist[t]+mid*wt[i]-wf[t])
            {
                dist[j]=dist[t]+mid*wt[i]-wf[t];
                cnt[j]=cnt[t]+1;
                if(cnt[j]>=n) return true;
                if(!st[j])
                {
                    st[j]=true;
                    q.push(j);
                }
            }
        }
    }
    return false;

}


int main()
{
    cin>>n>>m;
    for(int i=1;i<=n;i++) cin>>wf[i];
    memset(h,-1,sizeof h);

    for(int i=0;i<m;i++)
    {
        int a,b,c;
        cin>>a>>b>>c;
        add(a,b,c);
    }
    double l=0,r=1e6;
    while(r-l>1e-4)
    {
        double mid=(l+r)/2;
        if(check(mid)) l=mid;
        else r=mid;
    }
    printf("%.2lf",l);
    return 0;


}

作者:yankai
链接:https://www.acwing.com/activity/content/code/content/4455633/
来源:AcWing
著作权归作者所有。商业转载请联系作者获得授权,非商业转载请注明出处。

1165. 单词环

 

 分析:

  1. 建图:最直观的建图方式,是把每个单词看做节点,然后连接可接龙的单词。但是这样最多有1e5个点,有1e10条边,存不下。因此要换方法。考虑一个对偶的建图方式,将每一个单词看作一条边,其开头两个字符和结尾两个字符为它两边的点。这样建图的话,节点数就缩小到了676个(26∗26)边数为1e5条。
  2. 01分数规划:我们要求\frac{\sum len}{s}=k的最大值,设二分答案mid。当k>mid时:\sum len-s*mid>0,因此我们把边权设置为len[i]−mid。
  3. 优化:spfa求负环时,可能一直跑不出结果。可以采取一种比较取巧的方法:当求最长路时,经过的点大于某一个数时,我们就可以武断地认为当前图中存在一个正环.
#include<iostream>
#include<algorithm>
#include<cstring>
#include<queue>
using namespace std;

const int N =700,M=100010;

int n;
double dist[N];
int cnt[N];
int h[N],e[M],ne[M],w[M],idx;
bool st[N];

void add(int a,int b,int c)
{
    e[idx]=b,w[idx]=c,ne[idx]=h[a],h[a]=idx++;
}

bool check(double mid)
{
    queue<int> q;
    memset(dist,0,sizeof dist);
    memset(cnt,0,sizeof cnt);
    memset(st,0,sizeof st);

    for(int i=0;i<676;i++)
    {
        st[i]=true;
        q.push(i);
    }

    int count=0;
    while(q.size())
    {
        int t=q.front();
        q.pop();
        st[t]=false;
        for(int i=h[t];i!=-1;i=ne[i])
        {
            int j=e[i];
            if(dist[j]>dist[t]+mid-w[i])
            {
                dist[j]=dist[t]+mid-w[i];
                if(++count>10000) return true;
                cnt[j]=cnt[t]+1;
                if(cnt[j]>=N) return true;
                if(!st[j])
                {
                    st[j]=true;
                    q.push(j);
                }
            }
        }
    }
    return false;
}

int main()
{
    char str[1010];
    while(cin>>n,n)
    {
        memset(h,-1,sizeof h);
        idx=0;

        for(int i=0;i<n;i++)
        {
            cin>>str;
            int len=strlen(str);
            if(len>=2)
            {
                int left=(str[0] - 'a') * 26 + str[1] - 'a';
                int right = (str[len - 2] - 'a') * 26 + str[len - 1] - 'a';
                add(left,right,len);
            }
        }

        if(!check(0)) puts("No solution");
        else
        {
            double l=0,r=1010;
            while(r-l>1e-5)
            {
                double mid=(l+r)/2;
                if(check(mid)) l=mid;
                else r=mid;
            }
            printf("%.2lf\n",r);
        }
    }
    return 0;
}

作者:yankai
链接:https://www.acwing.com/activity/content/code/content/4455934/
来源:AcWing
著作权归作者所有。商业转载请联系作者获得授权,非商业转载请注明出处。

三、差分约束

差分约束就是一组N元一次不等式。给出了多个约束条件形如:X_i-X_j\leqslant c_k。我们要解决的问题就是求出一组X,使得所有约束条件得到满足。

我们可以把X_i-X_j\leqslant c_k看做 j 节点向 i 节点连一条长度为ck的有向边。

则如果a1 a2....an是一组解,则对任意的常数delta,a+delta也是一组解。所以不妨先求一组负数解。假设所有ai小于等于0。我们再增加一个超级原点a0=0,则多了N个形如a_i-a_0\leqslant 0的约束条件。即我们把0号点向所有点连一条长度为0的边。

在我们建的图上跑最短路,如果存在负环,则说明有ai<ai,这是不可能满足的,所以差分系统无解,否则,ai=dist[i]就是一组解。

如果约束条件不等号反向,我们仍然可以不改变图,只不过改为计算单源最长路,如果有正环则无解。

 参考:AcWing 1169. 糖果 - AcWing

AcWing 1169. 糖果 - AcWing 

也就是说,如果是求最小值,则求最长路。求最大值,则求最短路。

四、例题

1169. 糖果

 

 注意超级原点。因此本题求最小值,跑最长路判断即可,判断有无正环

#include<iostream>
#include<algorithm>
#include<cstring>
#include<stack>
using namespace std;


const int N =1e5+10,M=3e5+10;

int n,m;
int h[N],e[M],ne[M],w[M],idx;
long long  dist[N];
bool st[N];
int cnt[N];


void add(int a,int b,int c)
{
    e[idx]=b,w[idx]=c,ne[idx]=h[a],h[a]=idx++;
}

int spfa()
{
    memset(dist,-0x3f,sizeof dist);
    dist[0]=0;
    stack<int> q;
    q.push(0);
    st[0]=true;

    while(q.size())
    {
        int t=q.top();
        q.pop();
        st[t]=false;

        for(int i=h[t];~i;i=ne[i])
        {
            int j=e[i];
            if(dist[j]<dist[t]+w[i])
            {
                dist[j]=dist[t]+w[i];
                cnt[j]=cnt[t]+1;

                if(cnt[j]>=n+1) return false;
                if(!st[j]) 
                {
                    q.push(j);
                    st[j]=true;
                }
            }
        }
    }
    return true;
}

int main()
{
    cin>>n>>m;
    memset(h,-1,sizeof h);
    for(int i=0;i<m;i++)
    {
        int a,b,x;
        cin>>x>>a>>b;
        if(x==1) add(a,b,0),add(b,a,0);
        else if(x==2) add(a,b,1);
        else if(x==3) add(b,a,0);
        else if(x==4) add(b,a,1);
        else add(a,b,0);
    }

    for(int i=1;i<=n;i++) add(0,i,1);

    if(!spfa()) cout<<"-1";
    else 
    {
        long long  res=0;
        for(int i=1;i<=n;i++) res+=dist[i];
        cout<<res;
    }


}

作者:yankai
链接:https://www.acwing.com/activity/content/code/content/4482434/
来源:AcWing
著作权归作者所有。商业转载请联系作者获得授权,非商业转载请注明出处。

362. 区间(前缀和)

 AcWing 362. 区间 - AcWing

 

#include<iostream>
#include<algorithm>
#include<cstring>
#include<queue>
using namespace std;

const int N =5e4+10,M=2e5;

int n;
int h[N],e[M],ne[M],w[M],idx;
int dist[N];
bool st[N];

void add(int a,int b,int c)
{
    e[idx]=b,w[idx]=c,ne[idx]=h[a],h[a]=idx++;
}

void spfa()
{
    memset(dist,-0x3f,sizeof dist);
    queue<int> q;
    dist[0]=0;
    q.push(0);
    st[0]=true;
    
    while(q.size())
    {
        int t=q.front();
        q.pop();
        st[t]=false;
        
        for(int i=h[t];~i;i=ne[i])
        {
            int j=e[i];
            if(dist[j]<dist[t]+w[i])
            {
                dist[j]=dist[t]+w[i];
                if(!st[j])
                {
                    q.push(j);
                    st[j]=true;
                }
            }
        }
    }
}

int main()
{
    cin>>n;
    memset(h,-1,sizeof h);
    
    for(int i=1;i<N;i++)
    {
        add(i,i-1,-1);
        add(i-1,i,0);
    }
    
    for(int i=0;i<n;i++)
    {
        int a,b,c;
        cin>>a>>b>>c;
        a++,b++;
        add(a-1,b,c);
    }
    
    spfa();
    cout<<dist[50001];
    
}

1170. 排队布局

 无穷即判断1号点和n号点有无约束关系即,即连不连通。

 

#include<iostream>
#include<algorithm>
#include<cstring>
#include<queue>

using namespace std;

const int N =1010,M=10000+10000+1010,INF=0x3f3f3f3f;

int n,m1,m2;
int dist[N];
bool st[N];
int h[N],e[M],ne[M],w[M],idx;
int cnt[N];

void add(int a,int b,int c)
{
    e[idx]=b,w[idx]=c,ne[idx]=h[a],h[a]=idx++;
}

bool spfa(int size)
{
    memset(dist,0x3f,sizeof dist);
    memset(st,0,sizeof st);
    memset(cnt,0,sizeof cnt);
    queue<int> q;

    for(int i=1;i<=size;i++)
    {
        q.push(i);
        st[i]=true;
        dist[i]=0;
    }
    while(q.size())
    {
        int t=q.front();
        q.pop();
        st[t]=false;

        for(int i=h[t];i!=-1;i=ne[i])
        {
            int j=e[i];
            if(dist[j]>dist[t]+w[i])
            {
                dist[j]=dist[t]+w[i];
                cnt[j]=cnt[t]+1;
                if(cnt[j]>=n) return true;
                if(!st[j])
                {
                    q.push(j);
                    st[j]=true;
                }
            }
        }
    }
    return false;
}

int main()
{
    cin>>n>>m1>>m2;
    memset(h,-1,sizeof h);
    for(int i=1;i<n;i++) add(i+1,i,0);

    for(int i=0;i<m1;i++)
    {
        int a,b,c;
        cin>>a>>b>>c;
        if(a>b) swap(a,b);
        add(a,b,c);
    }
    for(int i=0;i<m2;i++)
    {
        int a,b,c;
        cin>>a>>b>>c;
        if(a>b) swap(a,b);
        add(b,a,-c);
    }

    if(spfa(n)) puts("-1");//出现矛盾 有负环
    else
    {
        spfa(1);
        if(dist[n]==INF) puts("-2");
        else cout<<dist[n];
    }
    return 0;
}

作者:yankai
链接:https://www.acwing.com/activity/content/code/content/4494625/
来源:AcWing
著作权归作者所有。商业转载请联系作者获得授权,非商业转载请注明出处。

393. 雇佣收银员

 

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/185447.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

JZ65 不用加减乘除做加法

【答案解析】&#xff1a;十进制相加思想&#xff1a; 1507 &#xff0c; 先计算不考虑进位的相加结果 12 &#xff08;因为 57 的不考虑进位的结果是 2 &#xff0c;遇 10 进位嘛&#xff09;&#xff0c;然后计算进位 57 进位是 10 &#xff0c;则 10 与 12 再次相加&#xf…

2023年山东最新交安安全员考试题库及答案

百分百题库提供交安安全员考试试题、交安安全员考试真题、交安安全员证考试题库等&#xff0c;提供在线做题刷题&#xff0c;在线模拟考试&#xff0c;助你考试轻松过关。 35.建设行政主管部门或者其他有关部门可以将施工现场的监督检查委托给建设工程&#xff08;&#xff09;…

Visual Studio 17.5: 有关 C++ 的新特性预览

Visual Studio 17.5 在 C 代码编辑方面带来了一些新的改进&#xff0c;这些改进包括&#xff1a;C 代码中的大括号对着色&#xff0c;拼写检查&#xff0c;多合一搜索&#xff0c;重新设计的成员列表以及宏展开改进等。上述这些改进都融入到了 Visual Studio 的最新预览版本 17…

CTF-Horizontall HackTheBox渗透测试(一)

** 0X01 简介** #Horizontall#难度是一个相对“简单”的 CTF Linux盒子。该CTF环境涵盖了通过利用Strapi RCE 漏洞并使用内部应用程序 (Laravel) 将隧道传输到本地计算机&#xff0c;并且在 Laravel v 7.4.18 上运行 漏洞PoC最来提升权限拿到root用户权限。 ** 1.1信息收集**…

GuLi商城-人人开源搭建后台管理系统

参考&#xff1a; 谷粒商城-基础篇(详细流程梳理代码) 谷粒商城-day01 项目的基本搭建_周周写不完的代码的博客-CSDN博客_谷粒商城 若依 谷粒商城分布式基础篇1-个人版_断河愁的博客-CSDN博客_谷粒商城 https://blog.csdn.net/yudbl/category_11902060.html 学习路线 源代…

博客搭建教程2-hexo框架初搭建

文章目录1 前言2 准备工作3 hexo安装4 共享文件夹创建(额外)1 前言 本次教程选用hexo来搭建博客&#xff0c;hexo是一个开源的架构&#xff0c;只需要进行简单的操作就可以拥有自己的博客。 参考网站&#xff1a; hexo官网 注意:下面的命令在root下进行&#xff0c;在日常的工…

【雷丰阳-谷粒商城 】【分布式基础篇-全栈开发篇】【04】跨域_OSS_后端校验

持续学习&持续更新中… 学习态度&#xff1a;守破离 【雷丰阳-谷粒商城 】【分布式基础篇-全栈开发篇】【04】跨域问题解决实现逻辑删除文件存储普通上传云存储阿里云OSS简介术语简单使用使用SpringCloudAlibaba—oss服务端签名后直传普通上传方式&#xff1a;服务端签名后…

Win10安装ElasticSearch笔记

1、安装前准备条件因为ElasticSearch7.17需要JDK1.8的支持&#xff0c;首先确保你的win10已经提前安装好了jdk8的版本ElasticSearch支持的JDK最低版本是1.8.0。ElasticSearch7.17及以下的版本最低版本是JDK1.8.0ElasticSearch8.0及以上的版本最低版本是JDK162、官网下载ES安装包…

Java开发环境配置 “IntelliJ IDEA”(超详细整理,适合新手入门)

前言 &#x1f4dc; “ 作者 久绊A ” 专注记录自己所整理的Java、web、sql等&#xff0c;IT技术干货、学习经验、面试资料、刷题记录&#xff0c;以及遇到的问题和解决方案&#xff0c;记录自己成长的点滴 目录 前言 一、IDEA的介绍 1、大概介绍 2、详细介绍 二、Intelli…

Python中的logging模块

软件开发中通过日志记录程序的运行情况是一个开发的好习惯&#xff0c;对于错误排查和系统运维都有很大帮助。Python标准库自带日志模块&#xff0c;程序的日志功能直接调用标准库的日志模块即可通过日志&#xff0c;开发者可以清楚的了解发生了哪些事件&#xff0c;包括出现了…

Linux随记(五)

一、已用statefulset创了两个nginx副本 web-0 和 web-1 &#xff0c;目的将各自容器的hostname重定向到index.html。但显示的是master节点的hostname。 怎么写才是搞成 pod里面的主机名&#xff1f;#最终解决方法&#xff1a; 在sh -c 双引号里$符号前面加上\转义。 或者 sh -…

Redis原理篇(二)网络模型

一、用户空间和内核空间 应用需要通过Linux内核与硬件交互。 内核本质也是应用&#xff0c;运行的时候也需要CPU资源、内存资源。用户应用也在消耗这些资源。 为了避免用户应用导致冲突甚至内核崩溃&#xff0c;用户应用与内核是分离的&#xff1a; 进程的寻址空间会划分为两…

《安富莱嵌入式周报》第301期:ThreadX老大离开微软推出PX5 RTOS第5代系统,支持回流焊的自焊接PCB板设计,单色屏实现多级灰度播放视频效果

往期周报汇总地址&#xff1a;嵌入式周报 - uCOS & uCGUI & emWin & embOS & TouchGFX & ThreadX - 硬汉嵌入式论坛 - Powered by Discuz! 祝大家开工大吉 视频版&#xff1a; https://www.bilibili.com/video/BV1GT411o7zr 1、ThreadX老大离开微软&…

又一个开源工具搞完了,工作效率直接翻倍

&#x1f3e1; 博客首页&#xff1a;派 大 星 ⛳️ 欢迎关注 &#x1f433; 点赞 &#x1f392; 收藏 ✏️ 留言 &#x1f3a2; 本文由派大星原创编撰 &#x1f6a7; 系列专栏&#xff1a;《开源专栏》 &#x1f388; 本系列主要输出作者自创的开源项目 &#x1f517; 作品&…

JavaScript switch 语句

JavaScript switch 语句 switch 语句评估一个表达式&#xff0c;将表达式的值与case子句匹配&#xff0c;并执行与该情况相关联的语句。 举例说明 const expr Papayas; switch (expr) {case Oranges:console.log(Oranges are $0.59 a pound.);break;case Mangoes:case Papa…

力扣刷题|104.二叉树的最大深度、559.n 叉树的最大深度、111.二叉树的最小深度、222.完全二叉树的节点个数

文章目录LeetCode 104.二叉树的最大深度题目链接&#x1f517;思路LeetCode 111.二叉树的最小深度题目链接&#x1f517;思路LeetCode 222.完全二叉树的节点个数题目链接&#x1f517;普通二叉树求法针对完全二叉树解法LeetCode 104.二叉树的最大深度 题目链接&#x1f517; …

WebAssembly 助力 桌面端运行web

WebAssembly (abbreviated Wasm) is a binary instruction format for a stack-based virtual machine. Wasm is designed as a portable compilation target for programming languages, enabling deployment on the web for client and server applications.WebAssembly 是基…

缺陷分析测试人员绕不开的知识点

缺陷分析也是测试工程师需要掌握的一个能力&#xff0c;但是很多时候大家只记得要提交缺陷、统计缺陷情况&#xff0c;而忽视了缺陷分析。 其实每个项目的缺陷记录都是有很大价值的。在测试阶段分析当前缺陷情况&#xff0c;及时发现存在的问题并调整测试策略&#xff0c;才能…

区间一维dp史上最细总结(绝对干货,还不会的一定要进来)

那年初夏&#xff08;三&#xff09; 注&#xff1a;此部分仅为娱乐和引入用&#xff0c;与本文没有太大关联&#xff0c;可以跳过&#xff0c;阅读下面的正文部分。 上篇出现于&#xff1a;DFS&#xff08;深度优先搜索&#xff09;详解&#xff08;概念讲解&#xff0c;图片…

HashSet源码分析

一、HashSet继承关系 1、继承 public boolean equals(Object o) {if (o this)return true;// o没有实现Set接口&#xff0c;返回falseif (!(o instanceof Set))return false;// 向下转换Collection<?> c (Collection<?>) o;// 元素个数不相等&#xff0c;返回f…