LA 2238 Fixed Partition Memory Management 二分图最佳匹配

LA 2238

题意：输入m个内存区的大小，然后输入n个程序信息，每个程序输入一个k，表示情况总数，然后是按s升序的k对整数 s，t，表示占用内存大小和所花时间，在内存大小为s的内存区时，程序运行时间 ti 满足si<=s<si+1，求最小平均回转时间和调度方案。

思路：每个内存区可设置n个层，倒数第 i 层所贡献的回转时间为 i*t，这样可以把n个程序当做x节点，n*m个内存区层当做y节点构造二分图，把边变成负数求最佳匹配即可，不用担心会出现不合理情况，比如一个内存区有第一层贡献了n*t的回转时间，没有第二层，但是这种情况一定不是最优，最优的一定是合理的。

#include<cstdio>
#include<cstring>
#include<vector>
#include<algorithm>
using namespace std;
const int maxn=505;
const int inf=1e8;
int W[maxn][maxn],cost[maxn][maxn],vis[maxn],a[maxn];
int Lx[maxn],Ly[maxn],left[maxn],slack[maxn];
bool S[maxn],T[maxn];
int n,m;
void init()
{
	for(int i=0;i<maxn;i++)
	for(int j=0;j<maxn;j++)
	W[i][j]=cost[i][j]=inf;
	memset(left,0,sizeof(left));
}
int match(int i)
{
	S[i]=1;
	for(int j=1;j<=n*m;j++)
	{
		if(W[i][j]==-inf)
		continue;
		if(Lx[i]+Ly[j]==W[i][j])
		{
			if(!T[j])
			{
				T[j]=1;
				if(!left[j]||match(left[j]))
				{
					left[j]=i;
					return 1;
				}
			}
		}
		else
		slack[j]=min(slack[j],Lx[i]+Ly[j]-W[i][j]);
	}
	return 0;
}
void update()
{
	int a=inf;
	for(int i=1;i<=n*m;i++)
	if(!T[i])
	a=min(a,slack[i]);
	for(int i=1;i<=n*m;i++)
	{
		if(S[i])Lx[i]-=a;
		if(T[i])Ly[i]+=a;
	}
}
void KM()
{
	memset(Ly,0,sizeof(Ly));
	for(int i=1;i<=m;i++)
	{
		Lx[i]=-inf;
		for(int j=1;j<=n*m;j++)
		if(W[i][j]!=inf)
		Lx[i]=max(Lx[i],W[i][j]);
	}
	for(int i=1;i<=m;i++)
	{
		for(int j=1;j<=n*m;j++)
		slack[j]=inf;
		while(1)
		{
			for(int j=1;j<=n*m;j++)
			S[j]=T[j]=0;
			if(match(i))
			break;
			else
			update();
		}
	}
}
void print()
{
	int t0[100],t1[100],t2[100],ans=0;
	for(int i=1;i<=n;i++)
	{
		int pre=0;
		for(int j=m;j>=1;j--)
		{
			int cur=(i-1)*m+j;
			if(!left[cur])
			continue;
			int p=left[cur];
			t0[p]=i;
			t1[p]=pre;
			t2[p]=pre+cost[p][cur];
			pre=t2[p];
			ans+=t2[p];
		}
	}
	double f=1.0;
	printf("Average turnaround time = %.2lf\n",f*ans/m);
	for(int i=1;i<=m;i++)
	printf("Program %d runs in region %d from %d to %d\n",i,t0[i],t1[i],t2[i]);
	printf("\n");
}
int work(int a,int s[20],int k)
{
	int i;
	for(i=1;i<=k;i++)
	if(s[i]>a)
	break;
	return i-1;
}
int main()
{
	int k,s[20],t[20],kase=0;
	while(~scanf("%d%d",&n,&m)&&n+m)
	{
		init();
		for(int i=1;i<=n;i++)
		scanf("%d",&a[i]);
		for(int i=1;i<=m;i++)
		{
			scanf("%d",&k);
			for(int j=1;j<=k;j++)
			scanf("%d%d",&s[j],&t[j]);
			for(int j=1;j<=n;j++)
			{
				int tmp=work(a[j],s,k);
				if(tmp==0)
				continue;
				for(int p=1;p<=m;p++)
				{
					int cur=(j-1)*m+p;
					cost[i][cur]=t[tmp];
					W[i][cur]=t[tmp]*p;
				}
			}
		}
		for(int i=1;i<=m;i++)
		for(int j=1;j<=n*m;j++)
		W[i][j]=-W[i][j];
		KM();
		printf("Case %d\n",++kase);
		print();
	}
}