BZOJ ~ 2118 ~ 墨墨的等式 （建图 + 最短路）【经典好题】

思路：我们取{an}序列中的最小值MIN（其实取哪个都无所谓，只不过越小复杂度越低，先往后看一会就知道了）。假设有X%MIN=i，还有一个数字B%MIN=i，且B>=X，那么B一定能通过X加上若干个（>=0）MIN得到。所以我们求出用{an}数组中的数凑出来的最小的X使得X%MIN=i（0<=i<=MIN-1），我们就可以判断出来BMin~BMax中哪些数字可以被凑出来。

我们定义dis[]数组，dis[i]表示用{an}数组中的数凑出来的%MIN=i的最小值。所以现在需要求的就是dis数组了，假设A%MIN=y，(A+a[k])%=z，我们可以看作可以由y->z但需要付出a[k]的代价。所以我们建立0~(MIN-1)个点，每个点i建立n条边（即i到(A+a[k])%MIN的边，边权为a[k]，0<=k<=n-1），求一边0到各个点的最短路就是我们要的dis[]数组。

对于求BMin~BMax中有多少个数字可以被组出来，我们可以先求出从dis[]到BMax有多少个数字可以被组出来然后减去从dis[]到BMin-1有多少个数字可以被组出来就是答案。（类似前缀和的思想）

如果dis[i]小于等于B，那么【dis[i]，B】区间中共有(B-dis[i])/MIN+1个数字可以被组出来（+1是因为dis[i]一定能被组出来）。

这样算正好不重不漏，因为%MIN=1的数字+MIN只能得到%MIN=1的数字，%MIN=2的数字+MIN只能得到%MIN=2的数字。看到这大家也应该明白为啥要选{an}中的最小值了把，因为选越大的值点就越多，边也就越多，复杂度就越高。选哪个值都一样我不知道该怎么解释了，，，

我这篇博客只是解释这个题是怎么写的，而不是怎么想到的，谁知道大神们都是怎么想到的，真TM的NB。

#include<bits/stdc++.h>
using namespace std;
const int MAXN = 5e5 + 5;
const long long INF = 1e12 + 5;
struct Edge
{
    int from, to, dist;       //起点,终点,距离
    Edge(int u, int v, int w):from(u), to(v), dist(w) {}
};

struct Dijkstra
{
    int n, m;                 //结点数,边数(包括反向弧)
    vector<Edge> edges;       //边表。edges[e]和edges[e^1]互为反向弧
    vector<int> G[MAXN];      //邻接表，G[i][j]表示结点i的第j条边在edges数组中的序号
    int vis[MAXN];            //标记数组
    long long d[MAXN];              //s到各个点的最短路
    int pre[MAXN];            //上一条弧

    void init(int n)
    {
        this->n = n;
        edges.clear();
        for (int i = 0; i <= n; i++) G[i].clear();
    }

    void add_edge(int from, int to, int dist)
    {
        edges.push_back(Edge(from, to, dist));
        m = edges.size();
        G[from].push_back(m - 1);
    }

    struct HeapNode
    {
        int from; long long dist;
        bool operator < (const HeapNode& rhs) const
        {
            return rhs.dist < dist;
        }
        HeapNode(int u, long long w): from(u), dist(w) {}
    };

    void dijkstra(int s)
    {
        priority_queue<HeapNode> Q;
        for (int i = 0; i <= n; i++) d[i] = INF;
        d[s] = 0;
        memset(vis, 0, sizeof(vis));
        Q.push(HeapNode(s, 0));
        while (!Q.empty())
        {
            HeapNode x = Q.top(); Q.pop();
            int u = x.from;
            if (vis[u]) continue;
            vis[u] = true;
            for (int i = 0; i < G[u].size(); i++)
            {
                Edge& e = edges[G[u][i]];
                if (d[e.to] > d[u] + e.dist)
                {
                    d[e.to] = d[u] + e.dist;
                    pre[e.to] = G[u][i];
                    Q.push(HeapNode(e.to, d[e.to]));
                }
            }
        }
    }
};

long long n, a[MAXN], MIN, BMin, BMax;
Dijkstra solve;

long long query(long long x)
{
    long long cnt = 0;
    for (int i = 0; i < MIN; i++)
    {
        if (solve.d[i] <= x) cnt += (x - solve.d[i]) / MIN + 1;
    }
    return cnt;
}

int main()
{
    while(~scanf("%d%lld%lld", &n, &BMin, &BMax))
    {
        MIN = INF;
        for (int i = 0; i < n; i++)
        {
            scanf("%d", &a[i]);
            MIN = min(MIN, a[i]);
        }
        solve.init(MIN);
        for (int i = 0; i < MIN; i++)
        {
            for (int j = 0; j < n; j++)
            {
                solve.add_edge(i, (i + a[j]) % MIN, a[j]);
            }
        }
        solve.dijkstra(0);
        printf("%lld\n", query(BMax) - query(BMin - 1));
    }
    return 0;
}
/*
2 5 10
3 5
*/