参考:https://www.cnblogs.com/neighthorn/p/6705785.html

并不是黑白染色而是三色染色(还有四色的,不过是一个意思



仔细观察一下不合法情况,可以发现都是特殊边两边有格子并且两个黑格子都在的时候黄蓝不能同在,所以(黄---黑)(黑---蓝)(黑---黑)都是最大权闭合子图中的依赖边

直接按照模型建就行,把黄蓝当成黑白染色。

……但是为什么不能用struct代替map呢QAQ

  1. #include<iostream>
  2. #include<cstdio>
  3. #include<queue>
  4. #include<cstring>
  5. #include<vector>
  6. #include<map>
  7. #include<algorithm>
  8. using namespace std;
  9. const int N=100005,mv[2][3][2]={-1,0,1,0,0,1,-1,0,1,0,0,-1},inf=1e9;
  10. int c,r,n,s,t,le[N],h[N],cnt=1;
  11. vector<int>v[N];
  12. // struct que
  13. // {
  14. // int x,y;
  15. // que(int X=0,int Y=0)
  16. // {
  17. // x=X,y=Y;
  18. // }
  19. // bool operator < (const que &a) const
  20. // {
  21. // return x>a.x;
  22. // }
  23. // };
  24. // map<que,que>mp;
  25. map< pair<int,int>,pair<int,int> > mp;
  26. struct qwe
  27. {
  28. int ne,to,va;
  29. }e[N*20];
  30. int read()
  31. {
  32. int r=0,f=1;
  33. char p=getchar();
  34. while(p>'9'||p<'0')
  35. {
  36. if(p=='-')
  37. f=-1;
  38. p=getchar();
  39. }
  40. while(p>='0'&&p<='9')
  41. {
  42. r=r*10+p-48;
  43. p=getchar();
  44. }
  45. return r*f;
  46. }
  47. void add(int u,int v,int w)
  48. {
  49. cnt++;
  50. e[cnt].ne=h[u];
  51. e[cnt].to=v;
  52. e[cnt].va=w;
  53. h[u]=cnt;
  54. }
  55. void ins(int u,int v,int w)
  56. {//cout<<u<<" "<<v<<" "<<w<<endl;
  57. add(u,v,w);
  58. add(v,u,0);
  59. }
  60. int bfs()
  61. {
  62. queue<int>q;
  63. memset(le,0,sizeof(le));
  64. le[s]=1;
  65. q.push(s);
  66. while(!q.empty())
  67. {
  68. int u=q.front();
  69. q.pop();
  70. for(int i=h[u];i;i=e[i].ne)
  71. if(e[i].va>0&&!le[e[i].to])
  72. {
  73. le[e[i].to]=le[u]+1;
  74. q.push(e[i].to);
  75. }
  76. }
  77. return le[t];
  78. }
  79. int dfs(int u,int f)
  80. {
  81. if(!f||u==t)
  82. return f;
  83. int us=0;
  84. for(int i=h[u];i;i=e[i].ne)
  85. if(e[i].va>0&&le[e[i].to]==le[u]+1)
  86. {
  87. int t=dfs(e[i].to,min(e[i].va,f-us));
  88. e[i].va-=t;
  89. e[i^1].va+=t;
  90. us+=t;
  91. }
  92. if(!us)
  93. le[u]=0;
  94. return us;
  95. }
  96. int dinic()
  97. {
  98. int re=0;
  99. while(bfs())
  100. re+=dfs(s,inf);
  101. return re;
  102. }
  103. int main()
  104. {
  105. c=read(),r=read(),n=read();
  106. s=0,t=n+1;
  107. for(int i=1;i<=n;i++)
  108. {
  109. int y=read(),x=read(),w=read();
  110. // mp[que(x,y)]=que(i,w);
  111. // v[x].push_back(y);
  112. mp[make_pair(x,y)]=make_pair(i,w);
  113. v[x].push_back(y);
  114. }
  115. for(int i=1;i<=r;i++)
  116. sort(v[i].begin(),v[i].end());
  117. for(int i=1;i<=r;i++)
  118. for(int j=0;j<v[i].size();j++)
  119. {
  120. int x=i,y=v[i][j];
  121. // if((x&1)&&y%4==1)
  122. // {
  123. // if(j<v[i].size()-1&&v[i][j+1]==y+1)
  124. // ins(mp[que(x,y)].x,mp[que(x,y+1)].x,min(mp[que(x,y)].y,mp[que(x,y+1)].y));
  125. // }
  126. // else if((x&1)&&y%4==2)
  127. // {
  128. // for(int k=0;k<3;k++)
  129. // if(mp.find(que(x+mv[0][k][0],y+mv[0][k][1]))!=mp.end())
  130. // ins(mp[que(x,y)].x,mp[que(x+mv[0][k][0],y+mv[0][k][1])].x,inf);
  131. // }
  132. // else if((x&1)==0&&y%4==0)
  133. // {
  134. // if(j>0&&v[i][j-1]==y-1)
  135. // ins(mp[que(x,y)].x,mp[que(x,y-1)].x,min(mp[que(x,y)].y,mp[que(x,y-1)].y));
  136. // }
  137. // else if((x&1)==0&&y%4==3)
  138. // {
  139. // for(int k=0;k<3;k++)
  140. // if(mp.find(que(x+mv[1][k][0],y+mv[1][k][1]))!=mp.end())
  141. // ins(mp[que(x,y)].x,mp[que(x+mv[1][k][0],y+mv[1][k][1])].x,inf);
  142. // }
  143. // else if(((x+y)&1)&&(x&1))
  144. // {
  145. // for(int k=0;k<3;k++)
  146. // if(mp.find(que(x+mv[0][k][0],y+mv[0][k][1]))!=mp.end())
  147. // ins(mp[que(x,y)].x,mp[que(x+mv[0][k][0],y+mv[0][k][1])].x,inf);
  148. // ins(s,mp[que(x,y)].x,mp[que(x,y)].y);
  149. // }
  150. // else if((x&1)&&((x+y)&1)==0)
  151. // {
  152. // for(int k=0;k<3;k++)
  153. // if(mp.find(que(x+mv[1][k][0],y+mv[1][k][1]))!=mp.end())
  154. // ins(mp[que(x+mv[1][k][0],y+mv[1][k][1])].x,mp[que(x,y)].x,inf);
  155. // ins(mp[que(x,y)].x,t,mp[que(x,y)].y);
  156. // }
  157. // else if(((x+y)&1)&&(x&1)==0)
  158. // {
  159. // for(int k=0;k<3;k++)
  160. // if(mp.find(que(x+mv[1][k][0],y+mv[1][k][1]))!=mp.end())
  161. // ins(mp[que(x,y)].x,mp[que(x+mv[1][k][0],y+mv[1][k][1])].x,inf);
  162. // ins(s,mp[que(x,y)].x,mp[que(x,y)].y);
  163. // }
  164. // else
  165. // {
  166. // for(int k=0;k<3;k++)
  167. // if(mp.find(que(x+mv[0][k][0],y+mv[0][k][1]))!=mp.end())
  168. // ins(mp[que(x+mv[0][k][0],y+mv[0][k][1])].x,mp[que(x,y)].x,inf);
  169. // ins(mp[que(x,y)].x,t,mp[que(x,y)].y);
  170. // }
  171. if((x&1)&&y%4==1)
  172. {
  173. if(j<v[i].size()-1&&v[i][j+1]==y+1)
  174. ins(mp[make_pair(x,y)].first,mp[make_pair(x,y+1)].first,min(mp[make_pair(x,y)].second,mp[make_pair(x,y+1)].second));
  175. }
  176. else if((x&1)&&y%4==2)
  177. {
  178. for(int k=0; k<3; k++)
  179. if(mp.find(make_pair(x+mv[0][k][0],y+mv[0][k][1]))!=mp.end())
  180. ins(mp[make_pair(x,y)].first,mp[make_pair(x+mv[0][k][0],y+mv[0][k][1])].first,inf);
  181. }
  182. else if((x&1)==0&&y%4==0)
  183. {
  184. if(j>0&&v[i][j-1]==y-1)
  185. ins(mp[make_pair(x,y)].first,mp[make_pair(x,y-1)].first,min(mp[make_pair(x,y)].second,mp[make_pair(x,y-1)].second));
  186. }
  187. else if((x&1)==0&&y%4==3)
  188. {
  189. for(int k=0; k<3; k++)
  190. if(mp.find(make_pair(x+mv[1][k][0],y+mv[1][k][1]))!=mp.end())
  191. ins(mp[make_pair(x,y)].first,mp[make_pair(x+mv[1][k][0],y+mv[1][k][1])].first,inf);
  192. }
  193. else if(((x+y)&1)&&(x&1))
  194. {
  195. for(int k=0; k<3; k++)
  196. if(mp.find(make_pair(x+mv[0][k][0],y+mv[0][k][1]))!=mp.end())
  197. ins(mp[make_pair(x,y)].first,mp[make_pair(x+mv[0][k][0],y+mv[0][k][1])].first,inf);
  198. ins(s,mp[make_pair(x,y)].first,mp[make_pair(x,y)].second);
  199. }
  200. else if((x&1)&&((x+y)&1)==0)
  201. {
  202. for(int k=0; k<3; k++)
  203. if(mp.find(make_pair(x+mv[1][k][0],y+mv[1][k][1]))!=mp.end())
  204. ins(mp[make_pair(x+mv[1][k][0],y+mv[1][k][1])].first,mp[make_pair(x,y)].first,inf);
  205. ins(mp[make_pair(x,y)].first,t,mp[make_pair(x,y)].second);
  206. }
  207. else if(((x+y)&1)&&(x&1)==0)
  208. {
  209. for(int k=0; k<3; k++)
  210. if(mp.find(make_pair(x+mv[1][k][0],y+mv[1][k][1]))!=mp.end())
  211. ins(mp[make_pair(x,y)].first,mp[make_pair(x+mv[1][k][0],y+mv[1][k][1])].first,inf);
  212. ins(s,mp[make_pair(x,y)].first,mp[make_pair(x,y)].second);
  213. }
  214. else
  215. {
  216. for(int k=0; k<3; k++)
  217. if(mp.find(make_pair(x+mv[0][k][0],y+mv[0][k][1]))!=mp.end())
  218. ins(mp[make_pair(x+mv[0][k][0],y+mv[0][k][1])].first,mp[make_pair(x,y)].first,inf);
  219. ins(mp[make_pair(x,y)].first,t,mp[make_pair(x,y)].second);
  220. }
  221. }
  222. printf("%d\n",dinic());
  223. return 0;
  224. }

bzoj 4823: [Cqoi2017]老C的方块【最大权闭合子图】的更多相关文章

  1. bzoj 4823: [Cqoi2017]老C的方块 [最小割]

    4823: [Cqoi2017]老C的方块 题意: 鬼畜方块游戏不解释... 有些特殊边,有些四个方块组成的图形,方块有代价,删掉一些方块使得没有图形,最小化代价. 比较明显的最小割,一个图形中必须删 ...

  2. bzoj 4823 [Cqoi2017]老C的方块——网络流

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4823 一个不合法方案其实就是蓝线的两边格子一定选.剩下两部分四相邻格子里各选一个. 所以这个 ...

  3. BZOJ 4823: [Cqoi2017]老C的方块

    分析: 我觉得我的网络流白学了...QAQ... 其实数据范围本是无法用网络流跑过去的,然而出题者想让他跑过去,也就跑过去了... 看到题目其实感觉很麻烦,不知道从哪里入手,那么仔细观察所给出的有用信 ...

  4. BZOJ 4823 [Cqoi2017]老C的方块 ——网络流

    lrd的题解:http://www.cnblogs.com/liu-runda/p/6695139.html 我还是太菜了.以后遇到这种题目应该分析分析性质的. 网络流复杂度真是$O(玄学)$ #in ...

  5. BZOJ.1312.[Neerc2006]Hard Life(分数规划 最大权闭合子图)

    BZOJ 最大密度子图. 二分答案\(x\),转为求是否存在方案满足:\(边数-x*点数\geq 0\). 选一条边就必须选两个点,所以可以转成最大权闭合子图.边有\(1\)的正权,点有\(x\)的负 ...

  6. [BZOJ 1497][NOI 2006]最大获利(最大权闭合子图)

    题目:http://www.lydsy.com:808/JudgeOnline/problem.php?id=1497 分析: 这是在有向图中的问题,且边依赖于点,有向图中存在点.边之间的依赖关系可以 ...

  7. BZOJ.1497.[NOI2006]最大获利(最小割 最大权闭合子图Dinic)

    题目链接 //裸最大权闭合子图... #include<cstdio> #include<cctype> #include<algorithm> #define g ...

  8. bzoj 1565 [NOI2009]植物大战僵尸【tarjan+最大权闭合子图】

    一上来以为是裸的最大权闭合子图,上来就dinic -然后没过样例.不得不说样例还是非常良心的给了一个强连通分量,要不然就WA的生活不能自理了 然后注意到有一种特殊情况:每个植物向他保护的植物连边(包括 ...

  9. bzoj4823: [Cqoi2017]老C的方块(最小割)

    4823: [Cqoi2017]老C的方块 题目:传送门 题解: 毒瘤题ORZ.... 太菜了看出来是最小割啥边都不会建...狂%大佬强强强   黑白染色?不!是四个色一起染,四层图跑最小割... 很 ...

随机推荐

  1. 删除字符串中的"\U0000fffc"数据 textView添加图片 以及添加后属性失效的解决

    背景:在实现textView的富文本时,如果添加一张图片后,如果直接发送textView的内容时,图片会被字符串“\U0000fffc”替换. 问题:如何删除“\U0000fffc”字符串:如何替换t ...

  2. oracle 启动监听报错TNS-12547: TNS:lost contact

    https://blog.csdn.net/liqfyiyi/article/details/7534018

  3. POJ 2240 【这题貌似可以直接FLOYD 屌丝用SPFA通过枚举找正权值环 顺便学了下map】

    题意: 给了n种硬币的名称,给了m种硬币间的转换关系. 从任意兑换地点开始兑换,看是否能够通过兑换的方式增加金钱. 思路: 用SPFA不断对各个点进行松弛操作,寻找正权值的环.如果找到则输出Yes. ...

  4. react 使用 eslint 的三种代码检查方案总结,多了解点--让代码更完美....

    1.介绍 ESLint 是一个可扩展,每条规则独立,被设计为完全可配置的lint工具. 可以用来检测代码,避免低级错误 可以用来规范代码的开发风格,统一代码习惯. 2.为什么使用 ESLint ? 统 ...

  5. TeamCity - Docker创建

    // 创建Server docker run -it --name teamcity-server-instance \-v /home/tc_datadir:/data/teamcity_serve ...

  6. 进程(WINAPI),遍历并查找树状的进程信息,实现控制系统进程

    #include <TlHelp32.h> //检索系统全部进程 void showall() { PROCESSENTRY32 pe32 = {0}; pe32.dwSize = siz ...

  7. markdown 插入latex公式练习

    markdown 插入latex公式 $$公式$$表示行间公式,本来Tex中使用\(公式\)表示行内公式,但因为Markdown中\是转义字符,所以在Markdown中输入行内公式使用\\(公式\\) ...

  8. 对dispatch_async到主线程的逻辑封装成C/C++接口类型

    背景:代码里面有时候会把将要运行的内容放到主线程里面运行,但假设已经是主线程里面的代码调用dispatch_async的时候偶尔会出现crash,所以就须要推断是否已经在主线程里面了. 通常的做法类似 ...

  9. dva/dynamic

    1.安装: yarn add dva 2.引入: import dynamic from 'dva/dynamic'; * dva路由跳转 * dynamic(app, model, componen ...

  10. 快看Sample代码,速学Swift语言(2)-基础介绍 快看Sample代码,速学Swift语言(1)-语法速览

    快看Sample代码,速学Swift语言(2)-基础介绍 Swift语言是一个新的编程语言,用于iOS, macOS, watchOS, 和 tvOS的开发,不过Swift很多部分内容,我们可以从C或 ...