sdut 1961 Image Compression

http://acm.sdut.edu.cn/sdutoj/problem.php?action=showproblem&problemid=1961

      这是此次比赛中的一个题目，由于赛前，英文的题目并没有做，也没仔细看，就加上了，没想到那么简单， 本想着这道题目这么多图应该不会有人做出来，可是没想到，实际做了一下，会是这么简单， 完全就是一个递归的简单应用。这次比赛出的有些简单了。

#include <stdio.h>
#include <string.h>
#include <math.h>
char s1[100][100];
int m;
int main()
{
    void f(int x1,int x2, int y1, int y2);
    int i,j,n,s,t,tem=1;
    while(scanf("%d",&n)!=EOF)
    {
        if(n==0)
        {
            break;
        }
        scanf("%d",&m);
        for(i=0;i<=n-1;i++)
        {
            scanf("%s",s1[i]);
        }
        f(0,n-1,0,n-1);
        printf("Image %d:\n",tem++);
        for(i=0;i<=n-1;i++)
        {
            printf("%s\n",s1[i]);
        }
    }
    return 0;
}
void f(int x1,int x2,int y1,int y2)
{
    int i,j,sum1,sum2;
    double t;
    for(i=x1,sum1=0,sum2=0;i<=x2; i++)
    {
        for(j=y1; j<=y2; j++)
        {
            if(s1[i][j]=='0')
            {
                sum1++;
            }else
            {
                sum2++;
            }
        }
    }
    t=(double)(sum1)/(double)(sum2+sum1)*100-(double)m;
    if(t>0||fabs(t-0)<=1e-7)
    {
        for(i=x1; i<=x2; i++)
        {
            for(j=y1;j<=y2;j++)
            {
                s1[i][j]='0';
            }
        }
        return ;
    }
    t=(double)(sum2)/(double)(sum2+sum1)*100-(double)m;
    if(t>0||fabs(t-0)<=1e-7)
    {
        for(i=x1; i<=x2; i++)
        {
            for(j=y1;j<=y2;j++)
            {
                s1[i][j]='1';
            }
        }
        return ;
    }
    f(x1,(x1+x2)/2, (y1+y2)/2+1, y2);
    f(x1,(x1+x2)/2, y1, (y1+y2)/2);
    f((x1+x2)/2+1,x2, y1, (y1+y2)/2);
    f((x1+x2)/2+1, x2, (y1+y2)/2+1,y2);
}