参考: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. ***iOS学习之Table View的简单使用和DEMO示例(共Plain普通+Grouped分组两种)

    Table View简单描述: 在iPhone和其他iOS的很多程序中都会看到Table View的出现,除了一般的表格资料展示之外,设置的属性资料往往也用到Table View,Table View ...

  2. [Poj3744]Scout YYF I (概率dp + 矩阵乘法)

    Scout YYF I Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 9552   Accepted: 2793 Descr ...

  3. Java实验--关于简单字符串回文的递归判断实验

    首先题目要求写的是递归的实验,一开始没注意要求,写了非递归的方法.浪费了一些时间,所谓吃一堑长一智.我学习到了以后看实验的时候要认真看实验中的要求,防止再看错. 以下是对此次的实验进行的分析: 1)递 ...

  4. Java后端技术书单

    写博客记录技术上使用的各种问题,这个只能算是一个打游击. 如果要把一个知识学透,最有效的方式就是系统学习,而系统学习就是看书,书本上有清晰的学习路线以及相应的技术栈. 下面是我收集的Java后端的技术 ...

  5. 【APUE】进程基础

    进程标识符:非负整数 ID为0的进程通常是是调度进程,常被称为交换进程.该进程是内核的一部分,它并不执行任何磁盘上的程序,因此也被称为系统进程 ID为1的进程是init进程,在自举过程结束时由内核调用 ...

  6. 机器学习技法总结(五)Adaptive Boosting, AdaBoost-Stump,决策树

    上一讲主要利用不同模型计算出来的g.採用aggregation来实现更好的g.假设还没有做出来g.我们能够採用bootstrap的方法来做出一系列的"diversity"的data ...

  7. VC++如何折叠代码

    工具-选项,然后在文本编辑器,C/C++中的格式设置,把大纲语句块设置为True   这样之后,还是不能像C#一样使用region折叠代码,但是可以方法和if语句都会自动显示可以折叠.   使用#pr ...

  8. SQL 视图(Views)

    SQL 视图(Views) 视图是可视化的表. 本章讲解如何创建.更新和删除视图. SQL CREATE VIEW 语句 在 SQL 中,视图是基于 SQL 语句的结果集的可视化的表. 视图包含行和列 ...

  9. hdoj 1203 I NEED A OFFER! 【另类01背包】【概率背包】

    题意:... 策略:动态规划. 由于是求至少能得到一个offer的概率,那我们能够反着求.求得不到一个offer的概率.最后用1减去就好了. 代码: #include<string.h> ...

  10. ldd

    ldd命令用于判断某个可执行的 binary 档案含有什么动态函式库 [diego@localhost ~/work/branch_dispatch_201511/rtqa_center/source ...