题目大意：定义点为零维元素，线为一维元素，面为二维元素，空间为三维元素，以此类推，求n维立方体中各维元素都有多少
令f[i][j]为i维立方体内j维元素的个数
考虑n维立方体中的i维元素，将n维立方体拓展至n+1维空间时(觉得抽象的可以想象平面扩展成立方体)
原先的i维元素增加了一倍，同时原先的i-1维元素变为了i维元素
故有f[i][j]=f[i-1][j]*2+f[i-1][j-1]
经过一系列的推导(我不会怎么推，我是打表之后斜着找规律的)，可以得到f[i][j]=2^(i-j)*C(i,j)
然后就有f[n][i]=f[n][i+1]*2*(i+1)/(n-i) 线性求出逆元 从后往前推即可
Update:尼玛谁这么蛋疼加强了数据- -||| 加强后的数据出现了n>=p的情况(其实这个题本来是n<p的,出题人忘记写在题面上了- -)
对当前的值记录p的次数 乘或除一个数的时候对次数进行加减 异或的时候如果次数不为零则异或0
时间复杂度不对啊- - 不管了能卡过去
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define M 10001000
using namespace std;
long long inv[M],ans=1;
int n,p;
struct Long_Long{
long long left;
int pow;
Long_Long(long long _):left(_),pow(0) {}
void operator *= (long long x)
{
while(x%p==0)
x/=p,pow++;
left*=x,left%=p;
}
void operator /= (long long x)
{
while(x%p==0)
x/=p,pow--;
left*=inv[x%p],left%=p;
}
};
void Linear_Shaker()
{
int i;
inv[1]=1;
for(i=2;i<=n&&i<p;i++)
inv[i]=(p-p/i)*inv[p%i]%p;
}
int main()
{
int i;
Long_Long temp=1;
cin>>n>>p;
Linear_Shaker();
for(i=n-1;~i;i--)
{
temp*=static_cast<long long>(i+1<<1);
temp/=static_cast<long long>(n-i);
ans^=temp.pow?0:temp.left;
//cout<<(temp.pow?0:temp.left)<<endl;
}
cout<<ans<<endl;
}
