spoj2817 线段树+离散化+动态规划

题目和2815差不多，只不过对于两个完全相同的序列，即使相对位置不完全相同，那也不能算作两个序列，问这样的序列的个数。

相较于上一题，这道题加入一个mark数组记录相同的一个数在之前是否出现过，如果出现过，那么对于以其为结尾的长度为k的序列，以当前位置的数为结尾肯定包含了以之前的数为结尾的序列，即例如1 2  5 4 5，那么第一个出现的5长度为3的序列为1 2 5，第二次出现的5长度为3的序列为1 2 5和1 4 5包含了1 2 5，那么只需省略掉之间的结果便可以得出除去重复计算部分的解。具体解法看spoj2815，此题只给出与2815不同部分的解法。

1A代码：

#include <iostream>
#include <cstring>
#include <algorithm>
#include <map>
#include <cstdio>

using namespace std;

#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define delf int m=(l+r)>>1
int sum[51][40001];
int n,m,s;
const int mod=5000000;
int mark[10001];

void build(int k,int l,int r,int rt)
{
    sum[k][rt]=0;
    if (l==r)
        return ;
    delf;
    build(k,lson);
    build(k,rson);
    return ;
}

void pushup(int k,int rt)
{
    sum[k][rt]=(sum[k][rt<<1]+sum[k][rt<<1|1])%mod;
}

int query(int k,int v,int l,int r,int rt)
{
    if (v<1)
        return 0;
    if (r<=v)
        return sum[k][rt]%mod;
    delf;
    if (m>=v)
        return query(k,v,lson)%mod;
    else
        return (query(k,v,lson)+query(k,v,rson))%mod;
}

void update(int p,int k,int c,int v,int l,int r,int rt)
{
    if (l==r)
    {
        if (p==1)
            sum[k][rt]=0;
        sum[k][rt]+=c;
        return ;
    }
    delf;
    if (m>=v)
        update(p,k,c,v,lson);
    else
        update(p,k,c,v,rson);
    pushup(k,rt);
    return ;
}

int main()
{
    while (scanf("%d%d",&n,&m)!=EOF)
    {
        int a[10001],b[10001];
        for (int i=1;i<=n;i++)
        {
            scanf("%d",&a[i]);
            b[i]=a[i];
        }
        sort(b+1,b+n+1);
        s=unique(b+1,b+n+1)-(b+1);
        /*
        map <int,int> mymap;
        for (int i=1;i<=s;i++)
            mymap[b[i]]=i;
        */
        for (int i=1;i<=m;i++)
            build(i,1,s,1);
        memset(mark,0,sizeof(mark));
        for (int i=1;i<=n;i++)
        {
            int j=lower_bound(b+1,b+s,a[i])-b;
            int p=0;
            if (mark[j]==1)
                p=1;
            //cout<<j<<endl;
            update(p,1,1,j,1,s,1);
            for (int t=2;t<=m;t++)
            {
                //cout<<mymap[a[i]]<<endl;
                int c=query(t-1,j-1,1,s,1)%mod;
                update(p,t,c,j,1,s,1);
                //cout<<"ANS: "<<sum[m][1]<<endl;
            }
            mark[j]=1;
            //cout<<endl;
        }
        printf("%d\n",sum[m][1]);
    }
}