Go: package main

import (
	"bufio"
	"os"
	"strings"
	"strconv"
	"fmt"
)

var scanner = bufio.NewScanner(os.Stdin)

var (
	dx = []int{0, -1, 0, 1}
	dy = []int{1, 0, -1, 0}
)

func main() {
	for scanner.Scan() {
		buf := strings.Split(scanner.Text(), " ")
		w, _ := strconv.Atoi(buf[0])
		h, _ := strconv.Atoi(buf[1])
		if w == 0 && h == 0 {
			break
		}
		scanner.Scan()
		buf_ := strings.Split(scanner.Text(), " ")
		xs, _ := strconv.Atoi(buf_[0])
		ys, _ := strconv.Atoi(buf_[1])
		scanner.Scan()
		buf__ := strings.Split(scanner.Text(), " ")
		xg, _ := strconv.Atoi(buf__[0])
		yg, _ := strconv.Atoi(buf__[1])
		xs--
		ys--
		xg--
		yg--
		scanner.Scan()
		n, _ := strconv.Atoi(scanner.Text())
		field := make([][]int, 100)
		for i := range  field {
			field[i] = make([]int, 100)
		}
		visited := make([][]bool, 100)
		for i := range visited {
			visited[i] = make([]bool, 100)
		}
		for i := 0; i < n; i++ {
			scanner.Scan()
			buf___ := strings.Split(scanner.Text(), " ")
			c, _ := strconv.Atoi(buf___[0])
			d, _ := strconv.Atoi(buf___[1])
			x, _ := strconv.Atoi(buf___[2])
			y, _ := strconv.Atoi(buf___[3])
			x--
			y--
			switch d {
			case 0:
				for j := 0; j < 4; j++ {
					for k := 0; k < 2; k++ {
						field[x+j][y+k] = c
					}
				}
			case 1:
				for j := 0; j < 2; j++ {
					for k := 0; k < 4; k++ {
						field[x+j][y+k] = c
					}
				}
			}
		}
		if field[xs][ys] != field[xg][yg] || // スタートとゴールの色が等しくないと"NG"
			field[xs][ys] == 0 { // スタートの色が無職
			fmt.Println("NG")
			continue
		}
		Q := make([]Entry, 0)
		visited[xs][ys] = true
		Q = append(Q, Entry{xs, ys})
		for len(Q) > 0 {
			row := Q[0].x
			col := Q[0].y
			Q = Q[1:]
			for i := 0; i < 4; i++ {
				if 0 <= row+dx[i] && row+dx[i] < w &&
					0 <= col+dy[i] && col+dy[i] < h && // 移動可能な範囲内かどうか
					!visited[row + dx[i]][col + dy[i]] && // 訪問済みか
					field[row+dx[i]][col+dy[i]]==field[row][col] { // 色が同じかどうか
						visited[row+dx[i]][col+dy[i]] = true
						Q = append(Q, Entry{row+dx[i],col+dy[i]})
				}
			}
		}
		if visited[xg][yg] {
			fmt.Println("OK")
		} else {
			fmt.Println("NG")
		}
	}
}

type Entry struct {
	x, y int
}

C++: #include<iostream>
#include<stack>
#include<string.h>
using namespace std;
struct ST {
	int a;
	int b;
	int c;
};
long long brock[1000][1000];
long long ds[1000][1000];
int q[4]={0,-1,0,1};
int p[4]={1,0,-1,0};
int main(){
	stack <ST> S;
	ST ma;
	ST ba;
	int w,h;//???????????§??????
	int xs,ys;//??????????????§?¨?
	int xg,yg;//??´????????§?¨?
	int n;//????????????????????°
	int c,d,x,y;//????????????????????±
	int i,j,m,check=0;
	while(1){
		memset(ds,0,sizeof(ds));
		memset(brock,0,sizeof(brock));
		cin>>w>>h;//???????????§????????\???
		if(w==0&&h==0){
			break;
		}
		cin>>xs>>ys>>xg>>yg;//??????????????¨??´???????????§?¨???\???
		cin>>n;//????????????????????°
		for(i=0;i<n;i++){
			cin>>c>>d>>x>>y;
			if(d==0){
				for(m=y;m<y+2;m++){
					for(j=x;j<x+4;j++){
						brock[m][j]=c;
					}
				}
			}
			else{
				for(m=y;m<y+4;m++){
					for(j=x;j<x+2;j++){
						brock[m][j]=c;
					}
				}
			}
		}
		/*for(i=1;i<h;i++){
			for(j=1;j<w;j++){
				cout<<brock[i][j];
			}
			cout<<endl;
		}*/
		if(brock[ys][xs]!=brock[yg][xg]){
			cout<<"NG"<<endl;
		}
		else{
			ma.a=ys;
			ma.b=xs;
			ma.c=brock[ys][xs];
			ds[ys][xs]=1;
			S.push(ma);
			check=0;
			while(!S.empty()){
				ba=S.top();
				S.pop();
				ds[ba.a][ba.b]=1;
				if(ba.a==yg&&ba.b==xg){
					check=1;
					break;
				}
				else{
					for(i=0;i<4;i++){
						if(ba.a+p[i]>0&&ba.a+p[i]<=h&&ba.b+q[i]>0&&ba.b+q[i]<=w&&ds[ba.a+p[i]][ba.b+q[i]]==0&&brock[ba.a+p[i]][ba.b+q[i]]==brock[yg][xg]){
							ma.a=ba.a+p[i];
							ma.b=ba.b+q[i];
							ma.c=brock[ba.a+p[i]][ba.b+q[i]];
							S.push(ma);
						}
					}
				}
			}
			if(check==1){
				cout<<"OK"<<endl;
			}
			else{
				cout<<"NG"<<endl;
			}
		}
	}
	return 0;
}

C++: #include <iostream>
#include <algorithm>
#include <map>
#include <set>
#include <queue>
#include <stack>
#include <numeric>
#include <bitset>
#include <cmath>

static const int MOD = 1000000007;
using ll = long long;
using u32 = unsigned;
using u64 = unsigned long long;
using namespace std;

template<class T> constexpr T INF = ::numeric_limits<T>::max()/32*15+208;

class UnionFind {
    int n;
    vector<int> uni;
    int forest_size;
public:
    explicit UnionFind(int n) : n(n), uni(static_cast<u32>(n), -1), forest_size(n) {};

    int root(int a){
        if (uni[a] < 0) return a;
        else return (uni[a] = root(uni[a]));
    }

    bool unite(int a, int b) {
        a = root(a);
        b = root(b);
        if(a == b) return false;
        if(uni[a] > uni[b]) swap(a, b);
        uni[a] += uni[b];
        uni[b] = a;
        forest_size--;
        return true;
    }
    int size(){ return forest_size; }
    int size(int i){ return -uni[root(i)]; }
    bool same(int a, int b) { return root(a) == root(b); }
};

void solve(int h, int w){
    int sy, sx, gy, gx;
    cin >> sx >> sy >> gx >> gy;
    vector<vector<int>> grid(h+2, vector<int>(w+2, 0));
    int n;
    cin >> n;
    for (int i = 0; i < n; ++i) {
        int c, d, x, y;
        cin >> c >> d >> x >> y;
        if(d){
            for (int j = 0; j < 4; ++j) {
                for (int k = 0; k < 2; ++k) {
                    grid[y+j][x+k] = c;
                }
            }
        }else {
            for (int j = 0; j < 2; ++j) {
                for (int k = 0; k < 4; ++k) {
                    grid[y+j][x+k] = c;
                }
            }
        }
    }
    UnionFind uf((h+2)*(w+2));
    for (int i = 1; i <= h; ++i) {
        for (int j = 1; j <= w; ++j) {
            if(grid[i][j] == grid[i+1][j]){
                uf.unite(i*(w+2)+j, (i+1)*(w+2)+j);
            }
            if(grid[i][j] == grid[i][j+1]){
                uf.unite(i*(w+2)+j, i*(w+2)+j+1);
            }
        }
    }
    puts(uf.same(sy*(w+2)+sx, gy*(w+2)+gx) ? "OK" : "NG");
}

int main() {
    int w, h;
    while(cin >> w >> h, w){
        solve(h, w);
    }
    return 0;
}

C++: #include<iostream>
#include<cstring>
#include<queue>
#include<map>
using namespace std;
int C,goal,n,x,y,xs,ys,xg,yg,mas[110][110]={{0}};
int dx[4]={0,0,1,1};
int dy[4]={0,1,0,1};
int d[110][110];
typedef pair<int,int>P;

int bfs();
int main(){

  while(1){
    goal=0;
    memset(mas,0,sizeof(mas));
    cin >> x >> y;//ボードの大きさ。 ※：x=横。y=縦。
    if(x==0 && y==0){
      break;
    }
    cin >> xs >> ys;//スタート地点。
    cin >> xg >> yg;//ゴール地点。
    cin >> n;//ブロックの数。
    
    //色染め。
    int color,dire,xx,yy;
    for(int t=0;t<n;t++){
      cin >> color >> dire >> xx >> yy;

      //横長。
      if(dire==0){
	for(int i=0;i<2;i++){
	  for(int j=0;j<4;j++){
	    mas[yy+i][xx+j]=color;
	  }
	}
      }
      
      //縦長。
      else{
	for(int i=0;i<4;i++){
	  for(int j=0;j<2;j++){
	    mas[yy+i][xx+j]=color;
	  }
	}
      }
      
    }
    
    //通る色　摘出。
    C=mas[ys][xs];
    
    bfs();
    
    if(goal==1){
      cout << "OK" << endl;
    }
    else{
      cout << "NG" << endl;
    }
    
  }
}

int bfs(){
  queue<P>que;
  memset(d,-1,sizeof(d));
  que.push(P(ys,xs));
  d[ys][xs]=0;
  while(!que.empty()){
    P p=que.front();
    que.pop();
    if(p.first==yg && p.second==xg){
      goal=1;
      break;
    }
    for(int i=0;i<4;i++){
      int ny=p.first+dy[i];
      int nx=p.second+dx[i];
      if(mas[ny][nx]==C && ny>=1 && ny<=y && nx>=1 && nx<=x && d[ny][nx]==-1){
	que.push(P(ny,nx));
	d[ny][nx]=0;
      }
    }
  }
  return (0);
}