AC通道:http://www.lydsy.com/JudgeOnline/problem.php?id=2879

这题codevs上也有,不过数据不同:http://codevs.cn/problem/1935/

如果你觉得你的程序跑得很快,那么你可以交codevs,如果你觉得你的程序可能是错的,那么请你交BZOJ

BZOJ数据貌似厉害一些?不过BZOJ时限也放得开些。

这题就是上题的升级版,在修车的基础上,数据很大,所以需要通过动态加边来使其变快。

因为上题也分析到了,只有厨师走完短的路,才会走长的,那么如果没走短的路,我长的路留在这里只是增加spfa的复杂度。

然后动态加边就可以了。

这个题刚开始的思路是错的,就是我在动态加边的同时,可以动态加点,但这样需要一个数组记下这是哪个厨师,也要记下这是厨师的第几个菜。

但是发现如果我的费用流需要退流的话,我记录第几个菜的数组需要-1,但是我做不到,这会导致这个点被重复算几次,然后就比答案小了。

所以必须确定下位置,这样的话,你就不会重复的添加点了,最多重复添加相同意义的边,而因为流量限制,这条边没有实际意义。

不过这样的话,速度就会变慢...

好慢啊好慢啊...

#include<queue>
#include<cstdio>
#include<cstring>
#include<algorithm> using namespace std; const int maxn=;
const int maxm=;
const int maxp=;
const int INF=0x3f3f3f3f; struct Node{
int data,next,low,cost;
}node[maxp*maxn*]; #define now node[point].data
#define www node[point].low
#define ccc node[point].cost
#define then node[point].next int n,m,cnt,ans;
int s,t,Idex,tot;
int head[maxp];
int c[maxp],ind[maxp],ink[maxp];
int dis[maxp],pre[maxp];
int a[maxm][maxm];
bool inq[maxp];
deque<int >q; void add(int u,int v,int w,int c){
node[cnt].data=v;node[cnt].next=head[u];node[cnt].low=w;node[cnt].cost=c;head[u]=cnt++;
node[cnt].data=u;node[cnt].next=head[v];node[cnt].low=;node[cnt].cost=-c;head[v]=cnt++;
} bool SPFA(){
int x,Low=INF,ta,tb;
memset(dis,0x3f,sizeof(dis));
q.push_back(s);dis[s]=;pre[s]=pre[t]=-;
while(!q.empty()){
x=q.front();q.pop_front();inq[x]=false;
for(int point=head[x];point!=-;point=then)
if(www && dis[now]>dis[x]+ccc){
dis[now]=dis[x]+ccc;pre[now]=point^;
if(!inq[now]){
inq[now]=true;
if(q.empty() || dis[q.front()]>dis[now])
q.push_front(now);
else
q.push_back(now);
}
}
}
if(pre[t]<) return false;
for(int point=pre[t];point!=-;point=pre[now]){
if(now==s){
x=node[point^].data;
ta=(x-)/tot+;tb=x%tot+;
}
Low=min(Low,min(Low,node[point^].low));
}
for(int point=pre[t];point!=-;point=pre[now])
node[point^].low-=Low,www+=Low;
for(int i=;i<=m;i++)
add((ta-)*tot+tb,n*tot+i,,tb*a[ta][i]);
ans+=Low*dis[t];
return true;
} int main(){
#ifndef ONLINE_JUDGE
freopen("2879.in","r",stdin);
freopen("2879.out","w",stdout);
#endif
scanf("%d%d",&m,&n);
for(int i=;i<=m;i++) scanf("%d",&c[i]),tot+=c[i];
for(int i=;i<=m;i++)
for(int j=;j<=n;j++)
scanf("%d",&a[j][i]);
t=n*tot+m+;
for(int i=s;i<=t;i++) head[i]=-;
for(int i=;i<=n*tot;i++) add(s,i,,);
for(int i=;i<=m;i++)
add(n*tot+i,t,c[i],);
for(int i=;i<=n;i++)
for(int k=;k<=m;k++)
add((i-)*tot+,n*tot+k,,a[i][k]);
while(SPFA());
printf("%d",ans);
return ;
}

BZOJ2879 [Noi2012]美食节的更多相关文章

  1. [BZOJ2879] [Noi2012] 美食节 (费用流 & 动态加边)

    Description CZ市为了欢迎全国各地的同学,特地举办了一场盛大的美食节.作为一个喜欢尝鲜的美食客,小M自然不愿意错过这场盛宴.他很快就尝遍了美食节所有的美食.然而,尝鲜的欲望是难以满足的.尽 ...

  2. [BZOJ2879][NOI2012]美食节(费用流)

    设sm为所有p之和,套路地对每道菜建一个点,将每个厨师拆成sm个点,做的倒数第i道菜的代价为time*i. S向每道菜连边<0,p[i]>(前者为代价后者为流量),i菜到j厨师的第k个点连 ...

  3. [BZOJ2879][Noi2012]美食节(费用流)

    题目描述 CZ市为了欢迎全国各地的同学,特地举办了一场盛大的美食节.作为一个喜欢尝鲜的美食客,小M自然不愿意错过这场盛宴.他很快就尝遍了美食节所有的美食.然而,尝鲜的欲望是难以满足的.尽管所有的菜品都 ...

  4. BZOJ2879 [Noi2012]美食节 【费用流】

    题目 CZ市为了欢迎全国各地的同学,特地举办了一场盛大的美食节.作为一个喜欢尝鲜的美食客,小M自然不愿意错过这场盛宴.他很快就尝遍了美食节所有的美食.然而,尝鲜的欲望是难以满足的.尽管所有的菜品都很可 ...

  5. [BZOJ2879][Noi2012]美食节(最小费用最大流动态加边)

    题目:http://www.lydsy.com:808/JudgeOnline/problem.php?id=2879 分析: 和bzoj1070一样,但这题的数据范围大了很多,如果直接建图就会TLE ...

  6. 【BZOJ2879】[Noi2012]美食节 动态加边网络流

    [BZOJ2879][Noi2012]美食节 Description CZ市为了欢迎全国各地的同学,特地举办了一场盛大的美食节.作为一个喜欢尝鲜的美食客,小M自然不愿意错过这场盛宴.他很快就尝遍了美食 ...

  7. BZOJ 2879: [Noi2012]美食节

    2879: [Noi2012]美食节 Time Limit: 10 Sec  Memory Limit: 512 MBSubmit: 1834  Solved: 969[Submit][Status] ...

  8. C++之路进阶——bzoj2879(美食节)

    2879: [Noi2012]美食节 Time Limit: 10 Sec  Memory Limit: 512 MBSubmit: 1304  Solved: 702[Submit][Status] ...

  9. 【BZOJ-2879】美食节 最小费用最大流 + 动态建图

    2879: [Noi2012]美食节 Time Limit: 10 Sec  Memory Limit: 512 MBSubmit: 1366  Solved: 737[Submit][Status] ...

随机推荐

  1. Cassandra 有限分页策略

    瀑布式分页 如果你的应用只需要瀑布式的分页,那么,Cassandra可以很好的支持,不过记得要指定好排序顺序. CLUSTERING ORDER BY (add_time DESC); 常见的分页,跳 ...

  2. ASP.NET的错误处理机制之一(概念)

    对Web应用程序来说,发生不可预知的错误和异常在所难免,我们必须为Web程序提供错误处理机制.当错误发生时,我们必须做好两件事情:一是将错误信息记录日志,发邮件通知网站维护人员,方便技术人员对错误进行 ...

  3. php教程一,变量

    php是一种动态脚本语言,比较适合web开发. php支持8种变量数据类型: 四种标量类型: boolean(布尔型) integer(整型) float(浮点型,也称作 double) string ...

  4. ASP.NET MVC4学习笔记之Controller激活的扩展

    一. 为什么要进行扩展 在前面的分析中,我们知道默认的Controller激活系统只能实例化无参构造函数的Controller类型,但在某些情况一下,我们希望某些服务的实例能够自动注入到Control ...

  5. C#之匿名类型与隐式局部变量

    一.匿名类型 下面一段代码展示了如何定义并且使用匿名类型: static void Main(string[] args) { var patent1 = new { Title = "Ne ...

  6. flask学习

    安装环境: centos 6.3 python2.6 使用easy_install安装方式: [root@localhost ~]# easy_install flask 简单的hello from  ...

  7. 远程连接数据库(通过pgAdmin)

    1.编辑/var/lib/pgsql/data/pg_hba.conf,增加语句  host all all 192.168.105.225/36 trust 让数据库接受网络 192.168.105 ...

  8. 《Prism 5.0源码走读》ModuleCatalog

    概念 ModuleCatalog 是Prism中主要概念之一,主要用来保存应用程序可用的modules(模块),每个module都是用ModuleInfo来定义(包含module的名称.类型和位置). ...

  9. JS中的函数,Array对象,for-in语句,with语句,自定义对象,Prototype

    一)函数 A)JS中的函数的定义格式: function add(a,b) { var sum = a+b; document.write("两个数的和是:" + sum); // ...

  10. SQL 跟踪方法相关介绍

    oracle sql跟踪方法:1.sql_trace打开跟踪:alter session set sql_trace=true;为跟踪文件做标记:alter session set tracefile ...