[bzoj2118] 图论SPFA

2118: 墨墨的等式

震惊！这是一道图论题！！

Time Limit: 10 Sec Memory Limit: 259 MB
Submit: 1283 Solved: 496
Description

墨墨突然对等式很感兴趣，他正在研究a1x1+a2x2+…+anxn=B存在非负整数解的条件，他要求你编写一个程序（钦定啊啊），给定N、{an}、以及B的取值范围，求出有多少B可以使等式存在非负整数解。

Input

输入的第一行包含3个正整数，分别表示N、BMin、BMax分别表示数列的长度、B的下界、B的上界。输入的第二行包含N个整数，即数列{an}的值。

Output

输出一个整数，表示有多少b可以使等式存在非负整数解。

Sample Input

2 5 10
3 5
Sample Output

5
HINT

对于100%的数据，N≤12，0≤ai≤5*10^5，1≤BMin≤BMax≤10^12。

先找出最小的ai
设为minn
用一个dis[x]数组表示模minn为x的最小B的值
用adde( i, (a[j] + i)%minn, a[j]);转移边
然后spfa一遍
用求出的dis值，即可给出答案

#include<cstdio>
#include<iostream>
#include<queue>
#include<cstring>
using namespace std;
const int N = 500000 + 5;
const int M = 10000000 + 5;
typedef long long ll;
ll head[N], next[M], to[M], w[M], tail;
ll a[N];
void adde( ll u, ll v, ll w1 ){
    next[++tail] = head[u]; 
    head[u] = tail;
    w[tail] = w1;
    to[tail] = v;
}
ll dis[N], ans, minn, n, bmx, bmn;
bool b[N];
void spfa() {
    queue<ll> q;
    for( ll i = 1; i <= minn; i++ ) dis[i] = 1e15;
    q.push(0);b[0] = 1;
    while( !q.empty() ) {
        ll u = q.front(); q.pop(); b[u] = 0;
        for( ll i = head[u]; i; i = next[i])
            if( dis[to[i]] > dis[u] + w[i] ) {
                dis[to[i]] = dis[u] + w[i];
                if( !b[to[i]] ) {
                    b[to[i]] = 1;
                    q.push(to[i]);
                }
            }
    }
}
int main() {
    scanf("%lld%lld%lld", &n, &bmn, &bmx);
    scanf("%lld", &a[1]); minn = a[1];
    for( int i = 2; i <= n; i++ ) {
        scanf("%lld", &a[i]);
        minn = min( minn, a[i] );
    }
    for( ll i = 0; i < minn; i++ )
        for( int j = 1; j <= n; j++ )
            if( a[j] != minn )
                adde( i, (a[j] + i)%minn, a[j]);
    spfa();
    for( ll i = 0; i < minn; i++ )
        if( dis[i] <= bmx ) {
            ll sd = max( dis[i], bmn ) - 1;
            ll l = sd/minn,r = bmx/minn;
            if((bmx % minn) >= i) r++;
            if((sd % minn) >= i) l++;
            ans += r-l;
        }
    printf("%lld\n", ans);
    return 0;
}