回文子序列数 编程之美 hihocoder 资格赛 第二题 2015年

http://hihocoder.com/contest/msbop2015qual/problem/2

 (ps:该的代码在30行的时候，会出现负数。也就是有问题。优化见下方)

O(n^3)

可以进一步优化为O(n^2)->通过优化cal函数

 

 

#include <iostream>
#include <string>
#include <stdio.h>
#include <memory>
#include <string.h>
using namespace std;
long long map[1001][1001];
string input;
void cal(int start ,int end)
{
    if(start==end)
    {
        map[start][end]=1;
        return ;
    }
    if(end-start==1)
    {
        if(input[start]==input[end])
            map[start][end]=3;
        else
            map[start][end]=2;
        return ;
    }
    map[start][end]+=map[start][end-1]+1;
    for(int i=end-1;i>=start;i--)
    {
        if(input[i]==input[end])
            map[start][end]+=map[i+1][end-1]+1;
    }
    map[start][end]%=100007;
}
int main()
{
    int T;
    int Case=0;
    cin>>T;
    
    while(T--)
    {
        Case++;
        memset(map,0,sizeof(map));
        
        cin>>input;
        for(int i=0;i<input.size();i++)
        {
            for(int j=i;j>=0;j--)
            {
                cal(j,i);
            }
        }
        cout<<"Case #"<<Case<<": "<<map[0][input.size()-1]<<endl;
    }

}

 

潇湘缘说的很对，以上代码是O(n^3)的，在面对最大数据1000的时候基本就是超时了的：

以下是更正代码：优化了cal()函数，把map[start+1][end]充分利用起来，是O(n^2)的时间开销了

(ps:为什么耗时还是19ms。。)

(ps:之前的代码因为少了37行的过程，输出会有负数，也就是有问题。)

以下代码已通过大数据测试

#include <iostream>
#include <string>
#include <stdio.h>
#include <memory>
#include <string.h>
using namespace std;
long long map[1001][1001];
string input;



int main()
{
    int T;
    int Case=0;
    cin>>T;

    while(T--)
    {
        Case++;
        memset(map,0,sizeof(map));
        cin>>input;
        for(int i=0;i<input.size();i++)
        {
            for(int j=i;j>=0;j--)
            {
                if(j==i)
                {
                    map[j][i]=1;
                    continue;
                }
                map[j][i]=map[j+1][i]-map[j+1][i-1]+map[j][i-1];
                if(input[j]==input[i])
                {
                    map[j][i]+= map[j+1][i-1]+1;
                }
                map[j][i]=(map[j][i]+100007)%100007;
            }
        }
        cout<<"Case #"<<Case<<": "<<map[0][input.size()-1]<<endl;
    }

}