zoukankan      html  css  js  c++  java
  • DLUTOJ 1033 Matrix

    传送门

    Time Limit: 2 Sec  Memory Limit: 128 MB

    Description

    We often use the matrix to analyze reality models. There are lots of algorithm about matrix in Linear Algebra like matrix multiplication, matrix determinant and matrix inversion, etc.

    Recently, I should use matrix to do structural mechanics analysis. The element in the matrix indicating the mechanical properties of each unit in the structure. Stable sub-structure means a part with same mechanical properties. I want to find the largest stable sub-struture as it has good engineering applications. Reflected in the matrix, the problem above equals to find the largest sub-matrix whose members have the same value.

    To accomplish the task perfectly, I wish you can help me to design a good algorithm to solve this problem.

    Input

    There are multiple test cases.

    The first line contains two integers N and M, indicating the size of this N * M matrix A.

    The next N line, each line containing M integers. The j-th integer in the i-th line means the element A(i, j).

    1 <= N, M <= 800
    1 <= A(i, j) <= 1000

    Output

    For each test, output the size of the largest sub-matrix satisfied the requests.

    Sample Input

    3 3
    1 1 1
    1 2 1
    1 1 1
     
     
    2 2
    1 2
    3 4
     
    4 4
    1 1 1 2
    1 3 3 2
    5 3 3 2
    6 6 6 7

    Sample Output

    3 1 4

    HINT

    Source

    2013大连市赛


    Solution

    单调栈


    Implementation

    #include <cstdio>
    #include <stack>
    using namespace std;
    typedef long long LL;
    
    const int N(800+5);
    
    int h[N], L[N], R[N], a[N][N];
    
    stack<int> st;
    //[L[i], R[i])
    int mono_stack(int l, int r){
        for(; st.size(); st.pop());
        for(int i=l; i<r; i++){
            for(; !st.empty()&&h[st.top()]>=h[i]; st.pop());
            if(st.empty()) L[i]=l;
            else L[i]=st.top()+1;
            st.push(i);
        }
        for(; st.size(); st.pop());
        for(int i=r-1; i>=l; i--){
            for(; !st.empty() && h[st.top()]>=h[i]; st.pop());
            if(st.empty()) R[i]=r;
            else R[i]=st.top();
            st.push(i);
        }
        int res=0;
        for(int i=l; i<r; i++)
            res=max(res, h[i]*(R[i]-L[i]));
        return res;
    }
    
    void solve(int n, int m){
        int res=0;
        for(int i=0; i<n; i++){
            if(i==0) for(int j=0; j<m; j++) h[j]=1;
            else for(int j=0; j<m; j++)
                if(a[i][j]==a[i-1][j]) h[j]++; 
                else h[j]=1;
            //two-pointers
            for(int l=0, r; l<m; l=r){
                for(r=l+1; r<m && a[i][r]==a[i][l]; r++);
                res=max(res, mono_stack(l, r));
            }
        }
        printf("%d
    ", res);
    }
    
    int main(){
        for(int n, m; ~scanf("%d%d", &n, &m); ){
            for(int i=0; i<n; i++)
                for(int j=0; j<m; j++)
                    scanf("%d", a[i]+j);
            solve(n, m);
        }
        return 0;
    }
  • 相关阅读:
    从aptitude 换回 apt-get .
    ubuntu 通过安装TOMCAT7
    第一次ubuntu使用的故障排除 ...the public key is not available: NO_PUBKEY...
    DES加密算法应用:分组加密模式
    使用“牛顿迭代法”求解方程
    html当中如何引用js文件
    javascript当中火狐的firebug如何单步调试程序?
    给出一个javascript的Helloworld例子
    java中讲讲PrintWriter的用法,举例?
    卷积和神经网络有什么关系?
  • 原文地址:https://www.cnblogs.com/Patt/p/5067912.html
Copyright © 2011-2022 走看看