指数型母函数（Counting Sequences II）

泰勒展开

公式：

常见的公式：

$e^x=1+x+\frac{x^2}{2!}+\frac{x^3}{3!}+...$

$e^{-x}=1-x+\frac{x^2}{2!}-\frac{x^3}{3!}+...$

$e^{kx}=1+kx+\frac{k^2{x}^2}{2!}+\frac{k^3{x}^3}{3!}+...$

$\frac{e^x+e^{-x}}{2}=1+\frac{x^2}{2!}+\frac{x^4}{4!}+\frac{x^6}{6!}...$

$\frac{e^x-e^{-x}}{2}=x+\frac{x^3}{3!}+\frac{x^5}{5!}+\frac{x^7}{7!}...$

---

指数型母函数

普通的母函数解决的问题是组合问题，例如可以选择$[0,x_1]$个$A$，$[0,x_2]$个$B$，求选择共$n$个的方案为式子$(1+x+x^2+...x^{x_1})(1+x^1+x^2+...x^{x_1})$展开后$x^n$的系数。但是这样只能求出组合的方案数，而对于组合并排列就束手无策了。

比方说，题目求的是选择$[0,x_1]$个$A$，$[0,x_2]$个$B$，并排列的方案数。

我们这样构造生成函数：
$$(1+x+\frac{x^2}{2!}+\frac{x^3}{3!}+...+\frac{x^{x_1}}{x_1!})(1+x+\frac{x^2}{2!}+\frac{x^3}{3!}+...+\frac{x^{x_2}}{x_2!})$$

那么展开式中$\frac{x^n}{n!}$的系数就是答案。

---

如果有限项时，逆元加上$FFT$或者直接$double$暴力。考虑无限项时的情况。（当所有数的取值范围大于最大可行取值时，可以视为无限项）

我们可以先将两个式子逆泰勒展开，乘完后再泰勒展开。例如当上述$x_1\to \infty ,x_2\to \infty$时，有：

$$(1+x+\frac{x^2}{2!}+\frac{x^3}{3!}+...)(1+x+\frac{x^2}{2!}+\frac{x^3}{3!}+...)=e^x\ast e^x$$

$$e^{2x}=1+2x+\frac{2^2{x}^2}{2!}+\frac{2^3{x}^3}{3!}+...$$
其中$\frac{x^n}{n!}$的系数为：$2^n$，所以方案数为$2^n$。

---

例题：Counting Sequences II

题意：

构造一个序列$a_{[1,n]}$，满足$a_i\in [1,m]$，其中偶数出现的次数为偶数次。

解析：

假设$2|m$，$\frac{m}{2}$个偶数，有$\frac{m}{2}$个奇数。

构造出母函数：
$$(1+\frac{x^2}{2!}+\frac{x^4}{4!}+\frac{x^6}{6!}...{)}^{\frac{m}{2}}(1+x+\frac{x^2}{2!}+\frac{x^3}{3!}+...{)}^{\frac{m}{2}}$$

$$\to (\frac{e^x+e^{-x}}{2}{)}^{\frac{m}{2}}(e^x{)}^{\frac{m}{2}}$$

$$\to (e^{2x}+1{)}^{\frac{m}{2}}(\frac{1}{2}{)}^{\frac{m}{2}}$$

$$\to (1+\frac{\frac{m}{2}}{1!}{e}^{2x}+\frac{\frac{m}{2}(\frac{m}{2}-1)}{2!}{e}^{4x}+\frac{\frac{m}{2}(\frac{m}{2}-1)(\frac{m}{2}-2)}{3!}{e}^{6x}+...)(\frac{1}{2}{)}^{\frac{m}{2}}$$

而$e^{kx}$中$\frac{x^n}{n!}$的系数为$k^n$，所以整个式子中$\frac{x^n}{n!}$的系数为：
$$(\frac{\frac{m}{2}}{1!}\ast 2^n+\frac{\frac{m}{2}(\frac{m}{2}-1)}{2!}\ast 4^n+\frac{\frac{m}{2}(\frac{m}{2}-1)(\frac{m}{2}-2)}{3!}\ast 6^n+...)(\frac{1}{2}{)}^{\frac{m}{2}}$$

如果$2|̸m$，相当于多一个奇数，也就是$e^x$，那么式子变为：
$$\to (e^x+\frac{\frac{m}{2}}{1!}{e}^{3x}+\frac{\frac{m}{2}(\frac{m}{2}-1)}{2!}{e}^{5x}+\frac{\frac{m}{2}(\frac{m}{2}-1)(\frac{m}{2}-2)}{3!}{e}^{7x}+...)(\frac{1}{2}{)}^{\frac{m}{2}}$$

答案为：
$$(1+\frac{\frac{m}{2}}{1!}\ast 3^n+\frac{\frac{m}{2}(\frac{m}{2}-1)}{2!}\ast 5^n+\frac{\frac{m}{2}(\frac{m}{2}-1)(\frac{m}{2}-2)}{3!}\ast 7^n+...)(\frac{1}{2}{)}^{\frac{m}{2}}$$

代码：

/*
 *  Author : Jk_Chen
 *    Date : 2019-09-15-12.22.19
 */
#include<bits/stdc++.h>
using namespace std;
#define LL long long
#define rep(i,a,b) for(int i=(int)(a);i<=(int)(b);i++)
#define per(i,a,b) for(int i=(int)(a);i>=(int)(b);i--)
#define mmm(a,b) memset(a,b,sizeof(a))
#define pb push_back
#define pill pair<int, int>
#define fi first
#define se second
#define debug(x) cerr<<#x<<" = "<<x<<'\n';
const LL mod=1e9+7;
const int maxn=4e5+9;
LL rd(){ LL ans=0; char last=' ',ch=getchar();
    while(!(ch>='0' && ch<='9'))last=ch,ch=getchar();
    while(ch>='0' && ch<='9')ans=ans*10+ch-'0',ch=getchar();
    if(last=='-')ans=-ans; return ans;
}
/*_________________________________________________________begin*/

LL fac[maxn];
LL ifac[maxn];
LL C(int i,int j){return fac[i]*ifac[j]%mod*ifac[i-j]%mod; }

LL Pow(LL a,LL b,LL mod){
    LL res=1;
    while(b>0){
        if(b&1)res=res*a%mod;
        a=a*a%mod;
        b>>=1;
    }
    return res;
}
/*_________________________________________________________Pow*/

const LL inv2=Pow(2,mod-2,mod);

LL deal(LL n,LL m){
    LL res=0;
    LL st=2;
    if(m%2){
        res=1;
        st=3;
    }
    LL mul=1;
    for(int i=m/2,j=1;i>=1;i--,j++){
        mul=mul*i%mod;
        res=(res+mul*Pow(st,n,mod)%mod*ifac[j]%mod)%mod;
        st+=2;
    }
    return res*Pow(inv2,m/2,mod)%mod;
}

int q[maxn];
int main(){
    fac[0]=1;
    rep(i,1,maxn-1)fac[i]=fac[i-1]*i%mod;
    ifac[maxn-1]=Pow(fac[maxn-1],mod-2,mod);
    per(i,maxn-2,0){
        ifac[i]=ifac[i+1]*(i+1)%mod;
    }

    int t=rd();
    while(t--){
        LL n,m;
        cin>>n>>m;
        printf("%lld\n",deal(n,m));
    }


    return 0;
}