生产产品 [VIJOS 1243]

题目地址请点击——

---

生产产品

---

Description

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

---

Input

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

---

Output

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

---

Sample Input

3 2 0 2
2 2 3
1 3 1

---

Sample Output

4

---

LIMIT

1s

---

DATA SIZE

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

---

Solution

设 $f[i][j]$ 表示第 $j$ 个步骤由 $i$ 完成，而第 $j+1$ 个步骤不由 $i$ 完成。
则

$$f[i][j]=min^{n}_{x=1}(x\neq i)min^{j-1}_{Max(y=j-l,0)}f[k][l]+\sum_{m=l+1}^jT[i][m]$$
可用单调队列和前缀和进行优化。

---

Code

#include <iostream>
#include <cstdio>
#include <cstring>

#define Min(x,y) ((x)<(y)?(x):(y))
#define Max(x,y) ((x)>(y)?(x):(y))

using namespace std;

int m,n,k,l,ans=0x3f3f3f3f;

const int oo=0x3f3f3f3f;

int t[10][100010];
int f[10][100010];
int A[10][100010];

int que[6][6][100010][3];
int top[10][10];
int tail[10][10];

int main(){
    scanf("%d%d%d%d",&m,&n,&k,&l);
    for(int i=1;i<=n;i++){
        f[i][0]=0;
        for(int j=1;j<=m;j++){
            f[i][j]=oo;
            scanf("%d",&t[i][j]);
            A[i][j]=A[i][j-1]+t[i][j];
        }
        for(int j=1;j<=n;j++)if(i!=j)top[i][j]=tail[i][j]=1;
    }
    for(int j=1;j<=m;j++){
        for(int i=1;i<=n;i++){
            for(int x=1;x<=n;x++)
                if(!que[x][i][top[x][i]][0]||x!=i){
                    while(que[x][i][top[x][i]][0]<j-l&&top[x][i]<=tail[x][i])top[x][i]++;
                    if(top[x][i]<=tail[x][i]&&!que[x][i][top[x][i]][0])
                        f[i][j]=Min(f[i][j],que[x][i][top[x][i]][1]+A[i][j]);
                    else if(top[x][i]<=tail[x][i])
                        f[i][j]=Min(f[i][j],que[x][i][top[x][i]][1]+A[i][j]+k);
                }
        }
        for(int i=1;i<=n;i++)if(f[i][j]!=oo){
            for(int x=1;x<=n;x++)if(x!=i){
                while(tail[i][x]>top[i][x]&&que[i][x][tail[i][x]][1]>=f[i][j]-A[x][j])tail[i][x]--;
                if(que[i][x][tail[i][x]][1]>=f[i][j]-A[x][j]){
                    que[i][x][tail[i][x]][1]=f[i][j]-A[x][j];
                    que[i][x][tail[i][x]][0]=j;
                }
                else{
                    que[i][x][++tail[i][x]][1]=f[i][j]-A[x][j];
                    que[i][x][tail[i][x]][0]=j;
                }
            }
        }
    }
    for(int i=1;i<=n;i++)ans=Min(ans,f[i][m]);
    printf("%d\n",ans);
    return 0;
}