将矩阵hash。s[0]忘了弄成0,输出中间过程发现了。
hash。sort。判重。大概这样子的步骤吧。
#include#include #include #include using namespace std;#define ll unsigned long long#define rep(i,n) for(int i=1;i<=n;i++)#define clr(x,c) memset(x,c,sizeof(x))#define REP(i,s,t) for(int i=s;i<=t;i++)int read(){ int x=0;char c=getchar();bool f=true; while(!isdigit(c)) { if(c=='-') f=false;c=getchar(); } while(isdigit(c)) x=x*10+c-'0',c=getchar(); return f?x:-x;} ll n,f[2][55][55],mm[55],s[3000];void init(){ n=read(); clr(f,0);clr(s,0); rep(i,n) rep(j,n) f[0][i][j]=f[0][i][j-1]*103+read(); rep(i,n) rep(j,n) f[1][i][j]=f[1][i][j-1]*103+read(); mm[1]=1;REP(i,2,54) mm[i]=mm[i-1]*103;}ll get(int x,int y,int t,int op){ ll ans=1; REP(i,x,x+t-1) ans=(ans*127+f[op][i][y+t-1]-f[op][i][y-1]*mm[t+1]); return ans;}bool check(int x){ s[0]=0; rep(i,n-x+1) rep(j,n-x+1) s[++s[0]]=get(i,j,x,0); sort(s+1,s+s[0]+1); //rep(i,s[0]) printf("%d ",s[i]);printf("\n"); rep(i,n-x+1) rep(j,n-x+1) { ll tmp=get(i,j,x,1); //printf("=>%d\n",tmp); if(*lower_bound(s+1,s+s[0]+1,tmp)==tmp) return true; } return false;}void work(){ int l=1,r=n,ans=0,mid; while(l<=r){ mid=(l+r)>>1; if(check(mid)) ans=mid,l=mid+1; else r=mid-1; } printf("%d\n",ans);}int main(){ init();work();return 0;}
1567: [JSOI2008]Blue Mary的战役地图
Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 800 Solved: 462[][][]Description
Blue Mary最近迷上了玩Starcraft(星际争霸) 的RPG游戏。她正在设法寻找更多的战役地图以进一步提高自己的水平。 由于Blue Mary的技术已经达到了一定的高度,因此,对于用同一种打法能够通过的战役地图,她只需要玩一张,她就能了解这一类战役的打法,然后她就没有兴趣再玩儿这一类地图了。而网上流传的地图有很多都是属于同一种打法,因此Blue Mary需要你写一个程序,来帮助她判断哪些地图是属于同一类的。 具体来说,Blue Mary已经将战役地图编码为n*n的矩阵,矩阵的每个格子里面是一个32位(有符号)正整数。对于两个矩阵,他们的相似程度定义为他们的最大公共正方形矩阵的边长。两个矩阵的相似程度越大,这两张战役地图就越有可能是属于同一类的。
Input
第一行包含一个正整数n。 以下n行,每行包含n个正整数,表示第一张战役地图的代表矩阵。 再以下n行,每行包含n个正整数,表示第二张战役地图的代表矩阵。
Output
仅包含一行。这一行仅有一个正整数,表示这两个矩阵的相似程度。
Sample Input
3 1 2 3 4 5 6 7 8 9 5 6 7 8 9 1 2 3 4
Sample Output
2
HINT
样例解释:
子矩阵:5 68 9为两个地图的最大公共矩阵约定:n<=50