【LibreOJ】#6354. 「CodePlus 2018 4 月赛」最短路 异或优化建图+Dijkstra

【题目】#6354. 「CodePlus 2018 4 月赛」最短路
【题意】给定n个点，m条带权有向边，任意两个点i和j还可以花费(i xor j)*C到达（C是给定的常数），求A到B的最短距离。\(n \leq 10^5,m \leq 5*10^5\)。
【算法】异或优化建图+Dijkstra
正常建边O(n^2)，与其考虑特殊边的处理不如考虑优化n^2的建边方案。一个点x到另一个点y的代价是由每个改变的数位得到的，所以枚举所有点x的每个数位j，从x向\(x \ \ xor \ \ 2^j\)连边代价为\(2^j\)。
复杂度\(O((m+n \ \ log \ \ n) \ \ log \ \ n)\)。

#include<cstdio>
#include<queue>
#include<cstring>
#include<algorithm>
bool isdigit(char c){return c>='0'&&c<='9';}
int read(){
    int s=0,t=1;char c;
    while(!isdigit(c=getchar()))if(c=='-')t=-1;
    do{s=s*10+c-'0';}while(isdigit(c=getchar()));
    return s*t;
}
using namespace std;
const int maxn=100010,maxm=500010;
int tot,first[maxn],n,m,C,d[maxn];
struct edge{int v,w,from;}e[maxm+maxn*20];
void insert(int u,int v,int w){tot++;e[tot].v=v;e[tot].w=w;e[tot].from=first[u];first[u]=tot;}
struct cyc{
    int d,id;
    bool operator < (const cyc &a)const{
        return d>a.d;
    }
};
priority_queue<cyc>Q;
void dijkstra(int S){
    memset(d,0x3f,sizeof(d));
    d[S]=0;Q.push((cyc){0,S});
    while(!Q.empty()){
        cyc X=Q.top();Q.pop();
        int D=X.d,x=X.id;
        if(D!=d[x])continue;
        for(int i=first[x];i;i=e[i].from)if(d[x]+e[i].w<d[e[i].v]){
            d[e[i].v]=d[x]+e[i].w;
            Q.push((cyc){d[e[i].v],e[i].v});
        }
    }
}
int main(){
    n=read();m=read();C=read();
    for(int i=1;i<=m;i++){
        int u=read(),v=read(),w=read();
        insert(u,v,w);
    }
    for(int i=0;i<=n;i++){//
        for(int j=0;j<=18;j++)if((i^(1<<j))<=n){
            insert(i,i^(1<<j),(1<<j)*C);
        }
    }
    int A=read(),B=read();
    dijkstra(A);
    printf("%d",d[B]);
    return 0;
}