日常训练 Idiot 的方程

题意简述:求以下不定方程非负整数解的组数 : ni=1Xi=m  0Xili,i=1,2,3,,n
数据范围: n20,m1018,lim
一开始毫无想法,因为上限限制只能背包,不能隔板,就是没想到容斥。
问题如果是有下界限制很好解决,所以容斥一发(感觉容斥好像有条件反向的作用),下界限制用隔板即可。
智障的我又忘记+P再取模了。

#include<bits/stdc++.h>
typedef long long ll;
const int N = 25;
const int P = 1234567891;
int n;
ll m,l[N],inv[N],ans;
ll ksm(ll x, int k){
    ll ans = 1;
    for (; k; k>>=1, x = x * x % P)
        if (k & 1) ans = ans * x % P;
    return ans;
}
ll C(ll n, ll k){
    if (n < k) return 0;
    ll ans = 1;
    for (ll i=n; i>n-k; i--)
        ans = ans * (i % P) % P;
    for (int i=k; i>0; i--)
        ans = ans * inv[i] % P;
    return ans;
}
int main(){
    scanf("%d%lld",&n,&m);
    for (int i=1; i<=n; i++)
        scanf("%lld",&l[i]);
    for (int i=1; i<=n; i++)
        inv[i] = ksm(i,P - 2);
    int all = 1 << n;
    for (int S=0; S<all; S++){
        ll temp = m + n - 1, k = 1;
        for (int i=1; i<=n; i++)
            if (S & (1 << (i-1)))
                temp -= l[i] + 1,
                k *= -1;
        ans += C(temp,n-1) * k;
    }
    printf("%lld\n",(ans % P + P) % P);
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值