luogu4366 [Code+#4]最短路[优化建边最短路]

显然这里的$n^2$级别的边数不能全建出来，于是盯住xor这个关键点去 瞎猜 探究有没有什么特殊性质可以使得一些边没有必要建出来。

发现一个点经过一次xor $x$，花费$x$这么多代价（先不看$C$），到达另一个点$u ext{xor}x$。

结合异或性质，发现其实这个过程完全可以通过把$x$拆成一位一位去分别xor上$u$，也就是说，任何一个点到达另一个点只需要不断走$2^i$这种xor值就可以到达，于是每个点连出$logn$条边，分别和其序号二进制位每一位异或一个1的数相连。这样，如果要走一条xor路径，就可以拆成走若干条上述简化路径。于是建边就可以得到简化，总边数$m+n ext{log}n$，然后跑最短路即可。。

注意一个RE了无数发的detail：食用上述建边方法需要注意有部分超出$n$但小于$2^{log(n)+1}$的点建的边以及$0$号点连的边也是要考虑进去的，具体为什么自己想。。。于是乎这个数组大小不能照1e5来开了，开两倍2e5。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
#define dbg(x) cerr << #x << " = " << x <<endl
using namespace std;
typedef long long ll;
typedef double db;
typedef pair<ll,int> pii;
template<typename T>inline T _min(T A,T B){return A<B?A:B;}
template<typename T>inline T _max(T A,T B){return A>B?A:B;}
template<typename T>inline char MIN(T&A,T B){return A>B?(A=B,1):0;}
template<typename T>inline char MAX(T&A,T B){return A<B?(A=B,1):0;}
template<typename T>inline void _swap(T&A,T&B){A^=B^=A^=B;}
template<typename T>inline T read(T&x){
    x=0;int f=0;char c;while(!isdigit(c=getchar()))if(c=='-')f=1;
    while(isdigit(c))x=x*10+(c&15),c=getchar();return f?x=-x:x;
}
const int N=1e5+7,M=3000000+7;
struct thxorz{int to,nxt,w;}G[M];
int Head[N<<1],tot;
int n,m,c,len,s,t;
inline void Addedge(int x,int y,int z){G[++tot].to=y,G[tot].nxt=Head[x],Head[x]=tot,G[tot].w=z;}
ll dis[N<<1];
priority_queue<pii,vector<pii>,greater<pii> >q; 
#define y G[j].to
inline void dij(){
    memset(dis,0x3f,sizeof dis);q.push(make_pair(dis[s]=0,s));
    while(!q.empty()){
        ll d=q.top().first;int x=q.top().second;q.pop();
        if(t==x)break;
        if(d^dis[x])continue;
        for(register int j=Head[x];j;j=G[j].nxt)if(MIN(dis[y],d+G[j].w))q.push(make_pair(dis[y],y));
    }
}
#undef y
int main(){//freopen("test.in","r",stdin);//freopen("test.ans","w",stdout);
    read(n),read(m),read(c);
    for(register int i=1,x,y,z;i<=m;++i)read(x),read(y),read(z),Addedge(x,y,z);
    len=__lg(n);read(s),read(t);
    for(register int i=0;i<=(1<<len+1)-1;++i)for(register int j=len;~j;--j)Addedge(i,i^(1<<j),(1<<j)*c);
    dij();
    return printf("%lld
",dis[t]),0;
}

总结：对于边过多的图尝试发掘性质简化建边，去除没有必要的边，用如前缀/异或/线段树等方法来降低边数或者用少量边替代全部情况。