网站建设语言都有什么软件,建设网站是否等于网络营销,东莞seo网站优化方式,凡网站建设正题 题目链接:https://www.ybtoj.com.cn/problem/526 题目大意 一个nmn\times mnm的网格上有字母#xff0c;你每次可以沿平行坐标轴对折网格#xff0c;要求对折的对应位置字母相同。 询问有多少个可能对折出来的子矩阵。 1≤nm≤1061\leq n\times m\leq 10^61≤nm≤106…正题 题目链接:https://www.ybtoj.com.cn/problem/526 题目大意 一个n×mn\times mn×m的网格上有字母你每次可以沿平行坐标轴对折网格要求对折的对应位置字母相同。 询问有多少个可能对折出来的子矩阵。 1≤n×m≤1061\leq n\times m\leq 10^61≤n×m≤106 解题思路 首先行和列是独立的行的对折不会和列的对折有任何关联所以可以分开考虑行和列可以对折出的区间。 然后设每一行分开对每个轴求出一个最大对折距离这个用二分hashhashhash或者马拉车就可以求出来了然后同位置的所有行取最小值就好了。 之后对于每个轴的位置就有一个可以转移过来的区间而且左右的对折如果过头了不会影响答案可以自己画个图因为回文串的性质那么两边一定可以先对折出一个更小不会冲突的区间 维护一个前缀和就好了考场上犯病写了个树状数组 时间复杂度O(nlog⁡n)O(n\log n)O(nlogn) code #includecstdio #includecstring #includealgorithm #define ll long long #define ull unsigned long long #define lowbit(x) (x-x) using namespace std; const ll N1e610; const ull g131; ll n,m,t[N],ac[N],cr[N],dp[N],lim; ull h[N],f[N],pw[N]; char c[N],*s[N]; void Change(ll x,ll val){while(xlim){t[x]val;xlowbit(x);}return; } ll Ask(ll x){ll ans0;while(x){anst[x];x-lowbit(x);}return ans; } ll Query(ll l,ll r) {return Ask(r)-Ask(l-1);} ull geth(ll l,ll r) {return h[r]-h[l-1]*pw[r-l1];} ull getf(ll l,ll r) {return f[l]-f[r1]*pw[r-l1];} signed main() {freopen(paper.in,r,stdin);freopen(paper.out,w,stdout);scanf(%lld%lld,n,m);pw[0]1;for(ll i1;imax(n,m);i)pw[i]pw[i-1]*g;memset(ac,0x3f,sizeof(ac));memset(cr,0x3f,sizeof(cr));s[1]c-1;for(ll p1;pn;p){scanf(%s,s[p]1);for(ll i1;im;i)h[i]h[i-1]*gs[p][i]-a;for(ll im;i1;i--)f[i]f[i1]*gs[p][i]-a;for(ll i2;im;i){ll l0,rmin(i-2,m-i);while(lr){ll mid(lr)1;if(geth(i-mid-1,i-1)getf(i,imid))lmid1;else rmid-1;}ac[i]min(ac[i],r);}s[p1]s[p]m;}f[n1]0;for(ll p1;pm;p){for(ll i1;in;i)h[i]h[i-1]*gs[i][p]-a;for(ll in;i1;i--)f[i]f[i1]*gs[i][p]-a;for(ll i2;in;i){ll l0,rmin(i-2,n-i);while(lr){ll mid(lr)1;if(geth(i-mid-1,i-1)getf(i,imid))lmid1;else rmid-1;}cr[i]min(cr[i],r);}}limm;Change(1,1);dp[1]1;for(ll i2;im;i){bool tmp(Query(i-ac[i]-1,i-1)!0);dp[i]dp[i-1]tmp;if(tmp)Change(i,1);}memset(t,0,sizeof(t));Change(m,1);ll sumdp[m];for(ll im-1;i1;i--){bool tmp(Query(i1,iac[i1]1)!0);if(tmp)sumdp[i],Change(i,1);}memset(t,0,sizeof(t));limn;Change(1,1);for(ll i2;in;i){bool tmp(Query(i-cr[i]-1,i-1)!0);dp[i]dp[i-1]tmp;if(tmp)Change(i,1);}memset(t,0,sizeof(t));Change(n,1);ll ansdp[n];for(ll in-1;i1;i--){bool tmp(Query(i1,icr[i1]1)!0);if(tmp)ansdp[i],Change(i,1);}printf(%lld\n,ans*sum);return 0; }