天天看點

spfa+位運算 2017.10.9 t3

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=;i<=n;i++)
    {
        scanf("%d",&a[i]);
        add(i,a[i]+n,);add(a[i]+n,i,);
    }
    for(int i=;i<=m;i++)
    {
        scanf("%d%d",&x,&y);
        add(x,y,);
    }
    for(int i=;i<=(<<)+n;i++)
    {
        if(i!=) dis[i]=; 
    }
    queue<int>q;
    q.push();
    use[]=;
    while(!q.empty())
    {
        int pc=q.front();
        q.pop();use[pc]=;
        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]]=;
                    q.push(to[i]);
                }
            }
        }
        if(pc>n)
        {
            pc-=n;
            for(int k=;(<<k)<=pc;k++)
            {
                if((pc&(<<k))&&((<<k)<pc))
                {
                    int v=pc^((<<k)+n);//按位依次亦或求出可能的val[j]值并連邊
                    if(dis[v]>dis[pc+n])
                    {
                        dis[v]=dis[pc+n];
                        if(!use[v])
                        {
                            use[v]=;
                            q.push(v);
                        }
                    }
                }
            }
        }
    }
    for(int i=;i<=n;i++)
    {
        if(dis[i]==) printf("-1\n");
        else printf("%d\n",dis[i]);
    }
}