codeforces1096Eザ・トップスコアラー

codeforces

/ CY数学演習の組み合わせに基づい

記載された問題の面によれば、困難ではない列挙考えて\(A_1 \)サイズ\(X \)と値\(X \)のカウント数\(Yを\)当選確率を記述することは容易であるように、

\ [\ sum_ {X = R} ^ S \ sum_ {Y = 1} ^ P \ FRAC {1} {Y} \ dbinom {P-1} {Y-1} F(PY、S-X-Y、X) \]

前記\(F(N、M、 X)\) を表し\(N-を\)互いに異なるボックス\(m個\)ボール、各ボックスは、ボールの数未満である\(X \)はプログラムの数です。包含および排除の考慮の古典的な問題は、数に等しいより大きくなければならない\(X \)がある
\ [F(N、M、 X)= \ sum_ {i = 1} ^ N(-1)^ iは\ dbinom {N } {I} G(M- XI、N)\]

\(G(N、M)は、 \) を示している\(m個\)互いに異なるに配置されたボール\(N- \)のためのプログラム・ボックスの数\(\ dbinom {M + N -1} { N-1} \)

最後に、他に例の合計数を忘れないでください\(G(SR、P-)\)いくつかの国境問題に同時に注意を払うで、

#include<iostream>
#include<string.h>
#include<string>
#include<stdio.h>
#include<algorithm>
#include<vector>
#include<math.h>
#include<queue>
#include<set>
#include<map>
using namespace std;
typedef long long ll;
typedef long double db;
typedef pair<int,int> pii;
const int N=5100;
const db pi=acos(-1.0);
#define lowbit(x) (x)&(-x)
#define sqr(x) (x)*(x)
#define rep(i,a,b) for (register int i=a;i<=b;i++)
#define per(i,a,b) for (register int i=a;i>=b;i--)
#define fir first
#define sec second
#define mp(a,b) make_pair(a,b)
#define pb(a) push_back(a)
#define maxd 998244353
#define eps 1e-8
int p,s,r;
ll c[5120][5120];

int read()
{
    int x=0,f=1;char ch=getchar();
    while ((ch<'0') || (ch>'9')) {if (ch=='-') f=-1;ch=getchar();}
    while ((ch>='0') && (ch<='9')) {x=x*10+(ch-'0');ch=getchar();}
    return x*f;
}

void init()
{
    rep(i,0,N)
    {
        c[i][0]=1;
        rep(j,1,i) c[i][j]=(c[i-1][j-1]+c[i-1][j])%maxd;
    }
}

ll calc(int n,int m,int lim)
{
    if ((!n) && (!m)) return 1;
    ll ans=0;
    rep(i,0,n)
    {
        if (i*lim>m) break;
        ll tmp=c[n][i]*c[m-i*lim+n-1][n-1]%maxd;
        if (i&1) ans=(ans+maxd-tmp)%maxd;else ans=(ans+tmp)%maxd;
    }
    return ans;
}

ll qpow(ll x,int y)
{
    ll ans=1;
    while (y)
    {
        if (y&1) ans=ans*x%maxd;
        x=x*x%maxd;y>>=1;
    }
    return ans;
}

int main()
{
    p=read();s=read();r=read();
    init();ll ans=0;
    rep(x,r,s)
    {
        if (x*p<s) continue;
        rep(y,1,p)
        {
            if ((y*x>s) || ((p-y)*(x-1)+y*x<s)) continue;
            ans=(ans+calc(p-y,s-y*x,x)*qpow(y,maxd-2)%maxd*c[p-1][y-1]%maxd)%maxd;
        }
    }
    ans=ans*qpow(c[s-r+p-1][p-1],maxd-2)%maxd;
    printf("%lld",ans);
    return 0;
}

おすすめ

転載: www.cnblogs.com/encodetalker/p/11441240.html