bzoj2245 [SDOI2011]工作安排
2017-09-09 10:13
323 查看
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
分析:
费用流
我一开始yy的一种建图方法,如下:
结果2500ms的时候WA,估计是4、5个点的时候
结果看到一个前辈说
但是挺慢的。。卡时过
后来看到网上
用了一种拆边的方法,可以把边数大幅度减少
减少了循环边和点入队的时间
时间直接变成了开始的1/4
你的公司接到了一批订单。订单要求你的公司提供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
分析:
费用流
我一开始yy的一种建图方法,如下:
结果2500ms的时候WA,估计是4、5个点的时候
结果看到一个前辈说
答案要开ll
然后就A了但是挺慢的。。卡时过
这里写代码片 #include<cstdio> #include<cstring> #include<iostream> #define ll long long using namespace std; const int INF=0x33333333; const int N=100000; struct node{ int x,y,v,c,nxt; }; node way[N*10]; int st ,tot=-1,pre ,dis ,q ,tou,wei,cnt; int n,m,ci ,mp[300][300],S ,T[300][10],W[300][10],s,t; bool p ; void add(int u,int w,int v,int cc) { tot++; way[tot].x=u;way[tot].y=w;way[tot].v=v;way[tot].c=cc;way[tot].nxt=st[u];st[u]=tot; tot++; way[tot].x=w;way[tot].y=u;way[tot].v=0;way[tot].c=-cc;way[tot].nxt=st[w];st[w]=tot; } int spfa(int s,int t) { memset(pre,-1,sizeof(pre)); memset(p,1,sizeof(p)); memset(dis,0x33,sizeof(dis)); tou=wei=0; q[++wei]=s; dis[s]=0; p[s]=0; do { int r=q[++tou]; for (int i=st[r];i!=-1;i=way[i].nxt) if (way[i].v&&dis[r]+way[i].c<dis[way[i].y]) { dis[way[i].y]=dis[r]+way[i].c; pre[way[i].y]=i; if (p[way[i].y]) { p[way[i].y]=0; q[++wei]=way[i].y; } } p[r]=1; } while (tou<wei); return dis[t]!=INF; } void doit() { ll ans=0; while (spfa(s,t)) { int sum=INF; for (int i=t;i!=s;i=way[pre[i]].x) sum=min(sum,way[pre[i]].v); ans+=(ll)sum*dis[t]; for (int i=t;i!=s;i=way[pre[i]].x) way[pre[i]].v-=sum, way[pre[i]^1].v+=sum; } printf("%lld",ans); } void lianbian() { int i,j,k,tt=0; s=0; t=2500; for (i=1;i<=n;i++) //everyone { ++tt; add(s,tt,INF,0); //tt 当前点 int now=tt; for (j=1;j<=S[i];j++){ tt++; add(now,tt,T[i][j]-T[i][j-1],0); for (k=1;k<=m;k++) if (mp[i][k]) add(tt,k+cnt,INF,W[i][j]); } tt++; add(now,tt,INF,0); //最后一级 for (k=1;k<=m;k++) if (mp[i][k]) add(tt,k+cnt,INF,W[i][S[i]+1]); } for (i=1;i<=m;i++) add(i+cnt,t,ci[i],0); } int main() { scanf("%d%d",&n,&m); memset(st,-1,sizeof(st)); cnt=n; for (int i=1;i<=m;i++) scanf("%d",&ci[i]); for (int i=1;i<=n;i++) for (int j=1;j<=m;j++) scanf("%d",&mp[i][j]); for (int i=1;i<=n;i++) { scanf("%d",&S[i]); cnt+=S[i]+1; for (int j=1;j<=S[i];j++) scanf("%d",&T[i][j]); for (int j=1;j<=S[i]+1;j++) scanf("%d",&W[i][j]); } lianbian(); doit(); return 0; }
后来看到网上
用了一种拆边的方法,可以把边数大幅度减少
减少了循环边和点入队的时间
时间直接变成了开始的1/4
这里写代码片 #include<cstdio> #include<cstring> #include<iostream> #define ll long long using namespace std; const int INF=0x33333333; const int N=50000; struct node{ int x,y,v,c,nxt; }; node way[N*10]; int st ,tot=-1,pre ,dis ,q ,tou,wei; int n,m,ci ,mp[300][300],S ,T[300][10],W[300][10],s,t; bool p ; void add(int u,int w,int v,int cc) { tot++; way[tot].x=u;way[tot].y=w;way[tot].v=v;way[tot].c=cc;way[tot].nxt=st[u];st[u]=tot; tot++; way[tot].x=w;way[tot].y=u;way[tot].v=0;way[tot].c=-cc;way[tot].nxt=st[w];st[w]=tot; } int spfa(int s,int t) { memset(pre,-1,sizeof(pre)); memset(p,1,sizeof(p)); memset(dis,0x33,sizeof(dis)); tou=wei=0; q[++wei]=s; dis[s]=0; p[s]=0; do { int r=q[++tou]; for (int i=st[r];i!=-1;i=way[i].nxt) if (way[i].v&&dis[r]+way[i].c<dis[way[i].y]) { dis[way[i].y]=dis[r]+way[i].c; pre[way[i].y]=i; if (p[way[i].y]) { p[way[i].y]=0; q[++wei]=way[i].y; } } p[r]=1; } while (tou<wei); return dis[t]!=INF; } void doit() { ll ans=0; while (spfa(s,t)) { int sum=INF; for (int i=t;i!=s;i=way[pre[i]].x) sum=min(sum,way[pre[i]].v); ans+=(ll)sum*dis[t]; for (int i=t;i!=s;i=way[pre[i]].x) way[pre[i]].v-=sum, way[pre[i]^1].v+=sum; } printf("%lld",ans); } void lianbian() { int i,j,k,tt=0; s=0; t=n+m+1; for (i=1;i<=n;i++) //everyone { for (j=1;j<=S[i];j++) add(s,i,T[i][j]-T[i][j-1],W[i][j]); add(s,i,INF,W[i][S[i]+1]); for (j=1;j<=m;j++) if (mp[i][j]) add(i,j+n,INF,0); } for (i=1;i<=m;i++) add(i+n,t,ci[i],0); } int main() { scanf("%d%d",&n,&m); memset(st,-1,sizeof(st)); for (int i=1;i<=m;i++) scanf("%d",&ci[i]); for (int i=1;i<=n;i++) for (int j=1;j<=m;j++) scanf("%d",&mp[i][j]); for (int i=1;i<=n;i++) { scanf("%d",&S[i]); for (int j=1;j<=S[i];j++) scanf("%d",&T[i][j]); for (int j=1;j<=S[i]+1;j++) scanf("%d",&W[i][j]); } lianbian(); doit(); return 0; }
相关文章推荐
- bzoj 2245 [SDOI2011]工作安排(最小费用最大流)
- 【费用流】【bzoj 2245】: [SDOI2011]工作安排
- BZOJ 2245 SDOI 2011 工作安排 费用流
- 【BZOJ 2245】[SDOI2011]工作安排
- BZOJ 2245 SDOI 2011 工作安排 费用流
- bzoj 2245 [SDOI2011]工作安排(最小费用最大流)
- 【BZOJ】【2245】【SDOI2011】工作安排
- BZOJ 2245: [SDOI2011]工作安排
- bzoj 2245 [SDOI2011]工作安排【最小费用最大流】
- BZOJ 2245: [SDOI2011]工作安排( 费用流 )
- Bzoj 2245: [SDOI2011]工作安排(费用流)
- bzoj 2245: [SDOI2011]工作安排(费用流)
- bzoj 2245: [SDOI2011]工作安排
- 2245: [SDOI2011]工作安排|费用流
- [SDOI2011][bzoj2245] 工作分配 [费用流]
- 2245: [SDOI2011]工作安排 (最小费用最大流)
- 2245: [SDOI2011]工作安排
- 2245: [SDOI2011]工作安排
- bzoj 2245 [SDOI2011]工作安排
- 2245: [SDOI2011]工作安排 费用流