dij+堆优化

写这个dij+堆优化的原因是有些地方卡SPFA，只能搞这个；

香甜的奶油：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<cstdlib>
#include<ctime>
#include<vector>
#include<algorithm>
#include<queue>
using namespace std;
#define LL long long
const int maxn=820;
const int inf=1000000000;
int d[maxn],vis[maxn],n,m,p,b[maxn];
struct node{
    int y,next,v;
}e[5000];
int linkk[maxn],len=0;
void insert(int x,int y,int v){
    e[++len].y=y;
    e[len].next=linkk[x];
    linkk[x]=len;
    e[len].v=v;
}
typedef pair<int,int> pii;
priority_queue<pii,vector<pii>,greater<pii> > q;
int dij(int s){
    while(!q.empty())q.pop();
    q.push(make_pair(0,s));
    memset(d,10,sizeof(vis));d[s]=0;
    memset(vis,0,sizeof(vis));
    while(!q.empty()){
        int v=q.top().first;
        int x=q.top().second;
        q.pop();
        vis[x]=1;
        for(int i=linkk[x];i;i=e[i].next){
            if(vis[e[i].y])continue;
            if(d[e[i].y]>v+e[i].v){
                d[e[i].y]=v+e[i].v;
                q.push(make_pair(d[e[i].y],e[i].y));
            }
        }
    }
    int ans=0;
    for(int i=1;i<=n;i++)ans+=d[b[i]];
    return ans;
}
void init(){
    scanf("%d%d%d",&n,&p,&m);//cow,muchang,daolu
    for(int i=1;i<=n;i++)scanf("%d",&b[i]);
    int x,y,v;
    for(int i=1;i<=m;i++){
        scanf("%d%d%d",&x,&y,&v);
        insert(x,y,v);insert(y,x,v);
    }
}
void work(){
    int minn=inf;
    for(int i=1;i<=p;i++)minn=min(minn,dij(i));
    cout<<minn<<endl;
}
int main(){
    //freopen("1.in","r",stdin);
    //freopen("1.out","w",stdout);
    init();
    work();
}

最后的结果是这种情况下的dij跑得比spfa慢不少；

当然可能是我用了stl，常数太大；

spfa和堆优化dij的唯一区别是一个用队列，需要多次迭代找出最优（也可能直接就找到最优值），一个用堆，可以直接找到一个不确定的点的最优值，但由于用了堆，复杂度要乘一个log；

其他的基本一样的；

spfa稀疏图好用，dij稠密图好用，上面的那道题就是个例子；