Bzoj2118: 墨墨的等式

题面

传送门

Sol

和这道题很像Luogu跳楼机
找到最小的那个 a[i] a [ i ] 设为 mn m n
SPFA S P F A 算出得到每个值为 i i 的是否能凑出来,设为f[i]

实际上 i i 变成i%mn,所有的这一类都可以通过 mn m n 这个桥梁得到
那么直接求到 i%mn i % m n 就好了

最后统计 [0,Bmax] [ 0 , B m a x ] 时就是所有的

Bmaxf[i]mn+1 ⌊ B m a x − f [ i ] m n ⌋ + 1

减去 Bmin B m i n 的就好了

# include <bits/stdc++.h>
# define RG register
# define IL inline
# define Fill(a, b) memset(a, b, sizeof(a))
using namespace std;
const int _(5e5 + 5);
typedef long long ll;

IL ll Input(){
    RG ll x = 0, z = 1; RG char c = getchar();
    for(; c < '0' || c > '9'; c = getchar()) z = c == '-' ? -1 : 1;
    for(; c >= '0' && c <= '9'; c = getchar()) x = (x << 1) + (x << 3) + (c ^ 48);
    return x * z;
}

int a[13], first[_], cnt, vis[_], mn = 2e9, n;
ll L, R, f[_];
struct Edge{
    int to, next, w;
} edge[_ * 15];
queue <int> Q;

IL void Add(RG int u, RG int v, RG int w){
    edge[cnt] = (Edge){v, first[u], w}, first[u] = cnt++;
}

IL ll Calc(RG ll x){
    RG ll ret = 0;
    for(RG int i = 0; i < mn; ++i) if(f[i] <= x) ret += (x - f[i]) / mn + 1;
    return ret;
}

int main(RG int argc, RG char* argv[]){
    n = Input(), L = Input(), R = Input();
    for(RG int i = 1; i <= n; ++i) a[i] = Input(), mn = min(mn, a[i]);
    for(RG int i = 0; i < mn; ++i) first[i] = -1;
    for(RG int i = 0; i < mn; ++i)
        for(RG int j = 1; j <= n; ++j) Add(i, (i + a[j]) % mn, a[j]);
    Fill(f, 127), f[0] = 0, vis[0] = 1, Q.push(0);
    while(!Q.empty()){
        RG int u = Q.front(); Q.pop();
        for(RG int e = first[u]; e != -1; e = edge[e].next){
            RG int v = edge[e].to, w = edge[e].w;
            if(f[u] + w < f[v]){
                f[v] = f[u] + w;
                if(!vis[v]) vis[v] = 1, Q.push(v);
            }
        }
        vis[u] = 0;
    }
    return printf("%lld\n", Calc(R) - Calc(L - 1)), 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值