• UVALive


    题目链接:

    http://acm.hust.edu.cn/vjudge/problem/129730

    The Big Painting

    Time Limit: 5000MS

    题意

    给你一个模板串和待匹配串,问模板串在待匹配串中出现的次数(这里的串是二维矩阵)

    题解

    每一行做前缀和哈希。
    统计的时候先按列再按行,这样在做行的话我们可以利用滚动的形式,计算纵向的哈希值(既总的哈希值)

    代码

    #include<map>
    #include<cmath>
    #include<queue>
    #include<vector>
    #include<cstdio>
    #include<string>
    #include<cstring>
    #include<iostream>
    #include<algorithm>
    using namespace std;
    #define X first
    #define Y second
    #define mkp make_pair
    #define lson (o<<1)
    #define rson ((o<<1)|1)
    #define mid (l+(r-l)/2)
    #define sz() size()
    #define pb(v) push_back(v)
    #define all(o) (o).begin(),(o).end()
    #define clr(a,v) memset(a,v,sizeof(a))
    #define bug(a) ;//cout<<#a<<" = "<<a<<endl
    #define rep(i,a,b) for(int i=a;i<(b);i++)
    
    typedef long long LL;
    typedef vector<int> VI;
    typedef pair<int,int> PII;
    typedef vector<pair<int,int> > VPII;
    
    const int INF=0x3f3f3f3f;
    const LL INFL=0x3f3f3f3f3f3f3f3fLL;
    const double eps=1e-8;
    
    //start----------------------------------------------------------------------
    
    const int maxn=2020;
    //P进制哈希
    const int P=123;
    
    unsigned long long H[maxn][maxn],xp[maxn*maxn],Val;
    
    int n1,m1,n2,m2;
    
    char s1[maxn][maxn],s2[maxn][maxn];
    //预处理出位权
    void pre() {
    	xp[0]=1;
    	rep(i,1,maxn*maxn) xp[i]=xp[i-1]*P;
    }
    //处理前缀和的哈希值
    void get_H() {
    	Val=0;
    	rep(i,1,n1+1) rep(j,1,m1+1) {
    		Val=Val*P+(s1[i][j]-'a'+1);
    	}
    
    	clr(H,0);
    	rep(i,1,n2+1) rep(j,1,m2+1) {
    		H[i][j]=H[i][j-1]*P+(s2[i][j]-'a'+1);
    	}
    }
    //枚举,匹配
    void solve() {
    	int ans=0;
    	for(int j=1; j<=m2-m1+1; j++) {
    		unsigned long long tmp=0,tmp2=0;
    		int lj=j-1,rj=j+m1-1;
    		int i;
    		for(i=1;i<=n1;i++){
    			tmp=tmp*xp[m1]+H[i][rj]-H[i][lj]*xp[m1];
    		}
    		do{
    			if(Val==tmp) ans++;
    			tmp-=(H[i-n1][rj]-H[i-n1][lj]*xp[m1])*xp[m1*(n1-1)];
    			tmp=tmp*xp[m1]+H[i][rj]-H[i][lj]*xp[m1];
    			i++;
    		}while(i<=n2+1);
    	}
    	printf("%d
    ",ans);
    }
    
    int main() {
    	pre();
    	while(scanf("%d%d%d%d",&n1,&m1,&n2,&m2)==4&&n1) {
    		rep(i,1,n1+1) scanf("%s",s1[i]+1);
    		rep(i,1,n2+1) scanf("%s",s2[i]+1);
    		get_H();
    		solve();
    	}
    	return 0;
    }
    
    //end-----------------------------------------------------------------------
    
    /*
    4 4 10 10
    oxxo
    xoox
    xoox
    oxxo
    xxxxxxoxxo
    oxxoooxoox
    xooxxxxoox
    xooxxxoxxo
    oxxoxxxxxx
    ooooxxxxxx
    xxxoxxoxxo
    oooxooxoox
    oooxooxoox
    xxxoxxoxxo
    */
  • 相关阅读:
    bash的for循环从命令读取值
    BFS-hdu-1226-超级密码
    计算机改名导致数据库链接的诡异问题
    There is insufficient system memory to run this query 错误
    SQL Server 2008 R2的发布订阅配置实践
    MS SQL 日常维护管理常用脚本(二)
    TreeSize工具介绍
    迁移Reporting Services的方法与WMI错误
    spring-session-data-redis解决session共享的问题
    elasticSearch6源码分析(12)DiscoveryModule
  • 原文地址:https://www.cnblogs.com/fenice/p/5774901.html
Copyright © 2020-2023  润新知