Codeforces677 E. Vanya and Balloons（对数转换+dp预处理）

题意：

解法：

显然选的数越多越好,但是不能选到0.

对于十字:
1.预处理行的前缀积,列的前缀积.
2.预处理每个位置向四个方向能扩展的最大距离.
最后枚举格子,取最大值即可.

对于叉:
和十字差不多,就是修改为预处理斜线.

但是乘积非常大,没办法存下来,
因为我们只需要比较大小,因此可以利用对数性质log(x*y)=log(x)+log(y),
将乘法转化为加法,这样的话用double就可以直接存了,大小也可以直接比较.

除了对数转化的部分,其他部分都很常规,直接写就行了.

code：

#include<bits/stdc++.h>
using namespace std;
const int maxm=1e3+5;
const int mod=1e9+7;
double ls[maxm][maxm];
double us[maxm][maxm];
double lus[maxm][maxm];
double rus[maxm][maxm];
int lp[maxm][maxm];
int rp[maxm][maxm];
int up[maxm][maxm];
int dp[maxm][maxm];
int lup[maxm][maxm];
int rup[maxm][maxm];
int ldp[maxm][maxm];
int rdp[maxm][maxm];
char a[maxm][maxm];
int n;
void solve(){
    cin>>n;
    for(int i=1;i<=n;i++){
        cin>>(a[i]+1);
        for(int j=1;j<=n;j++){
            a[i][j]-='0';
        }
    }
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            ls[i][j]=ls[i][j-1];
            us[i][j]=us[i-1][j];
            lus[i][j]=lus[i-1][j-1];
            rus[i][j]=rus[i-1][j+1];
            if(a[i][j]!=0){
                ls[i][j]+=log(a[i][j]);
                us[i][j]+=log(a[i][j]);
                lus[i][j]+=log(a[i][j]);
                rus[i][j]+=log(a[i][j]);
            }
        }
    }
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            if(!a[i][j])continue;
            lp[i][j]=lp[i][j-1]+1;
            up[i][j]=up[i-1][j]+1;
            lup[i][j]=lup[i-1][j-1]+1;
        }
        for(int j=n;j>=1;j--){
            if(!a[i][j])continue;
            rp[i][j]=rp[i][j+1]+1;
            rup[i][j]=rup[i-1][j+1]+1;
        }
    }
    for(int i=n;i>=1;i--){
        for(int j=1;j<=n;j++){
            if(!a[i][j])continue;
            dp[i][j]=dp[i+1][j]+1;
            ldp[i][j]=ldp[i+1][j-1]+1;
        }
        for(int j=n;j>=1;j--){
            if(!a[i][j])continue;
            rdp[i][j]=rdp[i+1][j+1]+1;
        }
    }
    double ma=-1;
    int f=-1,x=0,y=0,len1=0;
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            if(!a[i][j])continue;
            int len=min({lp[i][j],rp[i][j],up[i][j],dp[i][j]});
            double sum=0;
            sum+=ls[i][j+len-1]-ls[i][j-len];
            sum+=us[i+len-1][j]-us[i-len][j];
            sum-=log(a[i][j]);
            if(sum>ma){
                f=0;ma=sum;len1=len;x=i,y=j;
            }
            len=min({lup[i][j],rup[i][j],ldp[i][j],rdp[i][j]});
            sum=0;
            sum+=lus[i+len-1][j+len-1]-lus[i-len][j-len];
            sum+=rus[i+len-1][j-len+1]-rus[i-len][j+len];
            sum-=log(a[i][j]);
            if(sum>ma){
                f=1;ma=sum;len1=len;x=i,y=j;
            }
        }
    }
    int ans=1;
    if(f==0){//十字
        for(int i=x-len1+1;i<=x+len1-1;i++){
            ans=1LL*ans*a[i][y]%mod;
        }
        for(int j=y-len1+1;j<=y+len1-1;j++){
            if(j==y)continue;
            ans=1LL*ans*a[x][j]%mod;
        }
    }else if(f==1){//叉
        for(int i=x-len1+1,j=y-len1+1;i<=x+len1-1;i++,j++){
            ans=1LL*ans*a[i][j]%mod;
        }
        for(int i=x-len1+1,j=y+len1-1;i<=x+len1-1;i++,j--){
            if(i==x)continue;
            ans=1LL*ans*a[i][j]%mod;
        }
    }else if(f==-1){//矩阵中只有0的情况
        ans=0;
    }
    cout<<ans<<endl;
}
signed main(){
    ios::sync_with_stdio(0);cin.tie(0);
    solve();
    return 0;
}

---