acm常用技巧三 POJ 3279

一.POJ 3279

const int dx[]={-1,0,0,0,1};
const int dy[]={0,-1,0,1,0};

int M,N;
int tile[MXN][MXN];
int opt[MXM][MXN];
int flip[MXM][MXN];

int get(int x,int y){
    int c=tile[x][y];
    for(int d=0;d<5;++d){
        int x2=x+dx[d],y2=y+dy[d];
        if(0<=x2&&x2<M&&0<=y2&&y2<N)
            c+=flip[x2][y2];
    }
    return c%2;
}
int calc(){
    for(int i=1;i<M;++i)
        for(int j=0;j<N;++j)
            if(get(i-1,j))
            {
                flip[i][j]=1;
            }
    for(int j=0;j<N;++j)
        if(get(M-1,j)){
            return -1;
        }
    int res=0;
    for(int i=0;i<M;++i)
        for(int j=0;j<N;++j)
            res+=flip[i][j];
    return res;
}
void Fun(){
    int res=-1;
    for(int i=0;i<1<<N;++i){
        memset(flip,0,sizeof(flip));
        for(int j=0;j<N;++j)
            flip[0][N-j-1]=i>>j&1;
        int num=calc();
        if(num>=0&&(res<0||res>num)){
            res=num;
            memcpy(opt,flip,sizeof(flip));
        }
    }
    if(res<0)
        puts("IMPOSSIBLE");
    else{
        for(int i=0;i<M;++i)
            for(int j=0;j<N;++j)
                printf("%d%c",opt[i][j],j+1==N?'\n':' ');
    }
}

二.POJ 3684

const double g=10.0;
int N,H,R,T;
double y[MXN];
double calc(int T){
    if(T<0) return H;
    double t=sqrt(2*H/g);
    int k=(int)(T/t);
    if(k%2==0){
        double d=T-k*t;
        return H-g*d*d/2;
    }else{
        double d=k*t+t-T;
        return H-g*d*d/2;
    }
}
void Fun(){
    for(int i=0;i<N;++i)
        y[i]=calc(T-i);
    sort(y,y+N);
    for(int i=0;i<N;++i)
        printf("%.2f%c",y[i]+2*R*i/100.0,i+1==N?'\n':' ');
}

三.POJ 2785

void Fun(){
    for(int i=0;i<n;++i)
        for(int j=0;j<n;++j)
            CD[i*n+j]=C[i]+D[j];
    sort(CD,CD+n*n);
    LL res=0;
    for(int i=0;i<n;++i)
        for(int j=0;j<n;++j){
            int cd=-(A[i]+B[j]);
            res+=upper_bound(CD,CD+n*n,cd)-lower_bound(CD,CD+n*n,cd);
        }
    printf("%lld\n",res);
}