Problem:
Blur

図 1 のような 10×10 のマス目の「布」があり、（1，2）のようにX座標・Y座標の値の対でマス目を示すことにします。座標値は 0 から始まる整数とします。例えば、図 1 の◎の座標は (1, 2) です。
この「布」に染料を少しずつたらして染物を作ります。たらす染料の滴には「大」「中」「小」の 3 サイズがあり、染料の滴の落ちたマス目を中心に図 1 のように周囲も色がつきます。図 1 で☆が中心、○が色のにじむ範囲です。

「布」は、最初は「まっしろ」つまり、どのマス目も色の濃さを示す値が 0 とします。染料の滴が落ちるごとに値が 1 ずつ増えていきます。「小」が（1, 2）、「中」が（3, 2）に落ちた場合、各マス目の値は図 2 の左のようになります。染料がもったいないので、図 2 の右に示すようなにじむ範囲が布の外にはみ出すような染料の落とし方はしないことになっています。また、同じ場所に複数の染料をたらすこともあります。

このような作業を何回か繰り返した結果、布にすばらしい模様が浮かび上がりましたが、残念なことに、うっかり作業経過を記録し忘れてしまいました。どうにもさっぱり思い出せませんが、かろうじて、落とした染料の滴数は覚えていました。あなたは、すばらしい染物を再現しなければなりません。すばらしい染物のデータを読みこんで、どこに、どのような染料をたらしたかを出力するプログラムを作成してください。なお、落とした染料の滴数は 12 以下とします。

入力

入力の形式は以下のとおりです：

最初の 1 行に落とした染料の滴数 n が与えられます。次の行から 10 行に各座標の色の濃さが空白区切りで与えられます。

出力

出力は n 行からなります。染料の滴の「大」を 3、「中」を 2、「小」を 1 で表し、落とした染料ごとに、そのX座標、Y座標、滴の大きさを空白で区切って1行に出力してください。

染料のたらし方はどのような順番で出力しても良いものとします。

入力例1

2
0 0 0 0 0 0 0 0 0 0
0 0 0 1 0 0 0 0 0 0
0 0 1 1 1 0 0 0 0 0
0 0 0 1 0 0 0 1 1 1
0 0 0 0 0 0 0 1 1 1
0 0 0 0 0 0 0 1 1 1
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0

出力例1

3 2 1
8 4 2

入力例2

6
0 0 1 0 0 0 0 0 0 0
0 1 1 1 0 0 0 0 0 0
1 1 1 1 1 1 1 1 0 0
0 1 1 1 1 1 1 1 0 0
0 0 1 1 1 3 1 1 0 0
0 0 1 1 3 1 1 1 0 0
0 0 1 1 1 1 1 1 1 0
0 0 1 1 1 1 1 1 1 1
0 0 0 0 0 0 1 1 1 0
0 0 0 0 0 0 0 1 0 0

出力例2

2 2 3
7 7 3
6 3 2
3 6 2
4 4 1
5 5 1

Solution in C++ (judge verdict: Time Limit Exceeded)

#include <bits/stdc++.h>

using namespace std;

#define MAX_H 10
#define MAX_W 10
#define MAX 13*12+10
#define F first
#define S second
typedef pair<int,int> pii;
typedef pair<pii,int> Ink;

int mat[MAX_H][MAX_W];
const int size[] = {5,9,13};

const int dx[3][13] = {
    {-1,0,0,1,0},
    {-1,-1,-1,0,0,0,1,1,1},
    {-2,-1,-1,-1,0,0,0,0,0,1,1,1,2}
};

const int dy[3][13] = {
    {0,-1,0,0,1},
    {-1,0,1,-1,0,1,-1,0,1},
    {0,-1,0,1,-2,-1,0,1,2,-1,0,1,0}
};

bool inField(int x,int y){
    return (0 <= x && x < MAX_W &&
            0 <= y && y < MAX_H);
}

bool check(int x,int y,int s){
    for(int i = 0 ; i < size[s] ; i++){
        int nx = x + dx[s][i];
        int ny = y + dy[s][i];
        if(!inField(nx,ny)) return 0;
        if(mat[ny][nx] == 0) return 0;
    }
    return 1;
}

void removeBlur(int x,int y,int s,bool plus){
    for(int i = 0 ; i < size[s] ; i++){
        int nx = x + dx[s][i];
        int ny = y + dy[s][i];
        mat[ny][nx] = mat[ny][nx] + (plus ? +1 : -1);
    }
}

vector<Ink> v;
bool ok[MAX];

void solve(int x,int y,int rem,int n){
    if(n == 0){
        if(rem == 0) throw 0;
        return;
    }
    if(!ok[rem]) return;
    if(rem < n*5 || rem > n*13) return;
    for(; y < MAX_H-1 ; y++){
        for(; x < MAX_W-1 ; x++){
            for(int i = 0 ; i < 3 ; i++){
                if(rem >= size[i] && check(x,y,i)){
                    removeBlur(x,y,i,0);
                    v.push_back(Ink(pii(x,y),i+1));
                    solve(x,y,rem-size[i],n-1);
                    removeBlur(x,y,i,1);
                    v.pop_back();
                }
            }
        }
        x = 1;
    }
}

void init(){
    memset(ok,0,sizeof(ok));
    ok[0] = 1;
    for(int i = 0 ; i < MAX ; i++){
        if(!ok[i]) continue;
        for(int j = 0 ; j < 3 ; j++){
            if(i+size[j] < MAX){
                ok[i+size[j]] = 1;
            }
        }
    }
}

int main(){
    int n,cnt = 0;
    init();
    cin >> n;
    for(int i = 0 ; i < 10 ; i++){
        for(int j = 0 ; j < 10 ; j++){
            cin >> mat[i][j];
            cnt += mat[i][j];
        }
    }
    try{
        solve(1,1,cnt,n);
    }catch(...){}
    for(int i = 0 ; i < (int)v.size() ; i++){
        cout << v[i].F.F << " " << v[i].F.S << " " << v[i].S << endl;
    }
    return 0;
}