洛谷P4366 [Code+#4]最短路

题目
题目思路：
建边，dijkstra跑最短路，没了

建边
如果按照传统建边方式，将会需要开n²的数组来建邻接表
n最大为100000，显然如果开n²数组将会超空间，这个时候思路虽然很清晰，但是在建边的时候碰到了一个难题，此时我们有以下几种选择：
1.换一种算法（思路），如最小生成树，拓扑排序等，但是深入思考后仍然得不到最优解

2.找规律，这题有一个很明显的特点：异或。所以：
第一我们可以从整个题目的大规律出发，考虑任意两点间距离最小应该为多少（先不考虑有快速通道的情况），可以发现，其实两点的最小距离就是题目所给出的初始距离，但是因为有快速通道的存在，任意两点间的最短距离又不好计算了，所以难以找出大规律；
第二我们可以小规律出发，既刚刚特殊情况下的最短距离，但这还不够，我们继续往下想，为什么是这样的？于是我们可以多走几条其它的路线，能发现有几条路的距离竟然就等于最短路，因为是通过异或来确定距离，我们可以把这种规律往进制上转换，多试几个不同的例子，我们又能找到一个规律，一点到另一点都会通过一条特殊的路线，在经过的点中，每两点都编号转换为2进制都会只有1位不同，（实质上就像是平面向量中基底那样，只不过是每两点都可以由若干个基底表示出来），通过这个规律，我们发现，既然有很多种走法距离都相等，那么我们就可以把一种走法，拆成经过点最多的走法（每两点间编号二进制只有1位不同），接下来仔细想一想。最终，我们发现这个规律能帮助我们撤销掉很多条边，最终边的数量缩小到（nlogn）

对于边的预处理，我们可以通过循环来实现：
枚举点（0也要枚举）
枚举另一点与该点二进制下不同的那一位
建边

然后：愉快地跑dijsktra吧！A了它
关于建边，最好自己拿笔多举几个例子，耐心找规律

#include <bits/stdc++.h>
using namespace std;
const int manx=1e9;
struct edge
{
	int v,w,next;
}a[3000005];
struct dij
{
	int num,len;
	friend bool operator < (dij x,dij y)
	{
		return x.len>y.len;
	}
};
priority_queue<dij> q;
int cnt,n,m,c,head[100005],dis[100005];
bool bj[100005];
int read()
{
	int x=0,f=1;char ch=getchar();
	while (!isdigit(ch)){if (ch=='-')f=-1;ch=getchar();}
	while (isdigit(ch)){x=x*10+ch-'0';ch=getchar();}
	return x*f;
}
void add(int u,int v,int w)
{
	cnt++;
	a[cnt].v=v;
	a[cnt].w=w;
	a[cnt].next=head[u];
	head[u]=cnt;
}
int main()
{
	n=read(),m=read(),c=read();
	for (int i=1;i<=m;i++)
	{
		int x=read(),y=read(),z=read();
		add(x,y,z);
	}
	for (int i=0;i<=n;i++)
		for (int j=1;j<=n;j<<=1)
		{
			if ((i^j)<=n)
			{
				add(i,i^j,c*j);
			}
		}
	for (int i=0;i<=n;i++)
		dis[i]=1e9;
	int s=read(),t=read();
	q.push(dij{s,0});
	dis[s]=0;
	while (!q.empty())
	{
		dij x=q.top();
		q.pop();
		int u=x.num;
		if (bj[u]==0)
		{
			bj[u]=1;
			for (int i=head[u];i;i=a[i].next)
			{
				int v=a[i].v,w=a[i].w;
				if (dis[v]>dis[u]+w)
				{
					dis[v]=dis[u]+w;
					q.push(dij{v,dis[v]});
				}
			}
		}
	}
	printf("%d",dis[t]);
	return 0;
}