http://acm.hdu.edu.cn/showproblem.php?pid=3191

求次短路的长度和个数

相关分析在这里http://blog.csdn.net/u012774187/article/details/40681515

#pragma comment(linker, "/STACK:36777216")

#pragma GCC optimize ("O2")

#include <cstdio>

#include <cstdlib>

#include <cmath>

#include <cstring>

#include <string>

#include <queue>

#include <map>

#include <iostream>

#include <algorithm>

using namespace std;

#define RD(x) scanf("%d",&x)

#define RD2(x,y) scanf("%d%d",&x,&y)

#define RD3(x,y,z) scanf("%d%d%d",&x,&y,&z)

#define clr0(x) memset(x,0,sizeof(x))

#define clr1(x) memset(x,-1,sizeof(x))

#define eps 1e-9

const double pi = acos(-1.0);

typedef long long LL;

typedef unsigned long long ULL;

const int modo = 1e9 + 7;

const int INF = 0x3f3f3f3f;

const int inf = 0x3fffffff;

const LL _inf = 1e18;

const int maxn = 55,maxm = 10005;

struct edge{

    int v,w,next;

    edge(){};

    edge(int vv,int ww,int nnext):v(vv),w(ww),next(nnext){};

}e[maxn*maxn*2];

int head[maxn],inq[maxn][2],dist[maxn][2],cnt[maxn][2];//0最短1次短

int n,m,ecnt;

void init()

{

    clr1(head);

    ecnt = 0;

    for(int i = 1;i <= maxn;++i)

        dist[i][0] = dist[i][1] = inf;

    //fill(dist,dist+maxn*2,inf);

    clr0(inq),clr0(cnt);

}

void add(int u,int v,int w)

{

    e[ecnt] = edge(v,w,head[u]);

    head[u] = ecnt++;

//    e[ecnt] = edge(u,w,head[v]);

//    head[v] = ecnt++;

}

typedef pair<int,int> p2;

struct cmp {

     bool operator() (const p2 &a, const p2 &b)

     {

        if(dist[a.first][a.second] != dist[b.first][b.second])

            return dist[a.first][a.second] > dist[b.first][b.second];

        else

            return a.first > b.first;

     }

};

void spfa(int src,int dst)

{

    priority_queue<p2 , vector<p2> , cmp> q;

    q.push(make_pair(src,0));

    dist[src][0] = 0,cnt[src][0] = 1;

    while(!q.empty()){

        int u = q.top().first,flag = q.top().second;

        q.pop();

        if(inq[u][flag])    continue;

        inq[u][flag] = 1;

        for(int i = head[u];i != -1;i = e[i].next){

            int v = e[i].v,w = e[i].w;

            if(!inq[v][0] && dist[v][0] > dist[u][flag] + e[i].w){

                if(dist[v][0] != inf){

                    dist[v][1] = dist[v][0];

                    cnt[v][1] = cnt[v][0];

                    q.push(make_pair(v,1));

                }

                dist[v][0] = dist[u][flag] + e[i].w;

                cnt[v][0] = cnt[u][flag];

                q.push(make_pair(v,0));

            }else if(!inq[v][0] && dist[v][0] == dist[u][flag] + e[i].w){

                cnt[v][0] += cnt[u][flag];

            }else if(!inq[v][1] && dist[v][1] > dist[u][flag] + e[i].w){

                dist[v][1] = dist[u][flag] + e[i].w;

                cnt[v][1] = cnt[u][flag];

                q.push(make_pair(v,1));

            }else if(!inq[v][1] && dist[v][1] == dist[u][flag] + e[i].w){

                cnt[v][1] += cnt[u][flag];

            }

        }

    }

    printf("%d %d\n",dist[dst][1],cnt[dst][1]);

}

int main(){

    int u,v,w,s,t;

    while(~RD2(n,m)){

        RD2(s,t);s++,t++;

        init();

        while(m--){

            RD3(u,v,w);u++,v++;

            add(u,v,w);

        }

        spfa(s,t);

    }

    return 0;

}

上面是错误代码

经过筛查，发现是 

if(dist[a.first][a.second] != dist[b.first][b.second])

            return dist[a.first][a.second] > dist[b.first][b.second];

有问题，因为不同的路径到达同一个点的同一个状态的值不一样，所有比较和更新穿插在一起就会出错了

可是每次更新dist数组都是向小了更新，所以不应该有影响啊..?

以下是AC代码，谁能清楚一点地告诉为啥上面的程序就不对呢？

#pragma comment(linker, "/STACK:36777216")

#pragma GCC optimize ("O2")

#include <cstdio>

#include <cstdlib>

#include <cmath>

#include <cstring>

#include <string>

#include <queue>

#include <map>

#include <iostream>

#include <algorithm>

using namespace std;

#define RD(x) scanf("%d",&x)

#define RD2(x,y) scanf("%d%d",&x,&y)

#define RD3(x,y,z) scanf("%d%d%d",&x,&y,&z)

#define clr0(x) memset(x,0,sizeof(x))

#define clr1(x) memset(x,-1,sizeof(x))

#define eps 1e-9

const double pi = acos(-1.0);

typedef long long LL;

typedef unsigned long long ULL;

const int modo = 1e9 + 7;

const int INF = 0x3f3f3f3f;

const int inf = 0x3fffffff;

const LL _inf = 1e18;

const int maxn = 55,maxm = 10005;

struct edge{

    int v,w,next;

    edge(){};

    edge(int vv,int ww,int nnext):v(vv),w(ww),next(nnext){};

}e[maxn*maxn*2];

int head[maxn],inq[maxn][2],dist[maxn][2],cnt[maxn][2];//0最短1次短

int n,m,ecnt;

struct node{

    int v,flag,dis;

    node(int a,int b,int c):v(a),flag(b),dis(c){};

    bool operator < (const node & a) const{

        if(a.dis == dis)

            return a.v < v;

        return a.dis < dis;

    }

};

void init()

{

    clr1(head);

    ecnt = 0;

    for(int i = 0;i <= n;++i)

        dist[i][0] = dist[i][1] = inf;

    //fill(dist,dist+maxn*2,inf);

    clr0(inq),clr0(cnt);

}

void add(int u,int v,int w)

{

    e[ecnt] = edge(v,w,head[u]);

    head[u] = ecnt++;

//    e[ecnt] = edge(u,w,head[v]);

//    head[v] = ecnt++;

}

void spfa(int src,int dst)

{

    priority_queue<node> q;

    q.push(node(src,0,0));

    dist[src][0] = 0,cnt[src][0] = 1;

    while(!q.empty()){

        int u = q.top().v,flag = q.top().flag,dis = q.top().dis;

        q.pop();

        if(inq[u][flag])    continue;

        inq[u][flag] = 1;

        for(int i = head[u];i != -1;i = e[i].next){

            int v = e[i].v,w = e[i].w;

            if(!inq[v][0] && dist[v][0] > dis + e[i].w){

                if(dist[v][0] != inf){

                    dist[v][1] = dist[v][0];

                    cnt[v][1] = cnt[v][0];

                    q.push(node(v,1,dist[v][1]));

                }

                dist[v][0] = dis + e[i].w;

                cnt[v][0] = cnt[u][flag];

                q.push(node(v,0,dist[v][0]));

            }else if(!inq[v][0] && dist[v][0] == dis + e[i].w){

                cnt[v][0] += cnt[u][flag];

            }else if(!inq[v][1] && dist[v][1] > dis + e[i].w){

                dist[v][1] = dis + e[i].w;

                cnt[v][1] = cnt[u][flag];

                q.push(node(v,1,dist[v][1]));

            }else if(!inq[v][1] && dist[v][1] == dis+ e[i].w){

                cnt[v][1] += cnt[u][flag];

            }

        }

    }

    printf("%d %d\n",dist[dst][1],cnt[dst][1]);

}

int main(){

    int u,v,w,s,t;

    while(~RD2(n,m)){

        RD2(s,t);s++,t++;

        init();

        while(m--){

            RD3(u,v,w);u++,v++;

            add(u,v,w);

        }

        spfa(s,t);

    }

    return 0;

}

AC代码2

#pragma comment(linker, "/STACK:36777216")

#pragma GCC optimize ("O2")

#include <cstdio>

#include <cstdlib>

#include <cmath>

#include <cstring>

#include <string>

#include <queue>

#include <map>

#include <iostream>

#include <algorithm>

using namespace std;

#define RD(x) scanf("%d",&x)

#define RD2(x,y) scanf("%d%d",&x,&y)

#define RD3(x,y,z) scanf("%d%d%d",&x,&y,&z)

#define clr0(x) memset(x,0,sizeof(x))

#define clr1(x) memset(x,-1,sizeof(x))

#define eps 1e-9

const double pi = acos(-1.0);

typedef long long LL;

typedef unsigned long long ULL;

const int modo = 1e9 + 7;

const int INF = 0x3f3f3f3f;

const int inf = 0x3fffffff;

const LL _inf = 1e18;

const int maxn = 55,maxm = 10005;

struct edge{

    int v,w,next;

    edge(){};

    edge(int vv,int ww,int nnext):v(vv),w(ww),next(nnext){};

}e[maxn*maxn*2];

int head[maxn],inq[maxn][2],dist[maxn][2],cnt[maxn][2];//0最短1次短

int n,m,ecnt;

struct node{

    int v,flag,dis;

    node(int a,int b,int c):v(a),flag(b),dis(c){};

    bool operator < (const node & a) const{

        if(a.dis == dis)

            return a.v < v;

        return a.dis < dis;

    }

};

void init()

{

    clr1(head);

    ecnt = 0;

    for(int i = 0;i <= n;++i)

        dist[i][0] = dist[i][1] = inf;

    //fill(dist,dist+maxn*2,inf);

    clr0(inq),clr0(cnt);

}

void add(int u,int v,int w)

{

    e[ecnt] = edge(v,w,head[u]);

    head[u] = ecnt++;

//    e[ecnt] = edge(u,w,head[v]);

//    head[v] = ecnt++;

}

void spfa(int src,int dst)

{

    priority_queue<node> q;

    q.push(node(src,0,0));

    dist[src][0] = 0,cnt[src][0] = 1;

    while(!q.empty()){

        int u = q.top().v,flag = q.top().flag;//dis = q.top().dis;

        q.pop();

        if(inq[u][flag])    continue;

        inq[u][flag] = 1;

        for(int i = head[u];i != -1;i = e[i].next){

            int v = e[i].v,w = e[i].w;

            if(!inq[v][0] && dist[v][0] > dist[u][flag] + e[i].w){

                if(dist[v][0] != inf){

                    dist[v][1] = dist[v][0];

                    cnt[v][1] = cnt[v][0];

                    q.push(node(v,1,dist[v][1]));

                }

                dist[v][0] = dist[u][flag] + e[i].w;

                cnt[v][0] = cnt[u][flag];

                q.push(node(v,0,dist[v][0]));

            }else if(!inq[v][0] && dist[v][0] == dist[u][flag] + e[i].w){

                cnt[v][0] += cnt[u][flag];

            }else if(!inq[v][1] && dist[v][1] > dist[u][flag] + e[i].w){

                dist[v][1] = dist[u][flag] + e[i].w;

                cnt[v][1] = cnt[u][flag];

                q.push(node(v,1,dist[v][1]));

            }else if(!inq[v][1] && dist[v][1] == dist[u][flag] + e[i].w){

                cnt[v][1] += cnt[u][flag];

            }

        }

    }

    printf("%d %d\n",dist[dst][1],cnt[dst][1]);

}

int main(){

    int u,v,w,s,t;

    while(~RD2(n,m)){

        RD2(s,t);s++,t++;

        init();

        while(m--){

            RD3(u,v,w);u++,v++;

            add(u,v,w);

        }

        spfa(s,t);

    }

    return 0;

}

#pragma comment(linker, "/STACK:36777216")#pragma GCC optimize ("O2")#include <cstdio>#include <cstdlib>#include <cmath>#include <cstring>#include <string>#include
<queue>#include <map>#include <iostream>#include <algorithm>using namespace std;#define RD(x) scanf("%d",&x)#define RD2(x,y) scanf("%d%d",&x,&y)#define RD3(x,y,z) scanf("%d%d%d",&x,&y,&z)#define clr0(x) memset(x,0,sizeof(x))#define clr1(x) memset(x,-1,sizeof(x))#define
eps 1e-9const double pi = acos(-1.0);typedef long long LL;typedef unsigned long long ULL;const int modo = 1e9 + 7;const int INF = 0x3f3f3f3f;const int inf = 0x3fffffff;const LL _inf = 1e18;const int maxn = 55,maxm = 10005;struct edge{ int v,w,next; edge(){};
edge(int vv,int ww,int nnext):v(vv),w(ww),next(nnext){};}e[maxn*maxn*2];int head[maxn],inq[maxn][2],dist[maxn][2],cnt[maxn][2];//0最短1次短int n,m,ecnt;void init(){ clr1(head); ecnt = 0; for(int i = 0;i < maxn;++i) dist[i][0] = dist[i][1] = inf; //fill(dist,dist+maxn*2,inf);
clr0(inq),clr0(cnt);}void add(int u,int v,int w){ e[ecnt] = edge(v,w,head[u]); head[u] = ecnt++;// e[ecnt] = edge(u,w,head[v]);// head[v] = ecnt++;}typedef pair<int,int> p2;struct cmp { bool operator() (const p2 &a, const p2 &b) { if(dist[a.first][a.second]
!= dist[b.first][b.second]) return dist[a.first][a.second] > dist[b.first][b.second]; else return a.first > b.first; }};void spfa(int src,int dst){ priority_queue<p2 , vector<p2> , cmp> q; q.push(make_pair(src,0)); dist[src][0] = 0,cnt[src][0] = 1; while(!q.empty()){
int u = q.top().first,flag = q.top().second; q.pop(); if(inq[u][flag]) continue; inq[u][flag] = 1; for(int i = head[u];i != -1;i = e[i].next){ int v = e[i].v,w = e[i].w; if(!inq[v][0] && dist[v][0] > dist[u][flag] + e[i].w){ if(dist[v][0] != inf){ dist[v][1]
= dist[v][0]; cnt[v][1] = cnt[v][0]; q.push(make_pair(v,1)); } dist[v][0] = dist[u][flag] + e[i].w; cnt[v][0] = cnt[u][flag]; q.push(make_pair(v,0)); }else if(!inq[v][0] && dist[v][0] == dist[u][flag] + e[i].w){ cnt[v][0] += cnt[u][flag]; }else if(!inq[v][1]
&& dist[v][1] > dist[u][flag] + e[i].w){ dist[v][1] = dist[u][flag] + e[i].w; cnt[v][1] = cnt[u][flag]; q.push(make_pair(v,1)); }else if(!inq[v][1] && dist[v][1] == dist[u][flag] + e[i].w){ cnt[v][1] += cnt[u][flag]; } } } printf("%d %d\n",dist[dst][1],cnt[dst][1]);}int
main(){ int u,v,w,s,t; while(~RD2(n,m)){ RD2(s,t);s++,t++; init(); while(m--){ RD3(u,v,w);u++,v++; add(u,v,w); } spfa(s,t); } return 0;}

hdu 3191 次短路的长度和个数的更多相关文章

  1. HDU 3191 次短路长度和条数

    http://www.cnblogs.com/wally/archive/2013/04/16/3024490.html http://blog.csdn.net/me4546/article/det ...

  2. HDU3191-How many paths are there(次短路的长度及其个数)

    oooccc1 is a Software Engineer who has to ride to the work place every Monday through Friday. For a ...

  3. COJ 0579 4020求次短路的长度

    4020求次短路的长度 难度级别:C: 运行时间限制:1000ms: 运行空间限制:51200KB: 代码长度限制:2000000B 试题描述 在一个地图上共有N个路口(编号分别为1到N),R条道路( ...

  4. HDU - 2066 最短路+加一个节点

    一个图上,有M条边,Z个出发点,Y个终止点.求一条最短路,其中起点是Z中的任意一点,终点是Y中任意一点. Input 输入数据有多组,输入直到文件结束. 每组的第一行是三个整数M,Z,Y 接着有M行, ...

  5. ACM: HDU 2544 最短路-Dijkstra算法

    HDU 2544最短路 Time Limit:1000MS     Memory Limit:32768KB     64bit IO Format:%I64d & %I64u Descrip ...

  6. 【最长下降子序列的长度和个数】 poj 1952

    转自http://blog.csdn.net/zhang360896270/article/details/6701589 这题要求最长下降子序列的长度和个数,我们可以增加数组maxlen[size] ...

  7. UESTC 30 &&HDU 2544最短路【Floyd求解裸题】

    最短路 Time Limit: 5000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total Submiss ...

  8. hdu 5521 最短路

    Meeting Time Limit: 12000/6000 MS (Java/Others)    Memory Limit: 262144/262144 K (Java/Others)Total ...

  9. HDU 5726 GCD 区间GCD=k的个数

    GCD Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others)Total Submis ...

随机推荐

  1. cent 7 安装VNC

    1. yum install tigervnc-server; 2.cp /lib/systemd/system/vncserver@.service /etc/systemd/system/vncs ...

  2. NC 6系分配组织方法

    分配客户的方法 String[] customer = { cuVO.getPk_customer() }; ICustAssignService cs = NCLocator.getInstance ...

  3. Java中的四种内部类

    Java中有四种内部类: 成员内部类:定义在另一个类(外部类)的内部,而且与成员属性和方法平级,故称成员内部类.类比于外部类的非静态方法,如果用static修饰就变成了静态内部类 静态内部类:使用st ...

  4. [Hbase]Hbase章１　Hbase框架及基本概念

    Hbase框架介绍 HBase是一个分布式的.面向列的开源数据库. 不同点: l  和一般的关系数据库不同,hbase是一个适合于非结构化数据存储的数据库. l  Hbase是基于列而不是基于行的模式 ...

  5. python数据结构详解

    Python中常见的数据结构可以统称为容器(container).序列(如列表和元组).映射(如字典)以及集合(set)是三类主要的容器. 一.序列(列表.元组和字符串) 序列中的每个元素都有自己的编 ...

  6. Silverlight保存大数据到WCF出现Not Found问题

    Silveright调用WCF保存大数据时,直接在客户端报“ 远程服务器返回了错误:Not Found ”问题,少量是没问题,如下图: 后面在这个文章中找到解决办法: http://www.cnblo ...

  7. Tarjan 割点，桥

    /* ggg ggg ggggggg ggggggg ggggggggggggggggggg ggggggggggggggg ggggggggggg ggggggg ggg g */ /* gyt L ...

  8. vue.js实现购物车功能

    购物车是电商必备的功能,可以让用户一次性购买多个商品,常见的购物车实现方式有如下几种: 1. 用户更新购物车里的商品后,页面自动刷新. 2. 使用局部刷新功能,服务器端返回整个购物车的页面html 3 ...

  9. 52.tableViewCell重用机制避免重复显示问题

    表刷新超出页面显示的内容会重复出现 -(UITableViewCell*)tableView:(UITableView *)tableView cellForRowAtIndexPath:(NSInd ...

  10. SQL SERVER 如果判断text类型数据不为空

    一个字段Remark的数据类型设置先设置为varcharr(255),后来考虑到扩展性需要将其定义为TEXT类型,但是SQL 语句报错.      SQL 语句:      SELECT * FROM ...