【BZOJ1820】[JSOI2010]Express Service 快递服务

Description

「飞奔」快递公司成立之后，已经分别与市内许多中小企业公司签订邮件收送服务契约。由于有些公司是在同一栋大楼内，所以「飞奔」公司收件的地点（收件点）最多只有m点 (1, 2, …, m)，因此「飞奔」仅先行采购了三辆货車并聘用了三名司机，每天早上分别从收件地点「1」, 「2」及「3」出发。而在与客户的服务契约中有明确订约：「飞奔」必须在客户提出邮件寄送要求的隔天派人至该公司（地点）收件。为了能更有效率的服务客户并节省收件时间，该公司设立了收件服务登记网站，客户如有邮件需要寄送，必须在需要收件的前一天就先上网登记。为了节省油量，「飞奔」就利用晚上先行安排三位司机隔天的收件路线。每位司机至各地点收件的顺序应与各公司上网登记的顺序相符且必须能在最省油的情况下完成当天所有的收件服务。因此每位司机有可能需要在不同时间重复到同一地点收件，或不同的司机有可能需在不同的时间点前往同一地点收件。如下面范例二(收件公司地点依序为： 4 2 4 1 5 4 3 2 1)所示，虽然司机1一开始就已经在收件地点「1」了，但是他却不能先把后面第四个登记的公司（地点「1」）邮件先收了再前往第一、第二、或第三个登记收件地点（地点「4」,「2」,「4」）收件。但是如果前三个登记收件的服务是由司机2或3來负责，则司机1就可以在地点「1」收了第四个登记的邮件后再前往后面所登记的地点收件。此外，在某些情况下，不一定每辆车都要收到货，也就是說，最佳收件方式也有可能是只需出动一或兩辆车去收货。请写一个程序來帮「飞奔」公司计算每天依预约顺序至各收件地点收件的最少总耗油量。

Input

输入文件第一行有一个整数 m（3 < = m < = 200），代表「飞奔」公司收件的地点数，以1至m之间的整数代号來表示每个地点。接下來的m行（第2到第m+1行），每行有m个整数，代表一个矩阵D。第 i +1行的第 j 个整数是D(i, j)，D(i, j) 代表司机开车从收件点 i 到收件点 j 所需耗油量。最后有一行数串，数串之数字依序为前一天上网登记要求收件的公司地点代号，最多会有1000个收件请求。输入文件中任兩个相邻的整数都以一个空白隔开。 //D(i,j)<=maxint 注意：油量矩阵D满足三角不等式，也就是說 D(i, j) < = D(i, k) + D(k, j)，1 < = i, j, k < = m。因此，每辆车前往下一个收件地点时一定是直接前往，不必先绕道至其它地点再抵达下个收件地点。

Output

输出一个整数，代表收件所需最少总耗油量。

Sample Input

4
0 5 0 6
6 0 5 6
1 6 0 6
1 1 1 0
1 1 1 1 4 4 2 2 2 3

Sample Output

6
样例说明：到每个请求收件地点的司机分别为1 1 1 1 3 3 2 2 2 1，因此司机1只需从起使点1移动到地点3，司机2只需停留在地点2，司机3从起始点3移动到地点4。

题解：一开始以为是个费用流，建图都差不多了（可以从有上下界的费用流考虑，建出来的图长得跟餐巾计划差不多），发现边数有点爆炸。。。

看题解，nm是个暴力！用f[i][a][b][c]表示收完第i个件，1号司机在地点a，2号在b，3号在c的最小总耗油量。显然a,b,c中一定有一个是i件的所在位置，所以可以缩掉一维。并且发现可以用滚动数组，所以又缩掉一维。空间复杂度O(n^2)，时间复杂度O(n^2*m)，然后就卡过了。

#include <cstdio>

#include <cstring>

#include <iostream>

#include <queue>

using namespace std;

typedef long long ll;

int n,m;

int map[210][210],p[1010];

ll f[2][210][210];

ll ans;

int main()

{

	scanf("%d",&n);

	int i,j,k,d,a,b;

	for(i=1;i<=n;i++)	for(j=1;j<=n;j++)	scanf("%d",&map[i][j]);

	memset(f[0],0x3f,sizeof(f[0]));

	p[0]=1,f[0][2][3]=0;

	while(scanf("%d",&p[m+1])!=EOF)	m++;

	for(i=1;i<=m;i++)

	{

		d=i&1;

		memset(f[d],0x3f,sizeof(f[d]));

		for(j=1;j<=n;j++)	for(k=j;k<=n;k++)

		{

			a=min(p[i-1],j),b=max(p[i-1],j),f[d][a][b]=min(f[d][a][b],f[d^1][j][k]+map[k][p[i]]);

			a=min(p[i-1],k),b=max(p[i-1],k),f[d][a][b]=min(f[d][a][b],f[d^1][j][k]+map[j][p[i]]);

			a=min(j,k),b=max(j,k),f[d][a][b]=min(f[d][a][b],f[d^1][j][k]+map[p[i-1]][p[i]]);

		}

	}

	ans=1ll<<60;

	for(i=1;i<=n;i++)	for(j=i;j<=n;j++)	ans=min(ans,f[n&1][i][j]);

	printf("%lld",ans);

	return 0;

}