[HNOI/AHOI2018]寻宝游戏

题目大意：
　　$n(nle1000)$个$m(mle5000)$位的二进制数，第$0$个数为$0$。用$wedge$和$vee$将这些数连接起来。$q(qle1000)$次询问，每次给定一个$m$位二进制数$r$，问有多少种连接方案使得结果为$r$。

思路：
　　参考myy的官方题解：

如果第$i$个数之前的运算符是$wedge$，则这一位设为$1$，否则为$0$，得到的二进制数记为$x$。

对每一位分别考虑，对于第$i$位，如果第$j$个数是$1$，那么这一位设为$1$，否则为$0$，得到的二进制数记为$b_i$。

以左边为最低位，按前缀归纳容易证明，第$i$位的结果为$1$，当且仅当$x<b_i$。

我们将$b$从大到小排序，结果设为$c$，那么答案不为零仅当在$c$的顺序下，$r$中没有任何$0$在$1$的前面。找到$r$中第一个$0$的位置，假设是$k$，那么解$x$要满足$c_kle x<c_{k-1}$，于是答案是$c_{k-1}-c_k$。

#include<cstdio>
#include<cctype>
#include<climits>
#include<algorithm>
#include<sys/mman.h>
#include<sys/stat.h>
typedef long long int64;
class MMapInput {
    private:
        char *buf,*p;
        int size;
    public:
        MMapInput() {
            register int fd=fileno(stdin);
            struct stat sb;
            fstat(fd,&sb);
            size=sb.st_size;
            buf=reinterpret_cast<char*>(mmap(0,size,PROT_READ,MAP_PRIVATE,fileno(stdin),0));
            p=buf;
        }
        char getchar() {
            return (p==buf+size||*p==EOF)?EOF:*p++;
    }
};
MMapInput mmi;
inline int getint() {
    register char ch;
    while(!isdigit(ch=mmi.getchar()));
    register int x=ch^'0';
    while(isdigit(ch=mmi.getchar())) x=(((x<<2)+x)<<1)+(ch^'0');
    return x;
}
inline int getdigit() {
    register char ch;
    while(!isdigit(ch=mmi.getchar()));
    return ch^'0';
}
const int N=1001,M=5000,mod=1e9+7;
int pow[M],rank[M],tmp[M],a[M],sum[M];
int main() {
    const int n=getint(),m=getint(),q=getint();
    for(register int i=0;i<m;i++) rank[i]=i;
    for(register int i=pow[0]=1;i<=n;i++) {
        pow[i]=pow[i-1]*2%mod;
    }
    for(register int i=0;i<n;i++) {
        int cnt[2]={-1,m-1};
        for(register int j=0;j<m;j++) {
            if(!(a[j]=getdigit())) cnt[0]++;
            sum[j]=(sum[j]+(int64)a[j]*pow[i])%mod;
        }
        for(register int j=m-1;~j;j--) {
            tmp[cnt[a[rank[j]]]--]=rank[j];
        }
        std::swap(rank,tmp);
    }
    std::reverse(&rank[0],&rank[m]);
    for(register int i=0;i<q;i++) {
        for(register int i=0;i<m;i++) a[i]=getdigit();
        int last1=INT_MIN,first0=INT_MAX;
        for(register int i=m-1;~i;i--) {
            if(a[rank[i]]) {
                last1=i;
                break;
            }
        }
        for(register int i=0;i<m;i++) {
            if(!a[rank[i]]) {
                first0=i;
                break;
            }
        }
        if(last1>first0) {
            puts("0");
            continue;
        }
        const int sum1=last1==INT_MIN?pow[n]:sum[rank[last1]];
        const int sum2=first0==INT_MAX?0:sum[rank[first0]];
        printf("%d
",(sum1-sum2+mod)%mod);
    }
    return 0;
}