【CF643F】Bears and Juice

最新推荐文章于 2023-06-24 11:26:03 发布

FSYo

最新推荐文章于 2023-06-24 11:26:03 发布

阅读量245

点赞数

分类专栏：数学

本文链接：https://blog.csdn.net/sslz_fsy/article/details/106661123

版权

数学专栏收录该内容

33 篇文章 1 订阅

订阅专栏

传送门

~~成都七中外地生考试加强版~~
考虑喝过一桶酒或果汁的熊的集合，集合中的元素是一个二元组 $x,t_x)$ 表示熊的编号以及第几天和的，注意到如果一只熊和了这桶果汁没倒那么不会再喝，如果倒了那么也不会再和，那么每桶果汁或酒的集合就形如 ${(x,t_x)\}$ 其中 $x$ 只出现一次，若两个集合不同则可以区分它们，床的限制是对集合大小的限制，那么容易发现答案即为
$ans_t=\sum_{i=0}^{\min (n-1,p)}\binom{n}{i}t^i$
组合数考虑算 $\frac{n^{\underline i}}{i!}$ ，暴力约分，复杂度 $O(p^2\log +qp)$

#include<bits/stdc++.h>
#define cs const
#define pb push_back
using namespace std;
typedef unsigned int ui;
cs int N = 150;
ui q; int n, p, a[N], b[N]; ui c[N];
int main(){
	#ifdef FSYolanda
	freopen("1.in","r",stdin);
	#endif
	scanf("%d%d%d",&n,&p,&q);
	p = min(p,n-1); c[0]=1;
	for(int i=1; i<=p; i++) a[i] = i, b[i] = n-i+1;
	for(int i=1; i<=p; i++){
		for(int j=1; j<=i; j++){
			int t = __gcd(b[i],a[j]);
			a[j]/=t; b[i]/=t;	
		} for(int j=1; j<i; j++){
			int t = __gcd(b[j],a[i]);
			a[i]/=t; b[j]/=t;
		} c[i] = 1;
		for(int j=1; j<=i; j++) c[i] = c[i] * (ui)b[j];
	} ui ans = 0;
	for(int t=1; t<=q; t++){
		ui now = 0, mt = 1;
		for(int i=0; i<=p; i++)
		now += mt * c[i], mt *= (ui)t;
		ans ^= (ui)t * now;
	} cout << ans; return 0;
}