HDU 5637 Transform（bfs）

Description
给出n个数的序列a，对于一个整数x，有两种操作：
1.改变x二进制中任一位
2.将x变为x^a[i]，1<=i<=n
m次查询，每次查询输入两个整数x和y，问x最少经过多少次操作可以变成y
Input
第一行一整数T表示用例组数，每组用例首先输入两整数n和m表示序列长度和查询数，之后n个整数ai，最后m行每行两个整数x,y表示一次查询
(T<=20,1<=n<=15,1<=m<=1e5,1<=ai,x,y<=1e5)
Output
对于每组用例，输出sum{i*zi}，zi是第i次查询的答案，结果模1e9+7
Sample Input
1
3 3
1 2 3
3 4
1 2
3 9
Sample Output
10
Solution
每次改变x二进制位任一位其实是将x变成x^(1 << i)，鉴于x的范围，i最多17，考虑到异或的性质，x变成y最多17步，从这n+17个数中选出一些数异或x变成y换句话说也就是从这n+17个数中选出一些数使其异或和为x^y，故做一遍bfs预处理出从这n+17个数选出若干数使其异或和为i的最小步数dis[i]，对于每次查询O(1)输出dis[x^y]即可
Code

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
typedef long long ll;
#define mod 1000000007ll
#define maxn 222222
int T,n,m,a[16],dis[maxn];
queue<int>que;
void bfs()
{
    memset(dis,-1,sizeof(dis));
    while(!que.empty())que.pop();
    que.push(0),dis[0]=0;
    while(!que.empty())
    {
        int x=que.front();
        que.pop();
        for(int i=0;i<17;i++)
        {
            int y=x^(1<<i);
            if(dis[y]==-1)dis[y]=dis[x]+1,que.push(y);
        }
        for(int i=1;i<=n;i++)
        {
            int y=x^a[i];
            if(dis[y]==-1)dis[y]=dis[x]+1,que.push(y);
        }
    }
}
int main()
{
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&n,&m);
        for(int i=1;i<=n;i++)scanf("%d",&a[i]);
        bfs();
        ll ans=0;
        for(int i=1;i<=m;i++)
        {
            int x,y;
            scanf("%d%d",&x,&y);
            ans=(ans+1ll*i*dis[x^y]%mod)%mod;
        }
        printf("%I64d\n",ans);
    }
    return 0;
}