大意:给出一些工件,每个工件完成需要一些时间。所有的工件必须在特定的时间内完成。有无数台机器可以加工工件,加工两件工件之间必须间隔一小时,每台机器的费用等于总时间与在这台机器上加工的最后一个工件结束的时间的差的平方。求在不超过时间的情况下,把所有的工件加工完成需要花费的最小费用。
【题解】
一看就知道是dp题,但是没大有思路,百爪挠心。。。
设f(i,j)表示前i个工件且最后一台机器所用的时间是j的最优解,那么我们面对新的一个工件,只有两种选择:1、选择新开一台机器,那么这时就应该是f[i][a[i]]=min(f[i][a[i]],f[i-1][j]+(t-j)*(t-j)),代表的是上一个物品的最优值加上上一台机器空出来的部分的平方;2、选择用这台机器继续工作,那么这时就应该是f[i][j]=min(f[i][j],f[i-1][j-a[i]-1]),代表的只是上一个工件的最优值;就这样进行dp,需要注意的是最后算ans的时候要把最后一台机器的空出来的部分加进去。
输出的话,可以先看一下最优的方案的最后一台机器用了多少,然后倒着往前推:如果上一个工件的时间小于用了的时间,那么这台机器上一定还有一台机器在它前面;如果等于的话,说明这台机器上只有它一个,那么枚举一下上一台机器看看上一台机器用了多长时间,以此类推。
【代码】
来自BeiYu的代码
#include<cstdio>
#include<cstring>
#include<vector>
#include<iostream>
using namespace std;
#define N 505
#define M 30005
int t,n,ans,m,tot;int a[N],f[N][M];vector<int> k[N];
int in(){
int x=0;char ch=getchar();
while(ch>'9'||ch<'0') ch=getchar();
while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
return x;
}
void find(int i,int w){
if(i==0) return;
if(a[i]<w){
k[tot].push_back(a[i]);
find(i-1,w-a[i]-1);
}
else if(a[i]==w){
k[tot].push_back(a[i]);tot++;
for(int j=0;j<=t;j++)
if(f[i-1][j]+(t-j)*(t-j)==f[i][w]){
find(i-1,j);break;
}
}
}
void print(){
for(int i=tot-1;i>=0;i--){
for(int j=k[i].size()-1;j>=0;j--)
printf("%d ",k[i][j]);
putchar('\n');
}
}
int main(){
freopen("work.in","r",stdin);
freopen("work.out","w",stdout);
t=in();n=in();
for(int i=1;i<=n;i++) a[i]=in();
memset(f,127,sizeof(f));f[1][a[1]]=0;
for(int i=2;i<=n;i++){
for(int j=t;j>0;j--){
f[i][a[i]]=min(f[i][a[i]],f[i-1][j]+(t-j)*(t-j)); //添加机器
if(j-a[i]>1) f[i][j]=min(f[i][j],f[i-1][j-a[i]-1]); //继续工作
}
}
ans=0x7fffffff;
for(int i=0;i<=t;i++) if(ans>f[n][i]+(t-i)*(t-i)) ans=f[n][i]+(t-i)*(t-i),m=i;
printf("%d\n",ans);
find(n,m);
print();
return 0;
}