RMQ——ST算法

template<typename T>class RMQ{  T **f;  int n1;  int n2;  bool (*cmp)(const T &a,const T &b);  bool compare(const T &a,const T &b)  {    if(cmp==NULL)      return a < b;    return cmp(a,b);  }  void init()  {    for(int i=1;i<=n1;i++)    {      for(int j=0;j<n2-(1<<i)+1;j++)      {        f[i][j]=compare(f[i-1][j],f[i-1][j+(1<<(i-1))])?f[i-1][j]:f[i-1][j+(1<<(i-1))];      }    }  }  int log2(int p)  {    int ans=0;    while(p)    {      ans++;      p>>=1;    }    return ans-1;  }  public:    RMQ():n1(0),n2(0),f(NULL),cmp(NULL){}    RMQ(const T *p,const T *q)    {      n2=q-p;      n1=log2(n2);      cmp=NULL;      f=new T*[n1+1];      for(int i=0;i<=n1;i++)        f[i]=new T[n2];      for(int i=0;i<n2;i++)        f[0][i]=*(p+i);      init();    }    RMQ(const T *p,const T *q,bool (*cm)(const T &a,const T &b))    {      n2=q-p;      n1=log2(n2);      cmp=cm;      f=new T*[n1+1];      for(int i=0;i<=n1;i++)        f[i]=new T[n2];      for(int i=0;i<n2;i++)        f[0][i]=*(p+i);      init();    }    RMQ(int p,T mnt[]):n1(log2(p)),n2(p),cmp(NULL)    {      f=new T*[n1+1];      for(int i=0;i<=n1;i++)        f[i]=new T[n2];      for(int i=0;i<n2;i++)        f[0][i]=mnt[i];      init();    }    RMQ(int p,T mnt[],bool (*cm)(const T &a,const T &b)):n1(p),n2(log2(p)),cmp(cm)    {      f=new T*[n1+1];      for(int i=0;i<=n1;i++)        f[i]=new T[n2];      for(int i=0;i<n2;i++)        f[0][i]=mnt[i];      init();    }    RMQ(const RMQ &rm)    {      *this=rm;    }    RMQ &operator=(const RMQ &rm)    {        if(this==&rm)            return *this;      cmp=rm.cmp;      if(f==NULL||n2!=rm.n2)      {        n1=rm.n1;        n2=rm.n2;        f=new T*[n1+1];        for(int i=0;i<=n1;i++)          f[i]=new T[n2];      }      for(int i=0;i<=n1;i++)        for(int j=0;j<n2;j++)          f[i][j]=rm.f[i][j];      return *this;    }    ~RMQ()    {      clear();    }    void clear()    {      if(f==NULL)        return ;      for(int i=0;i<=n1;i++)        delete[] f[i];      delete[] f;    }    int size(){return n2;}    T query(int p,int q)    {      int len=q-p+1;      int l=log2(len);      return compare(f[l][p],f[l][q-(1<<l)+1])?f[l][p]:f[l][q-(1<<l)+1];    }};