题目链接:点击进入
题目
题意
规定长度为 n 的数组要是满足 a 1 a_1 a1& a 2 a_2 a2&… a n a_n an >= a 1 ⊕ a 2 ⊕ a_1 \oplus a_2 \oplus a1⊕a2⊕… ⊕ a n \oplus a_n ⊕an 则认为是胜利。在每个数组元素小于 2 k 2^k 2k 的前提下,问有多少种数组 a 满足胜利的条件。
思路
看到这个
2
k
2^k
2k ,当时就想到多半按位来看。若是固定每一位,然后看 n 个数这一位的与运算 与 n 个数这一位的异或运算。异或运算有个规律,只要有奇数个 1 那么异或结果为 1 ,否则为 0 。根据这个规律我们大概可以想到,对着一位来说,我们可以选择取偶数个 1 的情况( 0 个 1 、2 个 1 … ( n & 1 ) ? n - 1 : n 个 1 ),这样就满足异或结果为 0 ,但是注意这个时候与运算的值,n 为偶数的情况下与结果可能为 0 也可能为 1 ( 我一开始还以为都是 0 呢,我个傻憨憨 ), n 为奇数的情况下与结果可能为 0 也可能为 1 。所以分类讨论一下。
n 是奇数,那么选取选取偶数个 1 的情况,异或结果是 0 ,与结果也是 0 ,两者相等,还有一种情况就是 n 个数都是 1 ,这时候异或结果是 1 ,与结果也是 1 ,两者依旧相等。因此,单看一位,相等的情况 ( 没有大于只有相等 ) 有
s
u
m
=
C
n
0
+
C
n
2
+
.
.
.
+
C
n
n
−
1
+
1
sum=C_{n}^{0}+C_{n}^{2}+...+C_{n}^{n-1}+1
sum=Cn0+Cn2+...+Cnn−1+1,因为有 k 位,因此方案数为
s
u
m
k
sum^k
sumk
n 是偶数,那么选取选取偶数个 1 的情况,异或结果是 0 ,与结果不一定是 0 也可能是 1 ,因为选取偶数个 1 可能让 n 个数都是 1 的情况,那么这时候与结果为 1 ,否则其它情况与结果为 0 。也就是 1 选取 0 、 2 … n - 2 个的情况,异或结果是 0 ,与结果是 0 , 1 选取 n 个的情况,异或结果是 0 ,与结果是 1 。非全 1 的情况,结果是相等,全 1 的情况是与大于异或。我们可以想一下,若是存在一位与结果大于异或结果的情况,那么这一位大于了,我只要保证比这一位大的那些位是相等的,那么比这一位小的那些位,那不就是随便选了吗,这样最终结果一定是大于。一位上相等的情况共有
s
u
m
=
C
n
0
+
C
n
2
+
.
.
.
+
C
n
n
−
2
sum=C_{n}^{0}+C_{n}^{2}+...+C_{n}^{n-2}
sum=Cn0+Cn2+...+Cnn−2,最终结果大于的情况共有
∑
i
=
1
k
(
2
n
∗
(
i
−
1
)
∗
s
u
m
k
−
i
)
)
\sum_{i=1}^{k}(2^{n*(i-1)}*sum^{k-i}))
∑i=1k(2n∗(i−1)∗sumk−i)) (
s
u
m
k
−
i
sum^{k-i}
sumk−i就是比第 i 位 ( 从低位数第 i 位 ) 大的那些位相等的情况,
2
n
∗
(
i
−
1
)
2^{n*(i-1)}
2n∗(i−1) 就是比第 i 位小的那些位随便选的情况,注意
(
2
m
)
n
=
2
m
n
(2^m)^n=2^{mn}
(2m)n=2mn 而不是等于
2
m
+
n
2^{m+n}
2m+n ,我真是个傻× )。最后在加上相等的情况
s
u
m
k
sum^k
sumk 。
代码
// Problem: Rise of Shadows
// Contest: NowCoder
// URL: https://ac.nowcoder.com/acm/contest/11259/E
// Memory Limit: 524288 MB
// Time Limit: 2000 ms
//
// Powered by CP Editor (https://cpeditor.org)
//#pragma GCC optimize(3)//O3
//#pragma GCC optimize(2)//O2
#include<iostream>
#include<string>
#include<map>
#include<set>
//#include<unordered_map>
#include<queue>
#include<cstdio>
#include<vector>
#include<cstring>
#include<stack>
#include<algorithm>
#include<iomanip>
#include<cmath>
#include<fstream>
#define X first
#define Y second
#define best 131
#define INF 0x3f3f3f3f3f3f3f3f
#define pii pair<int,int>
#define lowbit(x) x & -x
#define inf 0x3f3f3f3f
#define int long long
//#define double long double
//#define rep(i,x,y) for(register int i = x; i <= y;++i)
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const double pai=acos(-1.0);
const int maxn=1e6+10;
const int mod=1e9+7;
const double eps=1e-9;
const int N=5e3+10;
/*--------------------------------------------*/
inline int read()
{
int k = 0, f = 1 ;
char c = getchar() ;
while(!isdigit(c)){if(c == '-') f = -1 ;c = getchar() ;}
while(isdigit(c)) k = (k << 1) + (k << 3) + c - 48 ,c = getchar() ;
return k * f ;
}
/*--------------------------------------------*/
int t,n,m,a[maxn],k;
ll qpow(ll a,ll b)
{
ll ans=1;
while(b)
{
if(b&1) ans=ans*a%mod;
a=a*a%mod;
b>>=1;
}
return ans;
}
ll f[maxn],inv[maxn];
ll init()
{
f[0]=1;
for(int i=1;i<maxn;i++)
f[i]=f[i-1]*i%mod;
inv[maxn-1]=qpow(f[maxn-1],mod-2);
for(int i=maxn-2;i>=0;i--)
inv[i]=inv[i+1]*(i+1)%mod;
return 0;
}
ll c(ll n,ll m)
{
if (n<m) return 0;
return f[n]*inv[m]%mod*inv[n-m]%mod;
}
signed main()
{
// ios::sync_with_stdio(false);
// cin.tie(0);cout.tie(0);
init();
cin>>t;
while(t--)
{
cin>>n>>k;
int tmp=1,len=(n&1)?(n-1):n,sum=0,ans=0;
if(n&1)
{
for(int i=0;i<=len;i+=2)
sum=(sum+c(n,i))%mod;
sum=(sum+1)%mod;
ans=qpow(sum,k)%mod;
}
else
{
for(int i=0;i<=len-2;i+=2)
sum=(sum+c(n,i))%mod;
ans=qpow(sum,k)%mod;
int sum2=0;tmp;
for(int i=1;i<=k;i++)
sum2=(sum2+qpow(2,n*(i-1))%mod*qpow(sum,k-i)%mod)%mod;
ans=(ans+sum2)%mod;
}
cout<<ans<<endl;
}
return 0;
}