Codeforces 713D Animals and Puzzle

题意：一个n*m的01矩阵，Q个询问，每次询问一个矩形区域内，最大的全1正方形的边长是多少？

题解：dp[0][0][i][j]表示以(i, j)为右下角的正方形的最长边长。RMQ后，二分答案即可。

 #include <cstdio>
 #include <iostream>
 #include <algorithm>
 using namespace std;
 typedef long long ll;
 ;
 int x1, y1, x2, y2, n, m;
 int a[N][N];
 ][][N][N];
 int Log[N];
 int ask(int x1, int y1, int x2, int y2){
     ], k2 = Log[y2-y1+];
     <<k2)+]);
     ans = max(ans, dp[k1][k2][x2-(<<k1)+][y1]);
     ans = max(ans, dp[k1][k2][x2-(<<k1)+][y2-(<<k2)+]);
     return ans;
 }
 bool test(int l){
     //dp[0][x1][y1] dp[0][x2-x+1][y2-x+1]
     , y1+l-, x2, y2);
     return ret >= l;
 }
 void init(){
     ; i <= Log[n]; i++)
         ; x+(<<i)- <= n; x++)
             ; y <= m; y++)
                 dp[i][][x][y] = max(dp[i-][][x][y], dp[i-][][x+(<<i-)][y]);
     ; i <= Log[n]; i++)
         ; j <= Log[m]; j++)
             ; x+(<<i)- <= n; x++)
                 ; y+(<<j)- <= m; y++)
                     dp[i][j][x][y] = max(dp[i][j-][x][y], dp[i][j-][x][y+(<<j-)]);
 }
 int main(){
     ; i < N; i++) Log[i] = Log[i>>]+;
     scanf("%d%d", &n, &m);
     ; i <= n; i++)
         ; j <= m; j++){
             scanf("%d", &a[i][j]);
             ) dp[][][i][j] = ;
             else{
                 ][][i-][j], dp[][][i][j-]);
                 dp[][][i][j] = len+a[i-len][j-len];
             }
         }
     init();

     int t; scanf("%d", &t);
     while(t--){
         scanf("%d%d%d%d", &x1, &y1, &x2, &y2);
         , R = min(x2-x1+, y2-y1+);
         while(L < R){
              >> ;
             if( test(M) )
                 L = M;
             else
                 R = M-;
         }
         printf("%d\n", L);
     }
     ;
 }