2019 上海icpc网络赛 E Counting Sequences II (生成函数)

最新推荐文章于 2024-06-11 19:36:08 发布

Toniht

最新推荐文章于 2024-06-11 19:36:08 发布

阅读量529

点赞数 1

本文链接：https://blog.csdn.net/qq_41746268/article/details/100858560

版权

题目Counting Sequences II

在这里插入图片描述

刚刚做这题，当成组合数学想了很久，一直做不动。挂机了一个多小时，最后一个小时，突然想到生成函数。早就想比赛做一道生成函数，万万没想到遇见了差点没认出来。
题意
构造长度为n的一个偶数出现偶数次的数列，其中 $a_i\leq m$ ,有多少种。
题解
首先n很大，所以组合数学枚举几个位置放偶数是做不到的。
考虑生成函数
我们设m个数里奇数个数为 $s o$ ，偶数个数为 $s e$
所以
$\begin{aligned} ans &= (1+x+\cfrac{x^2}{2!}+\cdots)^{so}(1+\cfrac{x^2}{2!}+\cfrac{x^4}{4!}+\cdots)^{se}\\ &=e^{x·so}(\cfrac{e^x+e^{-x}}{2})^{se}\\ &=\cfrac{1}{2^{se}}\sum_{i=0}^{se}e^{x·i}e^{-x·(se-i)}e^{x·so}\\ &=\cfrac{1}{2^{se}}\sum_{i=0}^{se}e^{(2·i-se+so)·x}\\ &=\cfrac{1}{2^{se}}\sum_{i=0}^{se}\sum_{n=0}^{\infty}(2·i-se+so)^n\cfrac{x^n}{n!} \end{aligned}$
故第n项系数为
$\cfrac{1}{2^{se}}\sum_{i=0}^{se}(2·i-se+so)^n$
se只有1e5，所以复杂度完全够，n很大，注意一下欧拉降幂即可。

// #include <bits/stdc++.h>
#include <iostream>
#include <algorithm>
#include <cstdio>
#include <queue>
#include <cmath>
#include <string>
#include <cstring>
#include <map>
#include <unordered_map>
#include <set>
#include <vector>
#include <assert.h>
#include <cmath>
#include <ctime>
using namespace std;
#define me(x,y) memset((x),(y),sizeof (x))
#define MIN(x,y) ((x) < (y) ? (x) : (y))
#define MAX(x,y) ((x) > (y) ? (x) : (y))
#define SGN(x) ((x)>0?1:((x)<0?-1:0))
#define ABS(x) ((x)>0?(x):-(x))
// #define int __int128

typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;

const int maxn = 1e5+10;
const int inf = __INT32_MAX__;
const ll INF = __LONG_LONG_MAX__;
const int MOD = 1e9+7;
const double eps = 1e-8;
const double pi = std::acos(-1);
const string cars[] = {"?","?","?"};

ll qpow(ll a,ll b,ll p){
    ll ans = 1;
    while(b){if(b&1) ans = ans*a%p; a=a*a%p; b>>=1;}
    return ans;
}

ll F[maxn],Finv[maxn],inv[maxn];
void init(){
    inv[1] = 1;
    for(ll i = 2; i < maxn; ++i){
        inv[i] = (MOD-MOD/i)*inv[MOD%i]%MOD;
    }
    F[0] = Finv[0] = 1;
    for(ll i = 1; i < maxn; ++i){
        F[i] = F[i-1]*i%MOD;
        Finv[i] = Finv[i-1]*inv[i]%MOD;
    }
}
ll comb(int n,int m){//C[n][m]
    if(m < 0 || m > n) return 0;
    return F[n]*Finv[n-m]%MOD*Finv[m]%MOD;
}

int main(){
    ios::sync_with_stdio(false);
#ifndef ONLINE_JUDGE
    freopen("1in.in","r",stdin);
    freopen("1out.out","w",stdout);
#endif    
    int t;
    init();
    cin>>t;
    while(t--){
        ll n,m;
        cin>>n>>m;
        ll se = m/2,so = m/2+m%2;
        ll inv2 = qpow(2,MOD-2,MOD);
        ll se2 = qpow(inv2,se,MOD);
        ll sum = 0;
        for(int i = 0; i <= se; ++i){
            sum += comb(se,i)*(qpow(2*i-se+so,n%(MOD-1),MOD))%MOD;
            sum%=MOD;
            // cout<<sum<<endl;
        }
        cout<<sum*se2%MOD<<endl;
    }
    

    return 0;
}