【CF643F】Bears and Juice

传送门

• 成都七中外地生考试加强版
• 考虑喝过一桶酒或果汁的熊的集合，集合中的元素是一个二元组 $(x,t_x)$ 表示熊的编号以及第几天和的，注意到如果一只熊和了这桶果汁没倒那么不会再喝，如果倒了那么也不会再和，那么每桶果汁或酒的集合就形如 { ( x , t x ) } \{(x,t_x)\} {(x,tx​)} 其中 $x$ 只出现一次，若两个集合不同则可以区分它们，床的限制是对集合大小的限制，那么容易发现答案即为
  $$an{s}_t=\sum _{i=0}^{\min (n-1,p)}\left(\genfrac{}{}{0ex}{}{n}{i}\right)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;
}