题意:输入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();
}
}