2118: 墨墨的等式
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。
n个物品,可以用0到+oo,问[bmin,bmax]区间内有多少价值可以凑出来。
任选一个ai>0,如果一个价值k∗ai+x(0≤x<ai,k≥0)可以被凑出来,那么显然(k+1)∗ai+x,(k+2)∗ai+x,...都可以被凑出来(这样x的范围就是小于ai了)
如果我们对于每个x都找到最小的k满足k∗ai+x可以被凑出来,这个问题就解决了,用SPFA求余d时最小的B,用dis[d]保存,首先队首是0,因为余0时,B为0是肯定可以的,所以dis[0]=0,通过加上不同的a,得到新的余数。
#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, L, R;
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]);
}
}
}
}
ll query(ll x) {
ll ans=0;
for(int i=0;i<minn;++i)
if(dis[i]<=x) ans+=(x-dis[i])/minn+1;
return ans;
}
int main() {
scanf("%lld%lld%lld", &n, &L, &R);
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();
printf("%lld\n", query(R)-query(L-1));
return 0;
}