BZOJ - 2118 墨墨的等式 (分析 + 最短路)

题目链接

题意：

给出每一个 $a_i$ ，$B$ 的取值范围，让你求这样一个等式 $a_1\ast x_1+a_2\ast y_2+\dots +a_n\ast x_n=B$，有多少组非负整数解。

分析：

乍一看好像很难，无从下手。

那么就百度一下（雾

首先，我们知道，对于一个可以求出来的 $B$，那么肯定会有 $B$ $\%$ $a_i$ = $x$，那么这样子每一个

$B+k\ast a_i$ $（k=1,2...）$

都是可以被构造出来的。

那么，这样我们就可以对于每一个 $x$ 找出最小可能的 $B$，然后就可以统计答案了。

那么找出最小可能的 $B$ 的过程其实就是一个最短路问题。

然后 $a_i$ 取多少呢，当然是最小非 0 的的那个，越小那么 $x$ 就越少，当然其他的 $a_i$ 也是可以的。

在这里我把 $v$ 作为我选的最小那个 $a_i$。

$dis[i]$ 表示 $(a_1\ast x_1+a_2\ast y_2+\dots +a_n\ast x_n)$ $\%$ $v=i$ 的最小路径 。

跑完之后统计答案就是对于每一个 $tag>=dis[i]$

$ans$ $+=(tag-dis[i])/v+1$

$+1$ 是 $dis[i]$ 本身也可以作为一个答案。

代码：

#include <bits/stdc++.h>
#include <ext/rope>
using namespace __gnu_cxx;
using namespace std;
#define mst(a,b) memset(a,b,sizeof(a))
#define ALL(x) x.begin(),x.end()
#define pii pair<int,int>
#define eps 1e-6
#define debug(a) cout << #a": " << a << endl;
#define eularMod(a, b) a < b ? a : a % b + b
inline int lowbit(int x){ return x & -x; }
const int N = 5e5 + 10;
const int mod = (int) 123456789;
const int INF = 0x3f3f3f3f;
const long long LINF = (1LL << 62);
typedef long long LL;
typedef unsigned long long ULL;
const double PI = acos(-1.0);
 
LL dis[N];
int v[20];
bool vis[N];
queue<int> Q;
int n, MI = INF;
 
void spfa () {
    for (;!Q.empty();) Q.pop();
    Q.push(0);
 
    for (;!Q.empty();) {
        int u = Q.front();
        Q.pop();
 
        vis[u] = false;
 
        for (int i = 1; i <= n; i++) {
            int to = (u + v[i]) % MI;
 
            if (dis[to] > dis[u] + v[i]) {
                dis[to] = dis[u] + v[i];
                if (!vis[to]) {
                    vis[to] = true;
                    Q.push(to);
                }
            }
        }
 
    }
}
LL query (LL x) {
    LL ans = 0;
 
    for (int i = 0; i < MI; i++) {
        if (x >= dis[i]) ans += (x - dis[i]) / MI + 1;
    }
 
    return ans;
}
 
int main() {
#ifdef purple_bro
    freopen("in.txt", "r", stdin);
//    freopen("out.txt","w",stdout);
#endif
    LL L, R;
 
    scanf("%d%lld%lld", &n, &L, &R);
 
    for (int i = 1; i <= n; i++) {
        scanf("%d", &v[i]);
        if (v[i] == 0) {
            i--;
            n--;
            continue;
        }
        MI = min(MI, v[i]);
    }
 
    mst(dis, 0x3f);
 
    dis[0] = 0;
 
    spfa();
 
    printf("%lld\n", query(R) - query(L - 1));
 
    return 0;
}