Problem:
平面に長方形のシートが何枚か置かれている. これらのシートで覆われている部分の面積や周長を求めるプログラムを作成せよ.

ただし, 平面を座標平面とみなしたとき, シートの配置は次の条件 (1), (2) を満たすものとする.

(1) 各シートの長方形の 4 頂点の x, y 座標は全て 0 以上 10000 以下の整数であり, 長方形の各辺は x 軸, または y 軸に平行である.

(2) シートの枚数は高々 10000 枚以下である.

入力データ の 1 行目に長方形の個数 n と問題の種類を表す整数 r が空白で区切って書かれている. 2 行目以降の各行には, 各シートの左下頂点座標 (x1, y1) と右上頂点座標 (x2, y2) の座標値が x1, y1, x2, y2 の順で空白で区切って書いてある.

出力には r = 1 のとき 1 行目に面積を, r = 2 のとき 1 行目に面積, 2 行目に周囲の長さを出力する. いずれの場合も最後に改行を入れること.

また, テストデータの 40% では長方形の頂点の座標が 0 以上 100 以下で, そのうちの 1/2 では面積だけを求める問題である. さらに, 全体の 1/2 は面積だけを求める問題である.

入力例１

入力例２

入力例３

入力例４

5 1

5 2

2 2

3 2

0 0 3 2

0 0 3 2

0 0 8 9

2 2 8 8

1 1 2 5

1 1 2 5

0 0 9 8

3 0 4 9

0 4 6 5

0 4 6 5

5 0 7 9

3 3 5 6

3 3 5 6

5 0 7 6

5 0 7 6

　

出力例１

出力例２

出力例３

出力例４

29

29

80

45

38

36

36

入力

入力は複数のデータセットからなる．n, r がともに 0 のとき入力が終了する．データセットの数は 10 を超えない．

出力

データセットごとに, r = 1 のとき 1 行目に面積を, r = 2 のとき 1 行目に面積, 2 行目に周囲の長さを出力する.

入力例

5 1
0 0 3 2
1 1 2 5
0 4 6 5
3 3 5 6
5 0 7 6
5 2
0 0 3 2
1 1 2 5
0 4 6 5
3 3 5 6
5 0 7 6
2 2
0 0 8 9
0 0 9 8
3 2
2 2 8 8
3 0 4 9
5 0 7 9
0 0

出力例

29
29
38
80
36
45
36

上記問題文と自動審判に使われるデータは、情報オリンピック日本委員会が作成し公開している問題文と採点用テストデータです。

Solution in C (judge verdict: Time Limit Exceeded)

#include <stdio.h>

#define M 10001

unsigned char m[M+2][(M+2)/8]={{0}},a;
const unsigned char mask[]={255,254,252,248,240,224,192,128,0};
unsigned char buf[(M+2)/8]={0};

int main(void){
    int i,j,k,n,r,x,y,t,u,s,c,p,mx,my,mt,mu;
    while((scanf("%d %d",&n,&r)>=2)&&n&&r){
        for(i=1;i<M+1;i++)for(j=1;j<(M+2)/8;j++)m[i][j]=0;
        mx=my=M+9;mt=mu=-1;
        for(i=0;i<n;i++){
            scanf("%d %d %d %d",&x,&y,&t,&u);
            x++;y+=8;t++;u+=8;
            mx=mx<x?mx:x;
            my=my<y?my:y;
            mt=mt>t?mt:t;
            mu=mu>u?mu:u;
            for(j=x;j<t;j++){
                for(k=y/8;k<=u/8;k++)buf[k]=0;
                k=y/8;
                buf[k]=mask[y&7];
                for(k++;k<=u/8;k++){
                    buf[k]=255;
                }
                k--;
                buf[k]&=255^mask[u&7];
                for(k=y/8;k<=u/8;k++)m[j][k]|=buf[k];
            }
        }
        /*
        for(i=0;i<10;i++){for(j=0;j<10;j++){printf("%02x ",m[i][j]);}printf("\n");}
        printf("%d %d %d %d\n",mx,my,mt,mu);
        */
        s=0;c=0;
        for(i=mx;i<mt;i++){
            p=0;
            for(j=my/8;j<=mu/8;j++){
                if(m[i][j]){
                    a=m[i][j];
                    a = (a & 0x55) + ((a>>1) & 0x55);
                    a = (a & 0x33) + ((a>>2) & 0x33);
                    a = (a & 0x0f) + ((a>>4) & 0x0f);
                    s+=a;
                    /*
                    for(k=1;k<256;k<<=1){
                        int h=!!(m[i][j]&k);
                        if(h){
                            if(!(m[i-1][j]&k))c++;
                            if(!(m[i+1][j]&k))c++;
                        }
                        if(p!=h){
                            c++;
                            p=h;
                        }
                    }
                    */
                    for(k=1;k<256;k<<=1){
                        int h=m[i][j]&k;
                        if(!!h!=p){
                            p=!!h;
                            c++;
                        }
                    }
                    a=m[i][j]&(255^m[i-1][j]);
                    a = (a & 0x55) + ((a>>1) & 0x55);
                    a = (a & 0x33) + ((a>>2) & 0x33);
                    a = (a & 0x0f) + ((a>>4) & 0x0f);
                    c+=a;
                    a=m[i][j]&(255^m[i+1][j]);
                    a = (a & 0x55) + ((a>>1) & 0x55);
                    a = (a & 0x33) + ((a>>2) & 0x33);
                    a = (a & 0x0f) + ((a>>4) & 0x0f);
                    c+=a;
                }else{
                    if(p){
                        c++;
                        p=0;
                    }
                }
            }
            if(p)c++;
        }
        printf("%d\n",s);
        if(r==2)printf("%d\n",c);
    }
    return 0;
}