spfa+位运算 2017.10.9 t3

#include<iostream>
#include<cstdio>
#include<queue>
#define Maxn (1<<20)+200005
#define Maxm 400001
using namespace std;
int n,m,to[Maxn],nxt[Maxn],head[Maxn],e,val[Maxn],a[Maxm],x,y,dis[Maxn];
bool use[Maxm];
void add(int x,int y,int z)
{
    e++;
    to[e]=y;
    nxt[e]=head[x];
    head[x]=e;
    val[e]=z;
}
int main()
{
    freopen("1.txt","r",stdin);
//  freopen("3.txt","w",stdout);
    cin>>n>>m;
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&a[i]);
        add(i,a[i]+n,1);add(a[i]+n,i,0);
    }
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d",&x,&y);
        add(x,y,1);
    }
    for(int i=1;i<=(1<<20)+n;i++)
    {
        if(i!=1) dis[i]=1e8; 
    }
    queue<int>q;
    q.push(1);
    use[1]=1;
    while(!q.empty())
    {
        int pc=q.front();
        q.pop();use[pc]=0;
        for(int i=head[pc];i;i=nxt[i])
        {
            if(dis[to[i]]>dis[pc]+val[i])
            {
                dis[to[i]]=dis[pc]+val[i];
                if(!use[to[i]])
                {
                    use[to[i]]=1;
                    q.push(to[i]);
                }
            }
        }
        if(pc>n)
        {
            pc-=n;
            for(int k=0;(1<<k)<=pc;k++)
            {
                if((pc&(1<<k))&&((1<<k)<pc))
                {
                    int v=pc^((1<<k)+n);//按位依次亦或求出可能的val[j]值并连边
                    if(dis[v]>dis[pc+n])
                    {
                        dis[v]=dis[pc+n];
                        if(!use[v])
                        {
                            use[v]=1;
                            q.push(v);
                        }
                    }
                }
            }
        }
    }
    for(int i=1;i<=n;i++)
    {
        if(dis[i]==1e8) printf("-1\n");
        else printf("%d\n",dis[i]);
    }
}