Week8 - 图和树的性质与应用(下)

最新推荐文章于 2022-11-12 11:02:53 发布

DoggyWei

最新推荐文章于 2022-11-12 11:02:53 发布

阅读量240

点赞数

本文链接：https://blog.csdn.net/DoggyWei/article/details/105494005

版权

这篇博客探讨了图和树的性质及其在区间选点和猫猫比赛排名问题中的应用。通过差分约束解决区间选点问题，采用SPFA跑最长路算法找到满足条件的最小点数。而在猫猫向前冲的问题中，利用拓扑排序确定猫猫的排名。文章详细介绍了如何构建图并应用SPFA和拓扑排序算法，以解决实际问题。

摘要由CSDN通过智能技术生成

A - 区间选点 II

题意

在这里插入图片描述
input
输入第一行一个整数 n 表示区间的个数，接下来的 n 行，每一行两个用空格隔开的整数 a，b 表示区间的左右端点。1 <= n <= 50000， 0 <= ai <= bi <= 50000 并且 1 <= ci <= bi - ai+1。

output
输出一个整数表示最少选取的点的个数

思路

区间选点，设sum[i]为区间[0,i]上所选点的数目总和，那么区间[a_i,b_i]中的c_i个点就可以转换为sum[b_i] - sum[a_i -1] >= c_i。差分约束公式为 x_i - x_j <= c_k，是小于等于，所以将SPFA中的松弛条件改为

dis[v] < dis[u] + e[i].weight

就能得到，当dis[v]小于dis[u] + e[i].weight的时候，松弛，可以将dis[v]更新为dis[u] + e[i].weight，更新完毕，这时dis[v]就是满足dis[v] >= dis[u] + e[i].weight的一个解，当图被SPFA跑最长路更新完毕，就可以的到一组解，并且是所以满足的解中最小的一组。
要保证解有意义，要满足0 <= sum[i] - sum[i - 1] <= 1 , 由于是跑最长路，所以>= 0的边正常添加就行，然后对右边 <=1取负号，得到sum[i -1] - sum[i] >= -1,插入反向的并且权重为-1的边，插入完毕后让解有意义。
将sum[i]转换为dis[i]，用maxb 来记录区间的最右边的端点，用 mina来记录区间中最左边的端点，用spfa以mina - 1作为起始点跑最长路，得到的dis[maxb]就是答案。

总结

差分约束题，就要看有没有满足x_i - x_j <= c_k的式子，如果这样的式子存在，就能用差分约束，如果题目给的是x_i - x_j => c_k,就跑最长路用dis[v] < dis[u] + e[i].weight松弛得到一组解的下界，如果是x_i - x_j <= c_k，就用dis[v] > dis[u] + e[i].weight的到一组解的上界。
若是在整数域里：x_i - x_j < c_k 可以转化为：x_i - x_j < c_k - 1
若出现了x_i / x_j <= c_k,可以取对数。
在这个题目中，当构建好图时，利用松弛条件dis[v] < dis[u] + e[i].weight，可以知道当dis[v] >= dis[u] + e[i].weight的情况不满足的时候，取到dis[v] = dis[u] + e[i].weight,正好满足条件。
当0到1之间 w=0 松弛的时候条件是dis[u]<dis[v]+e[i].weight 就不会进行下去就WA了，不能取spfa(0)

代码

#include <stdio.h>
#include <cstdio>
#include <queue>
#include <string.h>
using namespace std;

struct edge
{
   //链式前向星
    int to,weight,next;
}e[1000005];

const int inf = 1000000000;
int n,a,b,c,ans,tot;
int head[50005],inq[50005],dis[50005];

void add(int u,int v, int w)
{
   //加边
    ++tot;
    e[tot].next = head[u];
    e[tot].weight = w;
    e[tot].to = v;
    head[u] = tot;
}

void spfa(int s)
{
   //最长路
    for (int i = 0 ; i < 50005 ; i++)
    {
   
        dis[i] = -inf;
        inq[i] = 0;
    }
    queue<int> q;
    q.push(s);
    dis[s] = 0;
    inq[s] = 1;
    while (!q.empty())
    {
   
        int u = q.front();
        q.pop();
        inq[u] = 0;
        for(int i = head[u] ; i != 0 ; i = e[i].next)
        {
   
            int v = e[i].to;
            if(dis[v] < dis[u] + e[i].weight)
            {
   //变换松弛操作 找到dis[v] >= dis[u] + e[i].weight
                dis[v] = dis[u] + e[i].weight;
                if(inq[v] == 0)
                {
   
                    q.push(v);
                    inq[v] = 1;
                }
            }
        }
    }

}

int main()
{
   
    scanf("%d",&n);
    int maxb = 0;
    int mina = inf;
    for(int t = 0 ; t < n ; t++)
    {
   
        scanf("%d%d%d",&a,&b,&c);
        add