【问题描述】
克拉克今天看了一部电影《非常嫌疑犯》,然后他在梦境中穿越到了电影里。克拉克发现,嫌疑犯们被大卫抓住以后,大卫却感染上了一种细菌。
现在,在一个有n个点m条边的加权无向图中,大卫位于1号点,他得走到n号点,一开始大卫身上只有一个细菌,重量为1。当他每走一个单位长度,他身上的每个能分裂的细菌会分裂成两个重量相等的细菌,不过其中一个细菌会失去分裂能力,另一个细菌还能继续分裂。当他到达下一个点的时候,所有他身上的细菌会合成为一个能分裂的细菌,合成的细菌的重量等于合成前所有细菌重量之和。
现在大卫想知道,走到n号点时,自己身上的细菌重量最少是多少。为了方便计算,你只需要给出细菌重量对23333取模后的结果。
【输入格式】
输入的第一行包含两个正整数n,m。接下来m行,每行有三个正整数x,y和w(1<=x,y<=n, w<=1000),点x到点y有一条长度为w的边。保证输入的数据存在一条点1到点n的路径。
【输出格式】
输出一行,表示大卫走到n号点时细菌重量除以23333的余数。
【输入样例】
4 4
1 2 1
1 3 4
2 4 8
3 4 4
【输出样例】
18
【数据范围】
20%的数据,满足:2<=n,m<=10
50%的数据,满足:2<=n,m<=300
100%的数据,满足:2<=n,m<=30000
题解:
题目包装的很好,什么细菌啦分裂啦合成啦,让人看得云里雾里的。这里给出题目大意:
给定n个点,m条带权边的无向图G,求从1到n的最小边权积。
看一下数据规模,m<=30000,w<=1000,最后结果最大为10^90000,写高精度都会爆。。。顺便,如果在求最短路的时候取模的话绝对会WA得死死的
这里需要用到一点数学知识!我们对边的权值取一个log,则有:
log(a)+log(b)==log(a*b)
这不就变成边权和了吗!
最后得到了log(dist[n]),但是不能用exp反解出来(不然我们辛辛苦苦取log有什么意义?),这里标记一下从1到n的最短路径,再沿着路径老老实实求一遍边权积就行了
代码实现:
#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<queue>
using namespace std;
const int maxn=30005,maxm=30005,mo=23333;
const double inf=30000005; //这里的inf是经过粗略估计的值,有兴趣的话可以自己算一算
struct edge{
int to,next,w;
}e[maxm<<1];
int first[maxn]={0},np=0;
double dist[maxn]={0};
int n,m;
bool inq[maxn]={0},vis[maxn]={0};
inline int in(){
int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-') f=-f;ch=getchar();}
while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
return x*f;
}
void addedge(int u,int v,int w){
e[++np]=(edge){v,first[u],w},first[u]=np;
}
void SPFA(int s){
queue<int>q;
q.push(s);
for(int i=1;i<=n;i++) dist[i]=inf;
inq[s]=1,dist[s]=0;
while(!q.empty()){
int i=q.front();q.pop();
inq[i]=0;
for(int p=first[i];p;p=e[p].next){
int j=e[p].to,c=e[p].w;
if(dist[i]+log(c)>=dist[j]) continue;
dist[j]=dist[i]+log(c); //c++里log取的是自然对数
if(inq[j]) continue;
q.push(j);
inq[j]=1;
}
}
}
int main(){
n=in(),m=in();
int u,v,w;
for(int i=1;i<=m;i++){
u=in(),v=in(),w=in()+1;
addedge(u,v,w);
addedge(v,u,w);
}
SPFA(1);
int now=n,ans=1;
vis[now]=1;
while(now!=1){ //找一找1到n的路径
for(int p=first[now];p;p=e[p].next){
int j=e[p].to,c=e[p].w;
if(fabs(dist[j]+log(c)-dist[now])<0.0000001){
now=j,vis[j]=1;
break;
}
}
}
now=1,ans=1;
while(now!=n){
vis[now]=0;
for(int p=first[now];p;p=e[p].next){
int j=e[p].to,c=e[p].w;
if(vis[j]){
ans=ans*c%mo,now=j; //沿着路径计算
break;
}
}
}
printf("%d\n",ans);
return 0;
}