2021年中国大学生程序设计竞赛女生专场-F. 地图压缩

斐林试剂捏

于 2022-09-13 15:07:21 发布

阅读量298

点赞数

分类专栏： Fehling的ACM刷题笔记文章标签：算法 c++

本文链接：https://blog.csdn.net/m0_66624446/article/details/126833945

版权

Fehling的ACM刷题笔记专栏收录该内容

5 篇文章 0 订阅

订阅专栏

这一题在我初入ACM的时候集训队的学长曾经讲过。正好拿来复习一下KMP和哈希。

题意

一个n*n的地图。q次询问给出一个范围，求该范围内最小循环的部分。

思路

显然，这个二维地图横向与纵向之间的循环情况是互不影响，最小面积只需要分别求出横向与纵向后相乘即可。

那么寻找最小循环节的方法就是利用KMP算法了。而对于行和列的预处理只需要用的哈希算法即可。这一题给我的感觉是在于对这两个字符串算法的理解程度，等这篇题解写完之后我应该会来详细写写这两个字符串算法。

代码

#include<map>
#include<cmath>
#include<stack>
#include<queue>
#include <iomanip>  
#include<cstdio>
#include<vector>
#include<string>
#include<cstring>
#include<iostream>
#include<algorithm>
#define Endl "\n"
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
const int maxn = 2000007;
const int INF = 0x3f3f3f;
const int mod = 1e9 + 7;
const double pi = acos(-1); 
inline ll read(){
	ll x = 0, f = 1; char ch; ch = getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x*f;
}
const int base=131;
ll n,q;
ll power[2050];
char mp[2050][2050];
ll hashc[2050][2050],hashr[2050][2050];
int nxt[2050];
ll a[2050];
int kmp(int len)
{
    memset(nxt,0,sizeof(nxt));
    int j=0;
    for(int i=2;i<=len;i++)
    {
        while(j!=0&&a[i]!=a[j+1])
            j=nxt[j];
        if(a[i]==a[j+1])
            j++;
        nxt[i]=j;
    }
    return len-nxt[len];
}
void solve() 
{
    n=read(),q=read();
    power[0]=1;
    for(int i=1;i<=n;i++)
        power[i]=power[i-1]*base%mod;
    for(int i=1;i<=n;i++)
        cin>>(mp[i]+1);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
            hashr[i][j]=(hashr[i][j-1]*base%mod+(mp[i][j]-'a'+1))%mod;
    for(int j=1;j<=n;j++)
        for(int i=1;i<=n;i++)
            hashc[i][j]=(hashc[i-1][j]*base%mod+(mp[i][j]-'a'+1))%mod;
    while(q--)
    {
        int x1,x2,y1,y2;
        x1=read(),y1=read(),x2=read(),y2=read();
        int cnt=0;
        for(int i=x1;i<=x2;i++)
        {
            a[++cnt]=(hashr[i][y2]-hashr[i][y1-1]*power[y2-y1+1]%mod+mod)%mod;
        }
        int x=kmp(cnt);
        cnt=0;
        memset(a,0,sizeof(a));
        for(int i=y1;i<=y2;i++)
        {
            a[++cnt]=(hashc[x2][i]-hashc[x1-1][i]*power[x2-x1+1]%mod+mod)%mod;
        }
        int y=kmp(cnt);
        cout<<x*y<<Endl;
    }
}
 
int main() {
//for(int T = read(); T; T--)
		solve();
	return 0;
}