1249: 竞争性酶抑制剂和同工酶

Time Limit: 1 Sec  Memory Limit: 128 MB
Submit: 109  Solved: 49

Description

人体内很多化学反应都需要酶来催化。酶的功能可以简单理解为:将一种物质(底物)转化为另一种物质(目标产物)。
竞争性酶抑制剂会与底物竞争酶上的结合位点,当抑制剂达到一定剂量时,底物便竞争不过抑制剂,难以与酶结合,从而使反应无法进行。
结构不同但催化相同化学反应的酶称为一组同工酶。通常一种抑制剂只能抑制一种酶。当一种酶被它的抑制剂所抑制时,可以通过同工酶的催化使反应得以继续进行。如果一组同工酶全部被抑制,反应自然就无法再进行。但人体内的反应是千变万化的,一条反应途径被阻断,还可以通过其他反应途径,使底物经过多步转化,最终转化为目标产物。
现在已知各种物质之间的转化关系及抑制每种酶所需的抑制剂剂量,那么最少需要多少剂量的抑制剂,才能彻底阻断某种两种物质之间的转化呢?

Input

多组测试数据。对于每一组测试数据:
第一行两个整数:N、M,分别表示物质的种数、酶的种数(2<=N<=150)(0<=M<=5000)。N种物质分别编号为1到N。
接下来M行,每行描述一种酶。一行有三个整数A、B、C,表示这种酶可将A物质转化为B物质;若要抑制这种酶,需要相应的抑制剂C克(0<=C<=100000)。这M种酶中,有不少是同工酶,同工酶不超过250组。
最后一行,两个整数S、D,表示要彻底阻止S物质转化为D物质。

Output

每组测试数据输出一行。所需抑制剂的最小总量。

Sample Input

  1. 5 6
  2. 2 1 2
  3. 3 5 1
  4. 2 3 7
  5. 1 5 3
  6. 3 4 4
  7. 4 5 5
  8. 2 5
  9. 3 4
  10. 1 3 7
  11. 2 3 5
  12. 1 3 6
  13. 1 2 3
  14. 1 3
  15. 3 2
  16. 1 2 2
  17. 1 3 4
  18. 2 3
  19. 150 0
  20. 1 150

Sample Output

  1. 7
  2. 16
  3. 0
  4. 0

HINT

 

Source

CSU Monthly 2012 Apr.

解题:很明显最小割

  1. #include <iostream>
  2. #include <cstdio>
  3. #include <cstring>
  4. #include <cmath>
  5. #include <algorithm>
  6. #include <climits>
  7. #include <vector>
  8. #include <queue>
  9. #include <cstdlib>
  10. #include <string>
  11. #include <set>
  12. #include <stack>
  13. #define LL long long
  14. #define pii pair<int,int>
  15. #define INF 0x3f3f3f3f
  16. using namespace std;
  17. const int maxn = ;
  18. struct arc{
  19. int to,flow,next;
  20. arc(int x = ,int y = ,int z = -){
  21. to = x;
  22. flow = y;
  23. next = z;
  24. }
  25. };
  26. arc e[maxn<<];
  27. int head[maxn],d[maxn],cur[maxn];
  28. int tot,S,T,n,m;
  29. void add(int u,int v,int flow){
  30. e[tot] = arc(v,flow,head[u]);
  31. head[u] = tot++;
  32. e[tot] = arc(u,,head[v]);
  33. head[v] = tot++;
  34. }
  35. bool bfs(){
  36. queue<int>q;
  37. memset(d,-,sizeof(d));
  38. d[T] = ;
  39. q.push(T);
  40. while(!q.empty()){
  41. int u = q.front();
  42. q.pop();
  43. for(int i = head[u]; ~i; i = e[i].next){
  44. if(e[i^].flow && d[e[i].to] == -){
  45. d[e[i].to] = d[u] + ;
  46. q.push(e[i].to);
  47. }
  48. }
  49. }
  50. return d[S] > -;
  51. }
  52. int dfs(int u,int low){
  53. if(u == T) return low;
  54. int tmp = ,a;
  55. for(int &i = cur[u]; ~i; i = e[i].next){
  56. if(e[i].flow && d[e[i].to]+==d[u]&&(a=dfs(e[i].to,min(e[i].flow,low)))){
  57. e[i].flow -= a;
  58. e[i^].flow += a;
  59. low -= a;
  60. tmp += a;
  61. if(!tmp) break;
  62. }
  63. }
  64. if(!tmp) d[u] = -;
  65. return tmp;
  66. }
  67. int dinic(){
  68. int ans = ;
  69. while(bfs()){
  70. memcpy(cur,head,sizeof(head));
  71. ans += dfs(S,INF);
  72. }
  73. return ans;
  74. }
  75. int main() {
  76. int u,v,w;
  77. while(~scanf("%d %d",&n,&m)){
  78. memset(head,-,sizeof(head));
  79. for(int i = tot = ; i < m; ++i){
  80. scanf("%d %d %d",&u,&v,&w);
  81. add(u,v,w);
  82. }
  83. scanf("%d %d",&S,&T);
  84. printf("%d\n",dinic());
  85. }
  86. return ;
  87. }

CSU 1249 竞争性酶抑制剂和同工酶的更多相关文章

  1. csu 1812: 三角形和矩形 凸包

    传送门:csu 1812: 三角形和矩形 思路:首先,求出三角形的在矩形区域的顶点,矩形在三角形区域的顶点.然后求出所有的交点.这些点构成一个凸包,求凸包面积就OK了. /************** ...

  2. CSU 1503 点到圆弧的距离(2014湖南省程序设计竞赛A题)

    题目链接:http://acm.csu.edu.cn/OnlineJudge/problem.php?id=1503 解题报告:分两种情况就可以了,第一种是那个点跟圆心的连线在那段扇形的圆弧范围内,这 ...

  3. CSU 1120 病毒(DP)

    题目链接:http://acm.csu.edu.cn/OnlineJudge/problem.php?id=1120 解题报告:dp,用一个串去更新另一个串,递推方程是: if(b[i] > a ...

  4. CSU 1116 Kingdoms(枚举最小生成树)

    题目链接:http://acm.csu.edu.cn/OnlineJudge/problem.php?id=1116 解题报告:一个国家有n个城市,有m条路可以修,修每条路要一定的金币,现在这个国家只 ...

  5. CSU 1113 Updating a Dictionary(map容器应用)

    题目链接:http://acm.csu.edu.cn/OnlineJudge/problem.php?id=1113 解题报告:输入两个字符串,第一个是原来的字典,第二个是新字典,字典中的元素的格式为 ...

  6. CSU 1333 Funny Car Racing (最短路)

    题目链接: http://acm.csu.edu.cn/OnlineJudge/problem.php?id=1333 解题报告:一个图里面有n个点和m条单向边,注意是单向边,然后每条路开a秒关闭b秒 ...

  7. CSU 1337 搞笑版费马大定理(2013湖南省程序设计竞赛J题)

    题目链接:http://acm.csu.edu.cn/OnlineJudge/problem.php?id=1337 解题报告:虽然x和y的范围都是10^8,但是如果a 是大于1000的话,那么a^3 ...

  8. CSU 1328 近似回文词(2013湖南省程序设计竞赛A题)

    题目链接:http://acm.csu.edu.cn/OnlineJudge/problem.php?id=1328 解题报告:中文题题意就不说了.还好数据不大,只有1000,枚举回文串的中心位置,然 ...

  9. hdu 1249 三角形

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1249 part=3*s*(s-1)+2 #include<stdio.h> #includ ...

随机推荐

  1. Struts1、Struts2、Hibernate、Spring框架工作原理介绍

    Struts1工作原理 Struts1工作原理图 1.初始化:struts框架的总控制器ActionServlet是一个Servlet,它在web.xml中配置成自动启动的Servlet,在启动时总控 ...

  2. Mac 如何修改Mac系统的默认截图路径

    step 1 :打在桌面或者其他任意位置创建一个文件夹:截图图库.我创建的路径是:/Users/yilin/Documents/截图图库(仅供参考) step 2:打开终端,输入以下命令:defaul ...

  3. 运维派 企业面试题4&5 创建10个 用户 ; ping探测主机是否在线

    Linux运维必会的实战编程笔试题(19题) 企业面试题4: 批量创建10个系统帐号oldboy01-oldboy10并设置密码(密码为随机8位字符串). #!/bin/bash # ;i<=; ...

  4. 紫书 例题8-11 UVa 10954 (优先队列)

    解法和合并果子是一样的, 每次取最小的两个, 更新答案, 加入队列 #include<cstdio> #include<queue> #define REP(i, a, b) ...

  5. John Morgan:黎曼几何、曲率、Ricci流以及在三维流形上的应用二讲

    本文是笔者在线看Lektorium上John Morgan在圣彼得堡国立大学欧拉研究所的讲座做的笔记.第一讲以如下内容组成 1. 黎曼曲面上的联络 黎曼流形$(M^n,g)$中,$M$为$n$维流形, ...

  6. 制作PC端的安装程序

    一个多月不写博客了,不造大家有没有想我,(别自恋了,寥寥无几的粉丝,谁会想你),呜呜~~~ 好了,废话少叙,借用郭德纲老板的话,天儿不早了,干点正事儿吧! 一.序 Unity开发者都知道,打包出来的e ...

  7. 【codeforces 67A】Partial Teacher

    [题目链接]:http://codeforces.com/problemset/problem/67/A [题意] 给一个长度为n-1的字符串; 每个字符串是'L','R','='这3种字符中的一个; ...

  8. 2015 Multi-University Training Contest 1 Assignment

    Assignment Time Limit: 4000/2000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others)Total ...

  9. 如何将visual studio 2010编辑模式改为插入???

    按一下键盘上的 insert button 反之亦然

  10. 使用IR2101半桥驱动电机的案例

    作为一个电机驱动开发方面的菜鸟,近日研究了一下通过MOS管对整流后的电源斩波用以驱动直流电机进行调速的方案. 在驱动的过程中,遇到了很多问题,当然也有很多的收获. 写下来以供自己将来查阅,也为其他菜鸟 ...