Codeforces Round 925 (Div. 3)G. One-Dimensional Puzzle 组合数,隔板法,允许有空的先欠做法

本文主要是介绍Codeforces Round 925 (Div. 3)G. One-Dimensional Puzzle 组合数,隔板法,允许有空的先欠做法,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

Problem - G - Codeforces

当然得先分析题目的情况

类型1,2的数量差不能大于1:

样例3失败解释:

4 6 100 200

第六个型号2是没法再放进去的,所以可以分析出 类型1与类型2之间的差不能超过1

类型3,4可以自己拼自己;

3可以在1 2间,4 在2 1间:

1 2可以拼在一起,也可以在1 2 之间插入若干个3。

4同理。

所以题目问有多少种情况。

当1,2数目不同时:

//  1   2   1
//4   3   4   3  
//  2   1   2
//3   4   3   4

(设a,b,c,d分别是1,2,3,4类型的数目)

可以看出对3和4来说可以插入的位置数目是相同的,为max(a,b)

那么问题等价于若干个盒子,我们把一对球放入盒子中,有多少种放法。(3和4的方法互补干扰,放法之积即为答案。)

而且1 2 之间可以为空,“隔板法”去放3是不行的(隔板法把n个分成若干份,每份至少为1个)。

比如c个球放入a个盒子,这里我们可以先欠做法:

先欠做法:

我们放c+a个球,最后每个盒子取1个,不就把这a个取出来了吗。

这种就等价于c个球放入a个盒子,允许有空

当1,2数目相同时:

// 1 2 1 2
//4 3 4 3 4 
// 2 1 2 1 
//3 4 3 4 3  

要么3多一个,要么4多一个,两种情况加起来即可,同上做法。

代码:

注意阶乘要开到2e6,因为c虽然最大为1e6,但是我们先欠做法计算的是c + gap

ll a[2000005];
ll ksm(int x, int y, int mod) {if (x == 1) return 1;ll res = 1, base = x;while (y) {if (y & 1) res = (res * base) % mod;base = (base * base) % mod;y >>= 1;}return res;
}ll C(ll n, ll m, ll p)
{if (m > n)return 0;return ((a[n] * ksm(a[m], p - 2, p)) % p * ksm(a[n - m], p - 2, p) % p);
}void solve()
{int a, b, c, d;cin >> a >> b >> c >> d;if (abs(a - b) > 1 || a == 0 && b == 0 && c!=0&& d!=0){cout << 0 << endl;return;}if (a == 0&&b==0){if (c == 0 || b == 0)cout << 1 << endl;elsecout << 0 << endl;return;}//int ans = 0;int gap = max(a, b);if (a == b){ans = (C(c + gap, gap, mod) * C(d + gap - 1, gap - 1, mod) % mod + C(c + gap - 1, gap - 1, mod) * C(d + gap, gap, mod) % mod)%mod;}else{ans = C(c + gap - 1, gap-1, mod) * C(d + gap - 1, gap-1, mod) % mod;}cout << ans << endl;
}
signed main()
{ios::sync_with_stdio(false); cin.tie(0); cout.tie(0);a[0] = a[1] = 1;for (int i = 2; i <= 2000000; i++){a[i] = a[i - 1] * i%mod;}int t = 1;cin >> t;while (t--){solve();}return 0;
}

这篇关于Codeforces Round 925 (Div. 3)G. One-Dimensional Puzzle 组合数,隔板法,允许有空的先欠做法的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/726911

相关文章

豆包 MarsCode 不允许你还没有女朋友

在这个喧嚣的世界里,爱意需要被温柔地唤醒。为心爱的她制作每日一句小工具,就像是一场永不落幕的浪漫仪式,每天都在她的心田播撒爱的种子,让她的每一天都充满甜蜜与期待。 背景 在这个瞬息万变的时代,我们都在寻找那些能让我们慢下来,感受生活美好的瞬间。为了让这份浪漫持久而深刻,我们决定为女朋友定制一个每日一句小工具。这个工具会在她意想不到的时刻,为她呈现一句充满爱意的话语,让她的每一天都充满惊喜和感动

hdu4869(逆元+求组合数)

//输入n,m,n表示翻牌的次数,m表示牌的数目,求经过n次操作后共有几种状态#include<iostream>#include<algorithm>#include<cstring>#include<stack>#include<queue>#include<set>#include<map>#include<stdio.h>#include<stdlib.h>#includ

Codeforces Round #240 (Div. 2) E分治算法探究1

Codeforces Round #240 (Div. 2) E  http://codeforces.com/contest/415/problem/E 2^n个数,每次操作将其分成2^q份,对于每一份内部的数进行翻转(逆序),每次操作完后输出操作后新序列的逆序对数。 图一:  划分子问题。 图二: 分而治之,=>  合并 。 图三: 回溯:

Codeforces Round #261 (Div. 2)小记

A  XX注意最后输出满足条件,我也不知道为什么写的这么长。 #define X first#define Y secondvector<pair<int , int> > a ;int can(pair<int , int> c){return -1000 <= c.X && c.X <= 1000&& -1000 <= c.Y && c.Y <= 1000 ;}int m

Codeforces Beta Round #47 C凸包 (最终写法)

题意慢慢看。 typedef long long LL ;int cmp(double x){if(fabs(x) < 1e-8) return 0 ;return x > 0 ? 1 : -1 ;}struct point{double x , y ;point(){}point(double _x , double _y):x(_x) , y(_y){}point op

Codeforces Round #113 (Div. 2) B 判断多边形是否在凸包内

题目点击打开链接 凸多边形A, 多边形B, 判断B是否严格在A内。  注意AB有重点 。  将A,B上的点合在一起求凸包,如果凸包上的点是B的某个点,则B肯定不在A内。 或者说B上的某点在凸包的边上则也说明B不严格在A里面。 这个处理有个巧妙的方法,只需在求凸包的时候, <=  改成< 也就是说凸包一条边上的所有点都重复点都记录在凸包里面了。 另外不能去重点。 int

Codeforces 482B 线段树

求是否存在这样的n个数; m次操作,每次操作就是三个数 l ,r,val          a[l] & a[l+1] &......&a[r] = val 就是区间l---r上的与的值为val 。 也就是意味着区间[L , R] 每个数要执行 | val 操作  最后判断  a[l] & a[l+1] &......&a[r] 是否= val import ja

CSS实现DIV三角形

本文内容收集来自网络 #triangle-up {width: 0;height: 0;border-left: 50px solid transparent;border-right: 50px solid transparent;border-bottom: 100px solid red;} #triangle-down {width: 0;height: 0;bor

Go组合

摘要 golang并非完全面向对象的程序语言,为了实现面向对象的继承这一神奇的功能,golang允许struct间使用匿名引入的方式实现对象属性方法的组合 组合使用注意项 使用匿名引入的方式来组合其他struct 默认优先调用外层方法 可以指定匿名struct以调用内层方法 代码 package mainimport ("fmt")type People struct{}type Pe

组合c(m,n)的计算方法

问题:求解组合数C(n,m),即从n个相同物品中取出m个的方案数,由于结果可能非常大,对结果模10007即可。       共四种方案。ps:注意使用限制。 方案1: 暴力求解,C(n,m)=n*(n-1)*...*(n-m+1)/m!,n<=15 ; int Combination(int n, int m) { const int M = 10007; int