Timeline(时间线)

Timeline(时间线)

题目描述

Bessie 在过去的 MM 天内参加了 NN 次挤奶。但她已经忘了她每次挤奶是在哪个时候了。

对于第 ii 次挤奶，Bessie 记得它不早于第 Si天进行。另外，她还有 C条记忆，每条记忆形如一个三元组 (a,b,x)含义是第 b次挤奶在第 a次挤奶结束至少 x天后进行。

现在请你帮 Bessie 算出在满足所有条件的前提下，每次挤奶的最早日期。

保证 Bessie 的记忆没有错误，这意味着一定存在一种合法的方案，使得：

• 第 ii 次挤奶不早于第 Si天进行，且不晚于第 M天进行；
• 所有的记忆都得到满足

题目分析

差分约束

若将每次挤奶看作一个点，间隔时间看作边权，则可以得到一个有向无环图。可以创造一个虚拟原点，向i连边，边权为设定的不早于的天数。题目要求的是最早，但要同时满足三元组和Si两个限制即dist[j] >= max(Si, dist[i] + w)(i为j的前一天)，为了使结果最小取等即可。下面提供两种解题方法

• spfa：跑最长（最短路）的常用算法，不过有被卡的风险
• 拓扑排序：对于DAG来说，可以用拓扑排序跑最长（最短路），时间复杂度较低

code

spfa

#include<bits/stdc++.h>

using namespace std;

const int N = 1e6 + 10, M = 2 * N;
typedef long long ll;

ll n, m, k, t;
ll dist[N];
ll h[N], e[M], ne[M], w[M], idx;
ll q[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 ++;
}

int spfa()
{
    memset(dist, -0x3f, sizeof dist);
    dist[0] = 0;

    queue<int> q;
    q.push(0);
    st[0] = true;

    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];
                if (!st[j])
                {
                    q.push(j);
                    st[j] = true;
                }
            }
        }
    }

}

int main()
{
    memset(h, -1, sizeof h);

    cin >> n >> k >> m;
    for(int i = 1; i <= n; i ++)
    {
        int u;
        cin >> u;
        add(0, i, u);
    }

    for(int i = 0; i < m; i ++)
    {
        int a, b, c;
        cin >> a >> b >> c;
        add(a, b, c);
    }

    spfa();
    for(int i = 1; i <= n; i ++) cout << dist[i] << "\n";

    return 0;
}

拓扑排序

#include<bits/stdc++.h>

using namespace std;

const int N = 1e6 + 10, M = 2 * N;
typedef long long ll;

ll n, m, k, t;
ll dist[N];
ll h[N], e[M], ne[M], w[M], idx;
ll q[N], d[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 topsort()
{
    memset(dist, -0x3f, sizeof dist);
    int hh = 0, tt = 0;
    dist[0] = 0;
    q[0] = 0;

    while(hh <= tt)
    {
        auto t = q[hh ++];
        for(int i = h[t]; i != -1; i = ne[i])
        {
            int j = e[i];
            dist[j] = max(dist[j], dist[t] + w[i]);
            d[j] --;
            if(d[j] == 0)
                q[++ tt] = j;
        }
    }
}

int main()
{
    memset(h, -1, sizeof h);

    cin >> n >> k >> m;
    for(int i = 1; i <= n; i ++)
    {
        int u;
        cin >> u;
        add(0, i, u);
        d[i] ++;
    }

    for(int i = 0; i < m; i ++)
    {
        int a, b, c;
        cin >> a >> b >> c;
        add(a, b, c);
        d[b] ++;
    }

    topsort();
    for(int i = 1; i <= n; i ++) cout << dist[i] << "\n";

    return 0;
}