Codeforces 245H Queries for Number of Palindromes

题意：给你一个字符串，再给你一个q（询问个数1000000），问你这个区间内回文串的个数。

解题思路：

1）dp，先把 i 到j 是回文串的真值赋值给 dp[i][j] ，然后从后往前dp    dp[i][j] += dp[i+1][j] + dp[i][j-1]  -dp[i+1][j-1];

解题代码：

// File Name: 245h.dp.cpp
// Author: darkdream
// Created Time: 2015年03月30日 星期一 15时27分10秒

#include<vector>
#include<list>
#include<map>
#include<set>
#include<deque>
#include<stack>
#include<bitset>
#include<algorithm>
#include<functional>
#include<numeric>
#include<utility>
#include<sstream>
#include<iostream>
#include<iomanip>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<ctime>
#define LL long long
#define maxn 5005
using namespace std;
int dp[maxn][maxn];
char str[maxn];
int n; 
int main(){
    scanf("%s",&str[1]);
    n = strlen(&str[1]);
    for(int i = 1;i <= n;i ++)
    {
      if(str[i] == str[i+1]) 
          dp[i][i+1] = 1;
      dp[i][i] = 1; 
    }
    for(int i = n;i >= 1;i --)
    {
       for(int j = i+1;j <= n;j ++)
       {
          if(dp[i+1][j-1] == 1&& str[i] == str[j]) 
          {
            dp[i][j] = 1; 
          }
       }
    }
    for(int i = n;i >= 1;i --)
      for(int j = i ;j <= n ;j ++)
      {
        dp[i][j] += dp[i+1][j] + dp[i][j-1] -dp[i+1][j-1];
      }
    int q; 
    scanf("%d",&q);
    while(q--)
    {
      int l , r ;
      scanf("%d %d",&l,&r);
      printf("%d\n",dp[l][r]);
    }
return 0;
}

2）树状数组 + 离线  。。。 发现线段树过不了  毕竟我写的线段树常数太大。

解题代码：

// File Name: 245h.cpp
// Author: darkdream
// Created Time: 2015年03月30日 星期一 11时09分19秒

#include<vector>
#include<list>
#include<map>
#include<set>
#include<deque>
#include<stack>
#include<bitset>
#include<algorithm>
#include<functional>
#include<numeric>
#include<utility>
#include<sstream>
#include<iostream>
#include<iomanip>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<ctime>
#define LL long long

using namespace std;
char  str[6000];
map<int ,int >dpj[5005];
map<int ,int >dpo[5005];
map<int ,int >::iterator tp;
struct qu{
 int l , r,si; 
}a[1000005];
int ans[1000005];
int q, n ; 
int cmp(qu a, qu b)
{
  return a.l > b.l ;
}
int tree[6005];
int lowbit(int x)
{
  return x & -x;
}
void update(int x)
{
   while(x <= n)
   {
      tree[x] += 1; 
      x += lowbit(x);
   }
}
int getsum(int r)
{
  int ans  = 0 ; 
  while(r >= 1)
  {
     ans += tree[r];
     r -= lowbit(r);
  }
  return ans;
}
inline void read_int(int &ret){
   char ch = getchar();
   while(!isdigit(ch)){
     ch = getchar();
   }
   ret =0 ; 
   while(isdigit(ch)){
      ret = ret*10 + ch -'0' ;
      ch = getchar();
   }
}
int main(){
 //   freopen("input","r",stdin);
//    freopen("output","w",stdout);
    clock_t be ,en ;
    be = clock();
    scanf("%s",&str[1]);
     scanf("%d",&q);
     for(int i = 1;i <=q ;i ++)
     {
         read_int(a[i].l);
         read_int(a[i].r);
         a[i].si = i ; 
     }
     sort(a+1,a+1+q,cmp);
     n = strlen(&str[1]);
     int tt = 1;  
     for(int i = n;i >= 1;i --)
     {
        dpj[i][1] = i;
        update(i);
        if(str[i] == str[i+1])
        {
          dpo[i][2] = i+1;
          update(i+1);
        }
        for(tp = dpj[i+1].begin() ; tp != dpj[i+1].end();tp++)
        {
           if(str[i] == str[tp->second +1])
           {
               dpj[i][tp->first+2] = tp->second +1;
               update(tp->second+1);
           }
        }
        for(tp = dpo[i+1].begin() ; tp != dpo[i+1].end();tp++)
        {
           if(str[i] == str[tp->second +1])
           {
               dpo[i][tp->first+2] = tp->second +1;
               update(tp->second+1);
           }
        }
        while(a[tt].l == i && tt <= q)
        {
           ans[a[tt].si] = getsum(a[tt].r);
           tt++;
        }
        if(tt >q)
            break;
        dpj[i+1].clear();
        dpo[i+1].clear();
     }
     for(int i = 1;i <= q; i ++)
         printf("%d\n",ans[i]);
     en = clock();
//     printf("%lf\n",(en -be)*1.0/CLOCKS_PER_SEC);
return 0;
}