题解
感觉很有意思啊。
我们首先要考虑的事情就是，我们的策略应当对于每一种情况都适用，所以说我们不能安排超过 (p) 头熊去喝同一个酒桶中的饮料，同样的，我们也不能使得两个酒桶的饮酒方案相同，不然如果出现了对应此方案的醉倒情况，你无法判断到底是哪一个。
相当于题目就是在问我们有多少种醉倒情况不同的饮酒方案。
我们考虑首先醉倒的人数可以不同，同时醉倒的天数可以不同，醉倒的熊也可以不同。
所以我们对于饮酒天数为 (q) 的情况，我们可以列出下面的式子（在这里熊是可以区分的）。
[R_q=sum_{i=0}^{min(n-1,p)}C_n^icdot q^i
]
然后结合题目的式子，我们最后要求的就是：
[f(q)=igoplus _{i=1}^qi imes(sum_{j=0}^{min(n-1,p)}C_n^jcdot i^j)mod 2^{32}
]
发现暴力就可以了。
代码如下
#include<bits/stdc++.h>
using namespace std;
#define uint unsigned int
#define It map<int,int>::iterator
const int Q=2e6+5,P=135;
int n,p,q,cnt[P];
uint zuhe[P],res=0;
int gcd(int a,int b){return b==0?a:gcd(b,a%b);}
uint ksm(uint x,int k){
uint res=1;
for(;k;k>>=1,x=x*x)
if(k&1) res=res*x;
return res;
}
uint cal(int n,int m){
uint res=1;
for(int i=1;i<=m;++i) cnt[i]=1;
for(int i=n;i>=n-m+1;--i){
int now=i;
for(int j=m;j>=1;--j){
int tmp=gcd(now,j);
while(tmp>1&&cnt[j]){
now/=tmp,cnt[j]--;
cnt[j/tmp]++;
tmp=gcd(now,j);
}
}
res*=(uint)now;
}
// printf("%d
",cnt[1]);
return res;
}
int main(){
cin>>n>>p>>q;
for(int i=0;i<=min(p,n-1);++i) zuhe[i]=cal(n,i);
for(int i=1;i<=q;++i){
uint sum=0;
for(int j=0;j<=min(p,n-1);++j)
sum+=zuhe[j]*ksm(i,j);
sum*=(uint)i,res^=sum;
}
printf("%u
",res);
return 0;
}