2245: [SDOI2011]工作安排

2245: [SDOI2011]工作安排(傳送門點我QAQ)


**

Description

**

你的公司接到了一批訂單。訂單要求你的公司提供n類產品,產品被編號爲1~n,其中第i類產品共需要Ci件。公司共有m名員工,員工被編號爲1~m員工能夠製造的產品種類有所區別。一件產品必須完整地由一名員工製造,不可以由某名員工製造一部分配件後,再轉交給另外一名員工繼續進行製造。
我們用一個由0和1組成的m*n的矩陣A來描述每名員工能夠製造哪些產品。矩陣的行和列分別被編號爲1~m和1~n,Ai,j爲1表示員工i能夠製造產品j,爲0表示員工i不能製造產品j。
如果公司分配了過多工作給一名員工,這名員工會變得不高興。我們用憤怒值來描述某名員工的心情狀態。憤怒值越高,表示這名員工心情越不爽,憤怒值越低,表示這名員工心情越愉快。員工的憤怒值與他被安排製造的產品數量存在某函數關係,鑑於員工們的承受能力不同,不同員工之間的函數關係也是有所區別的。
對於員工i,他的憤怒值與產品數量之間的函數是一個Si+1段的分段函數。當他製造第1~Ti,1件產品時,每件產品會使他的憤怒值增加Wi,1,當他製造第Ti,1+1~Ti,2件產品時,每件產品會使他的憤怒值增加Wi,2……爲描述方便,設Ti,0=0,Ti,si+1=+∞,那麼當他製造第Ti,j-1+1~Ti,j件產品時,每件產品會使他的憤怒值增加Wi,j, 1≤j≤Si+1。
你的任務是制定出一個產品的分配方案,使得訂單條件被滿足,並且所有員工的憤怒值之和最小。由於我們並不想使用Special Judge,也爲了使選手有更多的時間研究其他兩道題目,你只需要輸出最小的憤怒值之和就可以了。


**

Input

**

第一行包含兩個正整數m和n,分別表示員工數量和產品的種類數;
第二行包含n 個正整數,第i個正整數爲Ci;
以下m行每行n 個整數描述矩陣A;
下面m個部分,第i部分描述員工i的憤怒值與產品數量的函數關係。每一部分由三行組成:第一行爲一個非負整數Si,第二行包含Si個正整數,其中第j個正整數爲Ti,j,如果Si=0那麼輸入將不會留空行(即這一部分只由兩行組成)。第三行包含Si+1個正整數,其中第j個正整數爲Wi,j。


**

Output

**

僅輸出一個整數,表示最小的憤怒值之和。


**

Sample Input

**

2 3

2 2 2

1 1 0

0 0 1

1

2

1 10

1

2

1 6

**

Sample Output

**

24


**

HINT

**


**

Source

**

第一輪day2


**

Solution

**

把源點向所有貨物連一條邊 容量是需求量 cost爲0

把貨物向可以工作的人連一條邊 容量爲inf cost爲0

把所有工人向匯點連好幾條邊……cost按照題目需求分段就好……

恩 就建圖麻煩點QAQ

然而好久不打了……調了好久QAQ


**

Code

**

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
using namespace std;
const int maxn=255,inf=(1<<30);
int a[maxn][maxn],dis[maxn+maxn],from[maxn+maxn],inq[maxn+maxn],head[maxn+maxn],cnt=1,n,m,S,T;
int fw[10];
long long ans;
struct data{int to,from,next,v,c;}e[maxn*maxn];
void ins(int u,int v,int w,int c)
{
    cnt++;
    e[cnt].to=v;e[cnt].from=u;
    e[cnt].next=head[u];head[u]=cnt;
    e[cnt].v=w;e[cnt].c=c;
}
void insert(int u,int v,int w,int c)
{
    ins(u,v,w,c);ins(v,u,0,-c);
}
bool spfa()
{
    queue<int> q;
    for(int i=1;i<=T;i++)dis[i]=inf;
    inq[S]=1;dis[S]=0;q.push(S);
    while(!q.empty())
    {
        int x=q.front();q.pop();
        for(int i=head[x];i;i=e[i].next)
            if(e[i].v>0&&dis[x]+e[i].c<dis[e[i].to])
            {
                dis[e[i].to]=dis[x]+e[i].c;
                from[e[i].to]=i;
                if(!inq[e[i].to])
                {
                    inq[e[i].to]=1;q.push(e[i].to);
                }
            }
        inq[x]=0;
    }
    return dis[T]!=inf;
}
void mcf()
{
    int x=inf;
    for(int i=from[T];i;i=from[e[i].from])
        x=min(x,e[i].v);
    for(int i=from[T];i;i=from[e[i].from])
    {
        e[i].v-=x;e[i^1].v+=x;
        ans+=e[i].c*x;
    }
}



int main()
{
    scanf("%d%d",&m,&n);S=n+m+1;T=S+1;
    for(int i=1,w;i<=n;i++)scanf("%d",&w),insert(S,i,w,0);
    for(int i=1,w;i<=m;i++)
        for(int j=1;j<=n;j++)
        {   
            scanf("%d",&w);
            if(w)insert(j,i+n,inf,0);
        }
    for(int i=1;i<=m;i++)
    {
        int s,w;
        scanf("%d",&s);
        for(int j=1;j<=s;j++)scanf("%d",&fw[j]);fw[++s]=inf;
        for(int j=1;j<=s;j++)scanf("%d",&w),insert(n+i,T,fw[j]-fw[j-1],w);
    }
    while(spfa())
        mcf();
    printf("%lld",ans);
    return 0;
}
//QQ:834026555 


倒要試南牆!

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章