codeforces 628e

好像是被我水过了    因为题目给了5s

题意：问图像中到底有多少个z型图像。

思路：对于每一个点都有可能是最左下角的点，那么把对于每一个 斜边 看成大区间  区间每一个点都是都有两个值 一个是能到右边的最远距离 一个是到左边的最远距离+i（加权）

只要对于每一个点的答案是把右边的最远距离作为区间长度 统计出区间 左边是否能满足

比如 3x3

zzz

.z.

zzz

对于 最左下角的点 统计区间长度为3  加权距离为行数+左边的最远距离 为4 然后区间内 只要大于4的数 就是满足条件的

就是转换成问这个数在区间内 是第几大的问题  就抄了别人模版=。=

ps：这个方法4s多才过，虽然我代码写搓 也有部分原因- -，不过  他们都是1s内过的  惭愧

#include<stdio.h>
#include<string.h>
#include<vector>
#include<algorithm>
#include<queue>
#include<math.h>
using namespace std;
#define LL long long
#define MOD 1000000007

int a[3005][3005];

int c[3005][3005];
int d[3005][3005];
char ch[3005][3005];
#define CLR(a) memset(a,0,sizeof(a))

const int MAXN = 3010;

struct Node
{
  int l,r;
  int len(){return r-l;}
  int mid(){return (l+r)/2;}
  bool in(int ll,int rr){return l>=ll && r<=rr;}
  void lr(int ll,int rr){l = ll; r=rr;}
}node[MAXN*4];

int num_left[15][MAXN],seg[15][MAXN],sa[MAXN];

void Init()
{
  //CLR(seg),CLR(num_left),CLR(node);
    memset(num_left,0,sizeof(num_left));
    memset(node,0,sizeof(node));
}

void PaBuild(int t,int l, int r,int d)
{
  node[t].lr(l,r);
  if(node[t].len() == 0)return;
  int mid=(l+r)/2,lsame=mid-l+1;
  for(int i=l;i<=r;i++)//首先确定分到每一侧的数的数目
    if(seg[d][i] < sa[mid])//因为相同的元素可能被分到两侧
      lsame--;
  int lpos=l,rpos=mid+1;
  for(int i=l;i<=r;i++)
  {
    if(i == l)
      num_left[d][i]=0;
    else
      num_left[d][i]=num_left[d][i-1];
    if(seg[d][i] < sa[mid])
    {
      num_left[d][i]++;
      seg[d+1][lpos++]=seg[d][i];
    }
    if(seg[d][i] > sa[mid])
      seg[d+1][rpos++] = seg[d][i];
    if(seg[d][i] == sa[mid])
      if(lsame > 0)// 如果大于0，说明左侧可以分和sa[mid]相同的数字
      {
        lsame--;
        num_left[d][i]++;
        seg[d+1][lpos++]=seg[d][i];
      }
      else	//反之，说明左侧数字已经分满了，就分到右边去
        seg[d+1][rpos++]=seg[d][i];
  }
  PaBuild(t*2,l,mid,d+1);
  PaBuild(t*2+1,mid+1,r,d+1);
}
void Build(int s,int t)
{
  sort(sa+s,sa+t+s);
  PaBuild(1,s,t,1);
}

int find_rank(int t,int l,int r,int d,int val)//val指的是k
{
  if(node[t].len() == 0)return seg[d][l];
  int s,ss;   //s表示区间[l,r]有多少个小于sa[mid]的数被分到左边
  if( l == node[t].l)//ss表示从当前区间的L到l-1有多少个小于sa[mid]的数被分到左边,L,R指的是树上当前节点的区间范围
    ss=0;
  else
    ss=num_left[d][l-1];
  s=num_left[d][r]-ss;
  if(s>=val)
    return find_rank(t*2, node[t].l+ss,node[t].l+ss+s-1,d+1,val);
  else
  {
    int mid = node[t].mid();
    int bb=l-node[t].l-ss;  //表示从当前区间L到l-1有多少个分到右边
    int b=r-l+1-s;		  //表示[l,r]有多少个分到右边
    return find_rank(t*2+1,mid+bb+1,mid+bb+b,d+1,val-s);
  }
}

int Query(int s,int t,int k)
{
  return find_rank(1,s,t,1,k);
}


int main()
{
    int n,m;
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
        scanf("%s",ch[i]+1);
    n=max(n,m);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
        {
            if(ch[i][j]=='z')
                a[i][j]=1;
        }
    for(int i=1;i<=n;i++)
    {
        int sum=0;
        for(int j=n;j>=1;j--)
        {
            if(a[i][j])
                sum++;
            else
                sum=0;
            c[i][j]=sum;
        }
    }
    for(int i=1;i<=n;i++)
    {
        int sum=0;
        for(int j=1;j<=n;j++)
        {
            if(a[i][j])
                sum++;
            else
                sum=0;
            d[i][j]=sum;
        }
    }
    LL ans=0;
    for(int i=1;i<2*n;i++)
    {
        int len=1;
        Init();
        for(int j=1;j<=n;j++,len++)
        {
            int x=i-j+1,y=j;
            if(i-j+1<=0)break;
            if(i-j+1>n){
                sa[len]=seg[1][len]=0;
                continue;
            }
            sa[len]=seg[1][len]=d[y][x]+j;
        }
        Build(1,len-1);
        int h=1;
        for(int j=1;j<=n;j++)
        {
            if(i-j+1<=0)break;
            if(i-j+1>n)continue;
            int x=i-j+1,y=j;
            if(d[y][x]==0){h=j+1;continue;}

            int t=max(j-c[y][x]+1,h);
            int l=1,r=j-t+1;
            int ret=r;

            while(l<=r)
            {
                int mid=(l+r)/2;
                int q=Query(t,j,mid);
                if(q>=j+1)
                {
                    r=mid-1;
                    ret=mid;
                }
                else
                {
                    l=mid+1;
                }
            }
            ans+=j-t+1-ret+1;
        }
    }
    printf("%lld\n",ans);
    return 0;
}