题目
题目思路:
建边,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;
}