SPFA: Power Transmission

来源:互联网 发布:mysql中删除列 编辑:程序博客网 时间:2024/05/29 04:40

Power transmission

Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 131072/131072 K (Java/Others)
Total Submission(s): 0    Accepted Submission(s): 0


Problem Description
The project West-East power transmission is famous around the world. It transmits the electricity from western areas to east China. There are many nodes in the power system. Each node is connected with several other nodes in the system by cable. Power can be only transmitted between two connected nodes. For each node, it can’t send power to two or more other nodes at the same time. As we have all known, power will be loss during the transmission. Bob is the chief engineer of the project. He wants to build a transmission line which send power from one node to another node and minimize the power loss at the same time. Now he asks you to help him solve the problem.
 

Input
There are several test cases. For each test case, the first line contains an integer N (0 < N ≤ 50000) which represents the number of nodes in the power system. Then there will be N groups of data following. For the i-th(0 < i ≤ N) group, the first line is an integer ki (ki ≤ 50), which means the node i is connected with ki nodes. The rest of the i-th group data are divided into ki lines. Each line contains an integer ai (0 < ai ≤ N, ai ≠ i) and an integer bi (0 ≤ bi ≤ 100), which represents power can be transmitted from node i to ai and will loss bi% while transmitting. The last line of input data contains three integers separated by single spaces. The first one is s, the second is t (0 < s, t ≤ N), and the third is the total power M (0 < M ≤ 10^6) at node s.
 

Output
For each test case, output the minimum of loss power while transmitting from node s to node t. The result should be printed with two digits to the right of the decimal point. If power cannot be transmitted from node s to node t, output “IMPOSSIBLE!” in a line.
 

Sample Input
422 503 7021 304 2021 104 4001 4 100
 
#include <vector>#include <list>#include <map>#include <set>#include <deque>#include <stack>#include <bitset>#include <iostream>#include <algorithm>#include <cstdio>#include <queue>#include <string>#include <cstring>#include <cmath>#include <cstdlib>#include <ctime>using namespace std;#define INF 1e9int beg, end, m;struct node                 //静态邻接表{    int to,next;          //next为下一条边的编号    double w;} e[3000000];int lastshow[100010];double dis[100010];bool inqueue[100010];int t,n,cnt;queue<int>q;void insert(int a,int b,double w)          //插入a指向b的边,权值为w{    e[cnt].to=b;    e[cnt].w=w;    e[cnt].next=lastshow[a];    lastshow[a]=cnt++;}bool spfa(){    q.push(beg);    while(!q.empty())    {        int x=q.front();        q.pop();        inqueue[x]=false;                               //清除在队列中的标志        int id=lastshow[x];        while(id!=-1)        {            if(dis[e[id].to] < dis[x]*e[id].w)            {                dis[e[id].to] = dis[x] * e[id].w;                if(!inqueue[e[id].to])                  //如果已经在队列中就不要重复加了                {                    inqueue[e[id].to]=true;                    q.push(e[id].to);                }            }            id=e[id].next;        }    }    return true;}void ini(){    for(int i=0; i<=n; ++i){        inqueue[i] = false;        dis[i] = 0;    }    dis[beg] = 1.0;    cnt=0;    while(!q.empty())    {        q.pop();    }}int main(){    int i, j, k;    double ww;    int a, b;    int ki;    double ans;    while(cin >> n){        int a,b;        double w;        for(i = 0 ; i <=n ; i ++)            lastshow[i] = -1;        for(i=1; i<= n; ++i){            scanf("%d", &ki);            for(j = 1; j <= ki; j ++){                scanf("%d", &b);                scanf("%lf", &ww);                w = 1.0 - ww / 100.0;                insert(i,b,w);            }        }        scanf("%d %d %d", &beg, &end, &m);        ini();        spfa();        if(0 == dis[end])            printf("IMPOSSIBLE!\n");        else{            ans = (double)m * (1.0 - dis[end]);            printf("%.2lf\n", ans);        }    }    return 0;}