众所周知，LZF很喜欢打一个叫Jubeat的游戏。这是个音乐游戏，游戏界面是4×4的方阵，会根据音乐节奏要求玩家按下一些指定方块（以下称combo）。LZF觉得这太简单了，于是自己仿了个游戏叫Jubeeeeeat，唯一不同之处就是界面大小，Jubeeeeeat的界面为n×n的方阵。
在某一刻，界面同时出现了若干个combo。LZF终于觉得有些困难了，但毕竟LZF不是普通人，他有很多只手。LZF的手分为m只“肉质手”和q只“意念手”。顾名思义，“肉质手”是实际存在的手，每只肉质手都有5根手指，每根手指能按一个combo，但每只手的速度都不同，受限于此，LZF的每只肉质手的控制范围是一个固定大小的正方形。“意念手”即虚无之手，每只手只有1根手指，但控制范围为全局。
现在LZF想知道，他最多能按下多少个combo。
- 输入
- 输入文件名为 jubeeeeeat.in。
第1行输入三个正整数n，m，q。
接下来是一个n×n的01矩阵，描述combo的位置，1为combo。
最后m行每行三个正整数xi，yi，ai，分别表示第i只肉质手掌控区域左上方块的行、列和边长。（行、列从1数起） - 输出
- 输出文件名为 jubeeeeeat.out。
输出一个正整数，表示最多能按下的combo数。 - 样例输入
3 1 3 1 0 1 1 1 1 1 0 1 1 1 2
- 样例输出
6
- 提示
- 【数据说明】
对于20%的数据，n=5，m=2，q=2；
对于50%的数据，1≤n≤20，1≤m, q≤50；
对于100%的数据，1≤n≤40，1≤m, q≤300，1≤xi, yi≤n，1≤xi+ai-1, yi+ai-1≤n。 - 网络流
- 源点向每一只肉质手连一条流量为5的边，然后每一只肉质手向他所能管到的格子连一条流量为1的边， 每一个格子向源点连一条流量为1的边，然后跑最大流
#include<queue> #include<cstdio> #include<cstring> #include<iostream> #include<algorithm> #define N 200100 #define inf 9999999 using namespace std; queue<int>q; ,a[][]; int to[N],lev[N],cap[N],cnt[N],head[N],nextt[N]; int read() { ,f=; char ch=getchar(); ;ch=getchar();} +ch-',ch=getchar(); return x*f; } int add(int x,int y,int z) { tot++;to[tot]=y,cap[tot]=z,nextt[tot]=head[x];head[x]=tot; tot++;to[tot]=x;cap[tot]=,nextt[tot]=head[y];head[y]=tot; } inline bool bfs() { while(!q.empty()) q.pop(); ;i<=n*n+m+;i++) { lev[i]=-; cnt[i]=head[i]; } q.push(s),lev[s]=; while(!q.empty()) { int x=q.front();q.pop(); for(int i=head[x];i;i=nextt[i]) { int t=to[i]; &&lev[t]==-) { lev[t]=lev[x]+; q.push(t); if(t==e) return true; } } } return false; } int dinic(int x,int flow) { if(x==e) return flow; ,delta; for(int &i=cnt[x];i;i=nextt[i]) { int t=to[i]; &&lev[t]==lev[x]+) { delta=dinic(t,min(cap[i],flow-rest)); if(delta) { rest+=delta; cap[i]-=delta; cap[i^]+=delta; if(rest==flow) break; } } } ; return rest; } int main() { n=read(),m=read(),p=read(); ;i<=n;i++) ;j<=n;j++) { a[i][j]=read(); if(a[i][j]) maxn++; } e=n*n+m+; ;i<=m;i++) add(s,i,); ;i<=m;i++) { x=read(),y=read(),z=read(); for(int j=x;j<x+z;j++) for(int k=y;k<y+z;k++) if(a[j][k]) add(i,m+(j-)*n+k,); } ;i<=n;i++) ;j<=n;j++) if(a[i][j]) add(m+(i-)*n+j,e,); while(bfs()) ans+=dinic(s,inf); ans=min(ans+p,maxn); printf("%d",ans); ; }