Codeforces 674F

先令$p=\min (p,n-1)$，这样可以理解为我们至多能损失$p$只熊。
设$F[i][j]$表示还有$i$天，至多还能损失$j$只熊能区分出的最多桶数目，所求即为$F[1][p]$~$F[q][p]$。那么我们考虑一次操作后会发生的情况，我们会将所有桶根据喝了它的熊的情况划分成$2^{n-p+j}$组，显然能立刻知道酒在哪个组中，不过我们只能保留那些喝了它的熊数目不超过$j$个的组（不然会暴毙）。于是有$F[i][j]={\sum }_{k=0}^j\left(\genfrac{}{}{0ex}{}{n-p+j}{k}\right)F[i-1][j-k]$，其中$F[0][i]=1$。
这样直接递推复杂度过高，运用一点常系数线性递推的知识容易做到$\mathcal{O}(p^3+pq)$。这里利用一下转移矩阵的性质，有$F[i][p]=n!\cdot ({\sum }_{j=0}^p\frac{1}{(n-j)!}\cdot [x^i]e^{jx})={\sum }_{j=0}^p\left(\genfrac{}{}{0ex}{}{n}{j}\right)i^j$。于是预处理组合数$\mathrm{m}\mathrm{o}\mathrm{d}{2}^{32}$的值即可简单计算。
时间复杂度$\mathcal{O}(p^2+pq)$。

#include <bits/stdc++.h>
 
using namespace std;
 
typedef unsigned int ui;
 
ui pow_mod(ui x,ui k) {
  ui ans=1;
  while (k) {
  	if (k&1) ans*=x;
  	x*=x;
  	k>>=1;
  }
  return ans;
}
 
ui C(ui n,ui m) {
  ui s1=1,s2=1;
  int cnt=0;
  for(int i=1;i<=m;i++) {
  	ui x=n-i+1;
  	int s=0;
  	while (x%2==0) {
  		x/=2;
  		s++;
	  }
	s1*=x;cnt+=s;
	x=i;s=0;
	while (x%2==0) {
		x/=2;
		s++;
	}
	s2*=x;cnt-=s;
  }
  return s1*pow_mod(s2,(1u<<31)-1)*pow_mod(2,cnt);
}
 
ui Cs[135];
 
int main() {
  int n,m,p;
  scanf("%d%d%d",&n,&p,&m);
  p=min(p,n-1);
  for(int i=0;i<=p;i++) Cs[i]=C(n,i);
  ui ans=0;
  for(int i=1;i<=m;i++) {
  	ui s=0,v=1;
  	for(int j=0;j<=p;j++) {
  		s+=v*Cs[j];
  		v*=i;
	  }
	ans^=(i*s);
  }
  printf("%u\n",ans);
  return 0;
}