今天我在某网站上稍微学习了一下组合数学,准确来讲,今天就看了看组合计数。像一些弱智的排列数、组合数大家肯定在小学奥数就已经精通了(只有我这种蒟蒻忘的精光)。当然,博主比较菜,连二项式定理、帕斯卡恒等式都不会,所以今天只能恶补一下啦。
排列数与组合数大家一定都十分的明白了,所以我在这里也不多说了。
然后二项式定理以及帕斯卡恒等式就是俩公式,自己看看就好了。
今天我做了两道题(貌似特别菜,专适我这种蒟蒻):
1.
解析:这道题看起来十分菜,直接用二项式定理就好了,把px和qy看成一个整体就好。但是写组合数的函数时千万不要直接搞或者用递归,不然。。。就会像我一样超时n次,最好直接用帕斯卡恒等式递推预处理出每组组合数的值(不要忘了边加边取模哦)
代码:
#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
typedef double DB;
const int MOD = 1e4 + 7;
int Max(int a, int b){return a > b ? a : b;}
int Min(int a, int b){return a < b ? a : b;}
int C(int n, int r){
//cout<<n<<" "<<r<<endl;
int ans = 0;
if(r == 0) return 1;
else if(n == r) return 1;
else {
if(r > n / 2 + 1) r = n - r;
ans = ((C(n - 1, r - 1) % MOD) + C(n - 1, r)) % MOD;
return ans;
}
}
int qpow(int l, int r){
int res = 1;
while(r){
if(r & 1) res = res * l % MOD;
l = l * l % MOD;
r >>= 1;
}
return res;
}
int C1(int m, int k){
int count = 0;
for(int i = m; i >= k; i --){
if (k > 1) C1(i - 1, k - 1);
else count ++;
}
return count;
}
inline int read(){
int r = 0, z = 1;
char ch = getchar();
while(ch < '0' || ch > '9'){if(ch == '-') z = -1; ch = getchar();}
while(ch >= '0' && ch <= '9'){r = r * 10 + ch - '0'; ch = getchar();}
return r * z;
}
void fre(){
freopen(".in", "r", stdin);
freopen(".out", "w", stdout);
}
void init(){
int p = read(), q = read(), k = read(), a = read(), b = read();
int ans = ((C1(k, b) * qpow(p, a)) % MOD * qpow(q, b)) % MOD;
int ans1 = qpow(p, a) * qpow(q, b);
cout<<ans<<endl;
}
int main(){
//fre();
init();
return 0;
}
2.
样例输入1:
5 2
1 2 3 4 5
样例输出1:
40
解析:这道题其实也是很简单的,只要仔细一想,就能知道正解。
我们首先先把这一列数从大到小排序,然后和最大的数在一组的k个数,其辣度值就是其本身,如果这k个数中定了这个最大的数,可以在n-1个数中调出k-1个数与这个最大的数组成一组,这样的方案数是
C(n−1,k−1)
,总辣度值就是
C(n−1,k−1)×最大的数a[1]
。
我们再找到次大的数。这是我们排除最大的数,从剩下的n-2个数中调出k-1个数与这个次大的数组成一组,这时所组成的组的辣度值都是这个次大的数这样的总方案数是
C(n−2,k−1)
,总辣度值是
C(n−2,k−1)×次大的数a[2]
。
以此类推,建议中间值开个long long,要不然有可能WA掉。
代码:
#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
typedef double DB;
const int MAXN = 1e5 + 30;
const int MOD = 1e9 + 7;
int c[MAXN][61], a[MAXN];
int Max(int a, int b){return a > b ? a : b;}
int Min(int a, int b){return a < b ? a : b;}
void C(int n, int k){
c[0][0] = 1;
for(int i = 1; i <= n; i ++){
c[i][0] = 1;
for(int j = 1; j <= Min(k, i); j ++)
c[i][j] = (c[i - 1][j] + c[i - 1][j - 1]) % MOD;
}
}
int cmp(int x, int y){return x > y;}
inline int read(){
int r = 0, z = 1;
char ch = getchar();
while(ch < '0' || ch > '9'){if(ch == '-') z = -1; ch = getchar();}
while(ch >= '0' && ch <= '9'){r = r * 10 + ch - '0'; ch = getchar();}
return r * z;
}
void fre(){
freopen(".in", "r", stdin);
freopen(".out", "w", stdout);
}
void init(){
int n = read(), k = read();
LL ans = 0;
C(n, k);
//out<<c[9][9]<<endl;
for(int i = 1; i <= n; i ++) a[i] = read();
sort(a + 1, a + n + 1, cmp);
//cout<<a[1]<<endl;
for(int i = 1; i <= n - k + 1; i ++){
ans = (ans + (LL)c[n - i][k - 1] * a[i]) % (LL)MOD;
}
printf("%lld\n", ans);
}
int main(){
//fre();
init();
return 0;
}