UVALive646 Deranged Exams [容斥+排列组合]

D e r a n g e d   E x a m s Deranged\ Exams Deranged Exams


D e s c r i p t i o n \mathcal{Description} Description
N N N 个位置, 求至少前 K K K 个数字不在自己位置上的方案数 .


正 解 部 分 正解部分
问题转化为: 求前 K K K 个数字存在 在自己位置上的数字 的方案数 .

分为 1 , 2 , 3... K 1,2,3...K 1,2,3...K 个数字在自己位置上的方案,

比如现在要求 i i i 个数字在自己位置上的方案数, 为 C K i ∗ ( N − i ) ! C_{K}^i*(N-i)! CKi(Ni)!,
但是这样会使得 i + 1 i+1 i+1 个数字在自己位置上的方案数多出 1 1 1 倍, 于是减去 C K i + 1 ∗ ( N − i − 1 ) ! C_{K}^{i+1}*(N-i-1)! CKi+1(Ni1)!,
此时又没有了 i + 2 i+2 i+2 个数字在自己位置上的方案数, 于是加上 C K i + 2 ∗ ( N − i − 2 ) ! C_{K}^{i+2}*(N-i-2)! CKi+2(Ni2)!,
以此类推, 奇加偶减 容斥 即可 .

症结 ↑

做了半天发现做错了题!! 本不该做这道题的…


实 现 部 分 实现部分

#include<cstdio>
#define reg register
typedef long long ll;

int Test_id;
int N;
int K;

ll Jc[25];
ll C[25][25];

void Work(){
        scanf("%d%d%d", &Test_id, &N, &K);
        ll Ans = 0;
        for(reg int i = 1; i <= K; i ++){
                if(i & 1) Ans += C[K][i] * Jc[N-i];
                else Ans -= C[K][i] * Jc[N-i];
        }
        printf("%d %lld\n", Test_id, Jc[N] - Ans);
}

int main(){
        C[0][0] = 1;
        for(reg int i = 1; i <= 24; i ++){
                C[i][0] = 1;
                for(reg int j = 1; j <= i; j ++)
                        C[i][j] = C[i-1][j] + C[i-1][j-1];
        }
        Jc[0] = 1;
        for(reg int i = 1; i <= 24; i ++) Jc[i] = i * Jc[i-1];
        int T;
        scanf("%d", &T);
        while(T --) Work();
        return 0;
}

/*
4
1 4 1
2 7 3
3 10 5
4 17 17
*/

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值