【数学基础】CQYZ_Vijos_P3768 热身赛 Suspect

【问题描述】

　　克拉克今天看了一部电影《非常嫌疑犯》，然后他在梦境中穿越到了电影里。克拉克发现，嫌疑犯们被大卫抓住以后，大卫却感染上了一种细菌。
　　现在，在一个有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;
}