【组合数学】[HNOI2008][HYSBZ/BZOJ1004]Cards

最新推荐文章于 2019-01-21 19:44:00 发布

outer_form

最新推荐文章于 2019-01-21 19:44:00 发布

阅读量1.6k

点赞数 1

CC 4.0 BY-SA版权

分类专栏：数论文章标签： c++ 组合数学 hnoi bzoj 算法

本文链接：https://blog.youkuaiyun.com/outer_form/article/details/50623347

数论专栏收录该内容

15 篇文章

订阅专栏

本文介绍了一种解决特定卡片染色问题的方法，该问题要求计算不同的染色方案数量，并考虑到通过特定洗牌方法可以获得相同的染色结果。文章详细解释了如何通过组合数学原理计算总的染色方案，并利用费马小定理来处理模运算。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

网上的题解都说用置换，然而我不会，所以，我按照自己的方法写了一份题解。

题目链接

分析

对题目的分析

题目说

可以使用多种洗牌法,而每种方法可以使用多次

然而

输入数据保证任意多次洗牌都可用这 m种洗牌法中的一种代替

说明我们只用考虑仅使用一种洗牌法的情况即可。
然后怎么做呢？搜索？ No!
首先，我们考虑一共有多少种染色的方法。
染色时，我们从n(

n=Sr+Sb+Sg $n=S_r+S_b+S_g$ )张牌中选

Sr $S_r$ 张染成红色，再从剩下的牌中选

Sb $S_b$ 张染成蓝色，剩下的染成绿色，所以根据乘法原理，总的方案数为

S = C (n S r) * C (n - S r S b) * C (n - S r - S b S g) = C (n S r) * C (n - S r S b)

$S=C{n\choose S_r}*C{n-S_r\choose S_b}*C{n-S_r-S_b\choose S_g}=C{n\choose S_r}*C{n-S_r\choose S_b}$
由于

两种染色方法相同当且仅当其中一种可以通过任意的洗牌法洗成另一种

有m种洗牌法，任何一种染色方案都可以根据洗牌法变成其他m种方案，假设这种方案为A。
我们把这m+1中方案看做一组，这一组中的方案通过洗牌能且只能变成其他m种方案，且不会有不属于这一组的方案能通过洗牌变成这一组方案的任何一种。

Why?

对每种洗牌法，都存在一种洗牌法使得能回到原状态。

对于这一组的任何方案，都可以通过洗牌变为A，再变为其他方案，由于输入数据保证任意多次洗牌都可用这 m种洗牌法中的一种代替，所以这一组的方案能通过洗牌变成这一组方案的任何一种。
对于这一组的任何一组方案T，如果能够变为不属于这一组的方案B，那B->T->A，那么B->A，A->B(->是变成的意思)，和假设B不属于这一组矛盾。
如果有不属于这一组的方案B能够变为这m+1种方案之一，证明方法同上。

所以，我们的染色方案被重复计算了m+1次。

a n s = S m + 1 mod p

$ans=\frac S{m+1}\mod p$

One more thing

求解组合数的过程中，可以用longlong存储。
求解S时就不得不mod p了
但是最后一步有除法，怎么办?
乘m+1的逆元就好了，因为p>m+1，且p是质数，所以gcd(m+1,p)=1，由费马小定理可得

(m+1)p−1≡1modp $(m+1)^{p-1}\equiv1\mod p$ ,所以

(m+1)p−2∗(m+1)≡1modp $(m+1)^{p-2}*(m+1)\equiv1\mod p$ ，所以m+1的逆元为

(m+1)p−2 $(m+1)^{p-2}$

代码

#include<cstdio>
int sr,sg,sb,m,p,ans=1,n;
void Read(int &x){
    char c;
    while(c=getchar(),c!=EOF)
        if(c>='0'&&c<='9'){
            x=c-'0';
            while(c=getchar(),c>='0'&&c<='9')
                x=x*10+c-'0';
            ungetc(c,stdin);
            return;
        }
}
int quick_pow(int a,int b){
    int ret=1;
    while(b){
        if(b&1)
            ret=ret*a%p;
        a=a*a%p;
        b>>=1;
    }
    return ret;
}
int main()
{
    Read(sr),Read(sb),Read(sg),Read(m),Read(p);
    n=sr+sb+sg;
    int i;
    long long t;
    for(t=i=1;i<=sr;i++)
        t=t*(n-i+1)/i;
    ans=ans*(t%p)%p;
    n-=sr;
    for(t=i=1;i<=sg;i++)
        t=t*(n-i+1)/i;
    ans=ans*(t%p)%p;
    ans=(ans*quick_pow(m+1,p-2))%p;
    printf("%d",ans);
}