【BZOJ 1375】 [Baltic2002]Bicriterial routing 双调路径

1375: [Baltic2002]Bicriterial routing 双调路径

Time Limit: 5 Sec  Memory Limit: 64 MB

Submit: 299  Solved: 115

[Submit][Status]

Description

来越多，因此选择最佳路径是很现实的问题。城市的道路是双向的，每条道路有固定的旅行时间以及需要支付的费用。路径由连续的道路组成。总时间是各条道路旅行时间的和，总费用是各条道路所支付费用的总和。同样的出发地和目的地，如果路径A比路径B所需时间少且费用低，那么我们说路径A比路径B好。对于某条路径，如果没有其他路径比它好，那么该路径被称为最优双调路径。这样的路径可能不止一条，或者说根本不存在。 给出城市交通网的描述信息，起始点和终点城市，求最优双条路径的条数。城市不超过100个，边数不超过300，每条边上的费用和时间都不超过100。

Input

第一行给出有多少个点,多少条边,开始点及结束点. 下面的数据用于描述这个地图

Output

有多少条最优双调路径

Sample Input

4 5 1 4

2 1 2 1

3 4 3 1

2 3 1 2

3 1 1 4

2 4 2 4

Sample Output

2

HINT

分层图上的最短路。

f[i][j]表示走到i结点，费用为j的最少时间，直接跑spfa即可。

#include <iostream>
#include <algorithm>
#include <cstring>
#include <cstdio>
#include <cstdlib>
#include <queue>
#define inf 0x3f3f3f3f
#define maxc (n-1)*100
using namespace std;
struct edge
{
	int y,ne,c,t;
}e[1000];
struct now
{
	int p,c;
};
int tot=0,n,m;
struct data
{
	int t,f;
}d[105][10005];
int s,t,h[105],inq[105][10005];
void Addedge(int x,int y,int co,int ti)
{
	tot++;
	e[tot].y=y;
	e[tot].ne=h[x];
	h[x]=tot;
	e[tot].c=co;
	e[tot].t=ti;
}
void spfa()
{
	for (int i=1;i<=n;i++)
		for (int j=0;j<=maxc;j++)
			d[i][j].f=0,inq[i][j]=0,d[i][j].t=inf;
	queue<now> q;
	now x;
	x.p=s,x.c=0;
	d[s][0].f=1,d[s][0].t=0;
	inq[s][0]=1;
	q.push(x);
	while (!q.empty())
	{
		x=q.front();
		q.pop();
		inq[x.p][x.c]=0;
                for (int i=h[x.p];i;i=e[i].ne)
		{
			int y=e[i].y;
			int co=e[i].c+x.c;
			if (co>maxc) continue;
			if (d[y][co].t>d[x.p][x.c].t+e[i].t)
			{
				d[y][co].t=d[x.p][x.c].t+e[i].t;
				d[y][co].f=1;
				if (!inq[y][co])
				{
					now aa;
					aa.p=y,aa.c=co;
					q.push(aa),inq[y][co]=1;
				}
			}
		}
	}
}
int main()
{
        scanf("%d%d%d%d",&n,&m,&s,&t);
	for (int i=1;i<=m;i++)
	{
		int x,y,ti,co;
		scanf("%d%d%d%d",&x,&y,&co,&ti);
		Addedge(x,y,co,ti);
		Addedge(y,x,co,ti);
	}
	spfa();
	int ans=0,minn=maxc+10;
        for (int i=0;i<=maxc;i++)
	{
		if (!d[t][i].f) continue;
		if (d[t][i].t>=minn) continue;
		minn=d[t][i].t;
		ans++;
	}
	cout<<ans<<endl;
	return 0;
}

wa是因为求答案时写了break，应该是continue。

其实这个题还可以进行优化：

如果存在f[i][k]，满足k<j且f[i][k]<f[i][j]，那么f[i][j]一定不是最优解。

因此我们可以给每一个结点建一棵线段树，维护最小值即可。

#include <iostream>
#include <algorithm>
#include <cstring>
#include <cstdio>
#include <cstdlib>
#include <queue>
#define inf 0x3f3f3f3f
#define maxc (n-1)*100+5
using namespace std;
struct edge
{
	int y,ne,c,t;
}e[1000];
struct now
{
	int p,c;
};
int tot=0,n,m;
struct data
{
	int t,f;
}d[105][10005];
int s,t,h[105],inq[105][10005];
struct segtree
{
	int l,r,minn;
}T[105][40005];
void Addedge(int x,int y,int co,int ti)
{
	tot++;
	e[tot].y=y;
	e[tot].ne=h[x];
	h[x]=tot;
	e[tot].c=co;
	e[tot].t=ti;
}
void Build(int now,int x,int l,int r)
{
	T[now][x].l=l,T[now][x].r=r;
	if (l==r)
	{
		T[now][x].minn=maxc;
		return;
	}
	int m=(l+r)>>1;
	Build(now,x<<1,l,m);
	Build(now,(x<<1)+1,m+1,r);
	T[now][x].minn=maxc;
}
void Update(int now,int x,int p,int k)
{
	if (T[now][x].l==T[now][x].r&&T[now][x].l==p)
	{
		T[now][x].minn=min(k,T[now][x].minn);
		return;
	}
	int m=(T[now][x].l+T[now][x].r)>>1;
	if (p<=m) Update(now,x<<1,p,k);
	else Update(now,(x<<1)+1,p,k);
	T[now][x].minn=min(T[now][x<<1].minn,T[now][(x<<1)+1].minn);
}
int Get(int now,int x,int l,int r)
{
	if (T[now][x].l>=l&&T[now][x].r<=r)
		return T[now][x].minn;
	int m=(T[now][x].l+T[now][x].r)>>1;
	if (r<=m) return Get(now,x<<1,l,r);
	if (l>m+1) return Get(now,(x<<1)+1,l,r);
	return min(Get(now,x<<1,l,r),Get(now,(x<<1)+1,l,r));
}
void spfa()
{
	for (int i=1;i<=n;i++)
		for (int j=0;j<=maxc;j++)
			d[i][j].f=0,inq[i][j]=0,d[i][j].t=inf;
	queue<now> q;
	now x;
	x.p=s,x.c=0;
	d[s][0].f=1,d[s][0].t=0;
	inq[s][0]=1;
	Update(s,1,0,0);
	q.push(x);
	while (!q.empty())
	{
		x=q.front();
		q.pop();
		inq[x.p][x.c]=0;
        for (int i=h[x.p];i;i=e[i].ne)
		{
			int y=e[i].y;
			int co=e[i].c+x.c;
			if (co>maxc) continue;
			if (d[y][co].t>d[x.p][x.c].t+e[i].t&&Get(y,1,0,co)>d[x.p][x.c].t+e[i].t)
			{
				d[y][co].t=d[x.p][x.c].t+e[i].t;
				d[y][co].f=1;
				Update(y,1,co,d[y][co].t);
				if (!inq[y][co])
				{
					now aa;
					aa.p=y,aa.c=co;
					q.push(aa),inq[y][co]=1;
				}
			}
		}
	}
}
int main()
{
        scanf("%d%d%d%d",&n,&m,&s,&t);
	for (int i=1;i<=m;i++)
	{
		int x,y,ti,co;
		scanf("%d%d%d%d",&x,&y,&co,&ti);
		Addedge(x,y,co,ti);
		Addedge(y,x,co,ti);
	}
	for (int i=1;i<=n;i++)
		Build(i,1,0,maxc);
	spfa();
	int ans=0,minn=maxc+10;
        for (int i=0;i<=maxc;i++)
	{
		if (!d[t][i].f) continue;
		if (d[t][i].t>=minn) continue;
		minn=d[t][i].t;
		ans++;
	}
	cout<<ans<<endl;
	return 0;
}

进入第一版~

其实没必要用线段树，因为是求前缀的最小值，直接树状数组就可以了。。