Vijos 生产产品

描述

在经过一段时间的经营后，dd_engi的OI商店不满足于从别的供货商那里购买产品放上货架，而要开始自己生产产品了！产品的生产需要M个步骤，每一个步骤都可以在N台机器中的任何一台完成，但生产的步骤必须严格按顺序执行。由于这N台机器的性能不同，它们完成每一个步骤的所需时间也不同。机器i完成第j个步骤的时间为T[i,j]。把半成品从一台机器上搬到另一台机器上也需要一定的时间K。同时，为了保证安全和产品的质量，每台机器最多只能连续完成产品的L个步骤。也就是说，如果有一台机器连续完成了产品的L个步骤，下一个步骤就必须换一台机器来完成。现在，dd_engi的OI商店有史以来的第一个产品就要开始生产了，那么最短需要多长时间呢？
某日Azuki.7对跃动说:这样的题目太简单,我们把题目的范围改一改
对于菜鸟跃动来说,这是个很困难的问题,他希望你能帮他解决这个问题

格式

输入格式

第一行有四个整数M, N, K, L
下面的N行，每行有M个整数。第I+1行的第J个整数为T[J,I]。

输出格式

输出只有一行，表示需要的最短时间。

样例1

样例输入1

3 2 0 2
2 2 3
1 3 1

Copy

样例输出1

4

Copy

限制

1s

提示

对于50%的数据，N<=5，L<=4，M<=10000
对于100%的数据，N<=5， L<=50000，M<=100000

#include<bits/stdc++.h>
using namespace std;
const int M=100010;
typedef long long ll;
struct node{
	ll pos,val;
}q[6][6][M];
ll f[M][6],sum[6][M],head[6][6],tail[6][6];
int main(){
	ll n,m,k,l,i,j,t,ok,ans=1e16;
	cin>>m>>n>>k>>l;
	for(i=1;i<=n;i++)
		for(j=1;j<=m;j++){
			cin>>sum[i][j];
			sum[i][j]+=sum[i][j-1];
		}
	for(i=1;i<=n;i++)
		for(j=1;j<=n;j++){
			head[i][j]=tail[i][j]=0;
			q[i][j][0]=(node){0,-k};
		}
	for(t=1;t<=m;t++){
		for(i=1;i<=n;i++){
			ok=0;
			for(j=1;j<=n;j++)
				if(i!=j){
					while(q[i][j][head[i][j]].pos<t-l)head[i][j]++;
 	                if(ok==0||f[t][i]>q[i][j][head[i][j]].val){
 	               		ok=1;
						f[t][i]=q[i][j][head[i][j]].val;
					}
				}
			f[t][i]+=k+sum[i][t];
            if(t==m)ans=min(ans,f[t][i]);
		}
		for(i=1;i<=n;i++)
			for(j=1;j<=n;j++)
				if(i!=j){
					while(head[j][i]<=tail[j][i]&&f[t][i]-sum[j][t]<=q[j][i][tail[j][i]].val)tail[j][i]--;
                	q[j][i][++tail[j][i]]=(node){t,f[t][i]-sum[j][t]};
				}
	}
	cout<<ans<<endl;
	return 0;
}