Openjudge jubeeeeeat

jubeeeeeat

题目链接

总时间限制: 

1000ms

 

内存限制: 

256000kB

描述

众所周知，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，肉质手与他所能按的键中有combe的键连一条容量为1的边
    由于按每个combe都只会有1的贡献，所以所有的combe向汇点连容量为1的边
    原点向每只肉质手连容量为INF的边 
*/
#include<iostream>
#include<cstdio>
#include<queue>
#include<algorithm>
#define maxn 2501
using namespace std;
int n,m,q,ans,num=1,head[maxn],cur[maxn],lev[maxn],map[41][41],s,t,sum;
struct node{
    int to,pre,cap;
}e[500001];
void Insert(int from,int to,int v){
    e[++num].to=to;
    e[num].cap=v;
    e[num].pre=head[from];
    head[from]=num;
}
bool bfs(){
    queue<int>q;
    q.push(s);
    for(int i=0;i<=t;i++)cur[i]=head[i],lev[i]=-1;
    lev[s]=0;
    while(!q.empty()){
        int now=q.front();q.pop();
        for(int i=head[now];i;i=e[i].pre){
            int to=e[i].to;
            if(e[i].cap>0&&lev[to]==-1){
                lev[to]=lev[now]+1;
                q.push(to);
                if(to==t)return 1;
            }
        }
    }
    return 0;
}
int dinic(int now,int flow){
    if(now==t)return flow;
    int rest=0,delta;
    for(int &i=cur[now];i;i=e[i].pre){
        int to=e[i].to;
        if(e[i].cap>0&&lev[to]==lev[now]+1){
            delta=dinic(to,min(flow-rest,e[i].cap));
            if(delta){
                e[i].cap-=delta;e[i^1].cap+=delta;
                rest+=delta;
                if(rest==flow)break;
            }
        }
    }
    if(rest!=flow)lev[now]=-1;
    return rest;
}
int main(){
    scanf("%d%d%d",&n,&m,&q);
    for(int i=1;i<=n;i++)for(int j=1;j<=n;j++){
        scanf("%d",&map[i][j]);
        if(map[i][j])sum++;
    }
    int x,y,z;
    for(int i=1;i<=m;i++){
        Insert(s,i,9999999);
        Insert(i,s,0);
    }
    for(int i=1;i<=m;i++){
        Insert(i,i+m,5);
        Insert(i+m,i,0);
    }
    s=0,t=n*n+m*2+1;
    for(int i=1;i<=m;i++){
        scanf("%d%d%d",&x,&y,&z);
        for(int j=x;j<x+z;j++)
            for(int k=y;k<y+z;k++)
                if(map[j][k]){
                    Insert(i+m,(j-1)*n+k+m*2,1);
                    Insert((j-1)*n+k+m*2,i+m,0);
                }
    }
    for(int i=1;i<=n*n;i++){
        Insert(2*m+i,t,1);
        Insert(t,2*m+i,0);
    }
    while(bfs())ans+=dinic(s,9999999);
    int w=min(q,sum-ans);
    ans+=w;
    printf("%d",ans);
}