「CodePlus 2018 4 月赛」组合数问题 2

最新推荐文章于 2021-11-09 15:53:28 发布

Freopen

最新推荐文章于 2021-11-09 15:53:28 发布

阅读量341

点赞数

分类专栏：数论奇巧淫技

本文链接：https://blog.csdn.net/qq_35950004/article/details/85635237

版权

数论同时被 2 个专栏收录

157 篇文章 0 订阅

订阅专栏

奇巧淫技

92 篇文章 0 订阅

订阅专栏

博客探讨了在CodePlus 2018年4月竞赛中遇到的一个组合数问题。文章指出当C(a, b)小于C(a, b+1)时，可以利用n个递增队列qi，每个队列的元素对应组合数C(i, j)，通过每次删除队头最大值并重复k次来解决该问题。比较大小的方法是取对数。" 98803408,6961049,Java中的Base64编码详解,"['Java', 'Base64编码', '数据传输', '字符编码', '网络安全']

摘要由CSDN通过智能技术生成

题面

发现 $C (a, b) < C (a, b + 1)$
那么就是有n个递增队列 $q_i$ ，满足 $q_{ij} = C(i,j)$
每次找出队头最大的，然后把队头删去，重复k次。
比大小可以取对数

#include<bits/stdc++.h>
#define maxn 1000006
#define mod 1000000007
#define eps 1e-8
#define LL long long
using namespace std;

LL C[505][505];
int fac[maxn]={1,1},inv[maxn]={1,1},invf[maxn]={1,1},ans;
int num[3000005],cnt,c[3000005];
double lgv[maxn],lgf[maxn],lgn[maxn];
inline bool cmp(const int &u,const int &v){ return lgn[u] > lgn[v] + eps; }

int main()
{	
	C[0][0] = 1;
	for(int i=1;i<505;i++)
	{
		C[i][0] = 1;
		for(int j=1;j<505;j++) C[i][j] = C[i-1][j-1] + C[i-1][j];
	}
	
	for(int i=1;i<505;i++) lgv[i] = log(i) , lgf[i] = lgv[i] + lgf[i-1];
	
	int n,k;
	scanf("%d%d",&n,&k);
	if(n <= 10)
	{
		priority_queue<LL>q;
		for(int i=0;i<=n;i++)
			for(int j=0;j<=i;j++) 
				q.push(C[i][j]);
		for(;k;k--)
			ans = (ans + q.top()) % mod , q.pop();
		printf("%d\n",(ans+mod)%mod);
	}
	else 
	{
		for(int i=2;i<=n;i++)
			fac[i] = 1ll * fac[i-1] * i % mod,
			inv[i] = 1ll * (mod - mod / i) * inv[mod % i] % mod,
			invf[i] = 1ll * invf[i-1] * inv[i] % mod;
		if(k==1) printf("%lld\n",1ll * fac[n] * invf[n/2] % mod * invf[n-(n/2)] % mod);
		else
		{
			for(int i=0;i<=n;i++)
				for(int j=0;j<=i;j++)
					num[++cnt] = C[i][j] , lgn[cnt] = lgf[i] - lgf[j] - lgf[i-j] , c[cnt] = cnt;
			sort(c+1,c+1+cnt,cmp);
			for(int i=1;i<=k;i++)
				ans = (ans + num[c[i]]) % mod;
			printf("%d\n",(ans+mod)%mod);
		}
	}
}