Moamen and XOR（dp，位运算，组合数）

链接

Moamen and Ezzat are playing a game. They create an array a of $n$ non-negative integers where every element is less than $2k$.

Moamen wins if $a_1\&a_2\&a_3\&\dots \&a_n\ge a_1\oplus a_2\oplus a_3\oplus \dots \oplus a_n$.

Here & denotes the bitwise AND operation, and ⊕ denotes the bitwise XOR operation.

Please calculate the number of winning for Moamen arrays $a$.

As the result may be very large, print the value modulo $1000000007(10^9+7)$.

Input

The first line contains a single integer $t(1\le t\le 5)$— the number of test cases.

Each test case consists of one line containing two integers $n$ and $k(1\le n\le 2\cdot 10^5,0\le k\le 2\cdot 10^5)$.

Output

For each test case, print a single value — the number of different arrays that Moamen wins with.

Print the result modulo $1000000007(10^9+7)$.

Example

input

3
3 1
2 1
4 0

output

5
2
1

Note

In the first example, $n=3,k=1$. As a result, all the possible arrays are $[0,0,0],[0,0,1],[0,1,0],[1,0,0],[1,1,0],[0,1,1],[1,0,1]$, and $[1,1,1]$.

Moamen wins in only $5$ of them: $[0,0,0],[1,1,0],[0,1,1],[1,0,1]$, and $[1,1,1]$.

思路

当n为奇数时，二进制下的每一位，都有 $C_n^0+C_n^2+...+C_n^{n-1}$ 以及全选 $C_n^n$ 这么多种选法。由杨辉三角的性质，上式等于 $2^{n-1}+1$。那么结果就是 $(2^{n-1}+1{)}^k$。

当 $n$ 为偶数时，二进制下的每一位，若是全选，那么 $n$ 个数剩下的 $n\times (k-1)$ 位二进制的任意组合都满足题意；此外，每一位还有 $C_n^0+C_n^2+...+C_n^{n-2}$ 种选择，可化解为 $2^{n-1}-1$。设 $dp[i]$ 表示 $k=1$ 时的选择数，那么$dp[i]=dp[i-1]\times (2^{n-1}-1)+2^{(i-1)\times n}$。

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int mod=1e9+7;
int T,n,k,dp;

int qpow(int a,int b){ int ret=1; while(b){ if(b&1) ret=ret*a%mod; a=a*a%mod,b>>=1; } return ret; }

void solve(){
	cin>>n>>k; dp=1;
	int x=(qpow(2,n-1)+1)%mod,y=(qpow(2,n-1)+mod-1)%mod;
	if(n&1) return (void)(cout<<qpow(x,k)<<"\n");
	for(int i=1;i<=k;i++) dp=(y*dp%mod+qpow(2,(i-1)*n))%mod;
	cout<<dp<<"\n";
}

signed main(){
	ios::sync_with_stdio(false);
	for(cin>>T;T;T--) solve();
}

代码种多次用到 $2$ 的 $i$ 次幂，可以不使用快速幂，而是用一些预处理、递推，可以降低复杂度。

今天写代码的时候，把dp=(y*dp%mod+qpow(2,(i-1)*n))%mod;写成了dp=(y*dp%mod+qpow(2,(i-1)*n%mod))%mod;（ $\text{qpow}$ 中多了个 $\text{mod}$ ）。结果一直 $wa$ ，也不容易找 $bug$ ，后来发现，$2^{11}\text{mod}10\ne 2^{11\text{mod}10}\text{mod}10$ ，把 $\text{mod}$ 去掉就可以了。