数学课

题目描述

wzy又来上数学课了…… 虽然他很菜，但是数学还是懂一丢丢的。老师出了一道题，给定一个包含n个元素的集合P=1,2,3,…,n，求有多少个集合A⊆P，满足任意x∈A有2x∉A，且对于A在P中的补集B，也满足任意x∈B有2x∉B。
wzy花费了1E100天终于算出来了这个答案，但是可恶的caoxia居然又加了一个条件！他要求A的大小恰好为m，这样又有多少个
A呢？这回wzy真的不会了，他找到了你，希望能够得到帮助。由于答案太大，你只需要输出答案mod10000019即可。

题解

我们可以证明对于一个(x,2x)对必须选1个，所以对于一条链我们可以要跳着选。因为每条链的长度不超过log，所以我们可以分块

代码

#include <bits/stdc++.h>
#define I inline
#define LL long long
#define mod 10000019
using namespace std;
LL n,m,q,f,jc[mod],ny[mod];
I LL read(){LL res,f=1; char c; while(!isdigit(c=getchar())) if(c=='-') f=-1; res=(c^48); while(isdigit(c=getchar())) res=(res<<3)+(res<<1)+(c^48); return res*f;}
I LL Pow(LL x,LL k){LL res=1; while(k){if(k&1) res=res*x%mod; x=x*x%mod; k>>=1ll;} return res;}
I LL C(LL a,LL b){if(a<b) return 0; return jc[a]*ny[b]%mod*ny[a-b]%mod;}
I LL lucas(LL a,LL b){if(a<b || b==0) return 1; return C(a%mod,b%mod)*lucas(a/mod,b/mod)%mod;}
I LL g (LL l,LL r){return (r-l+1>>1)+((l&1)&(r&1));}
LL tmp,num,chm;
int main(){
    n=read(); q=read(); chm=g(1,n); 
    jc[0]=1; for(LL i=1;i<mod;i++) jc[i]=jc[i-1]*i%mod;
    ny[mod-1]=Pow(jc[mod-1],mod-2); for(LL i=mod-1;i;i--) ny[i-1]=ny[i]*i%mod;
    for(LL i=1;n;i++){
        LL nx=n>>1,cnt=g(nx+1,n);
        if(i&1) tmp+=cnt; 
        num+=cnt*(i>>1); n=nx;
    }
    chm=Pow(2,chm-tmp);
    while(q--){
        m=read();
        if(m<num || m>num+tmp) puts("0");
        else printf("%lld\n",lucas(tmp,m-num)*chm%mod);
    }
    return 0;
}