元月初七hu测 T3.秩序（dp）

Description

秩序带来美感，秩序带来高效。
在高度发达的三体世界，秩序的意义非常明显。
为了追求秩序，一些事物会被规范化，同时失去了它原有的性质，变得不
那么自然。
为了便于管理，三体世界的奶牛们工作完后都需要排队去食堂打饭。
一共有 m种食物，每头奶牛固定只吃一种食物。现在，n头奶牛排好队来
食堂就餐，食堂工作人员每次可以放进来队伍前端任意数量的奶牛，然后给它们做食物吃。当放进来的奶牛吃 k种不同的食物时，工作人员需要$k^2$的时间让它们吃饱。
追求秩序的世界就会追求效率。工作人员必须用最少的时间让所有奶牛都吃饱。他们需要知道最少需要多长时间。

Input

第一行是两个数n、m。
接下来n行，其中的第i行有一个数，表示队伍的第i位上的奶牛吃的食物编号，用1-m的整数表示。

Output

一行一个数，表示最少的时间。

Sample Input

13 4
1
2
1
3
2
2
3
4
3
4
3
1
4

Sample Output

11

Hint

30%的数据，1≤m≤n≤20
60%的数据，1≤m≤n≤5000
100%的数据，1≤m≤n≤50000

分析：
60%的数据，$n^2$直接shui过即可（注意空间。。。）

正解：
思路在$n^2$的级别上卡住了
看来想要进一步降低复杂度需要的是质的飞跃

100%的数据
第二维不用n个全枚举
枚举到种类数超过$sqrt(n)$就不枚举了
那样肯定不会更优，因为每个都单独一组才需要n的时间

骗人，只有80分：

int o=sqrt(n); f[0]=0;
for (int i=1;i<=n;i++) f[i]=INF; 
for (int i=1;i<=n;i++)
{
    memset(p,0,sizeof(p));
    int cnt=0;
    for (int j=i;j>=1;j--)
    {
        if (!p[a[j]]) p[a[j]]=1,cnt++;
        if (cnt>=o) break;
        f[i]=min(f[i],f[j-1]+cnt*cnt);
    }
}

上面的写法还是不够优秀
可以注意到，每次转移的时候，还有memset的复杂度（mamset是很慢的）
实际上我们可以发现，cnt是单调递增的

     0 1 2 3 4 5 6 7 8 9
  a: 1 3 1 4 2 2 1 3 1 2
//从位置2开始计算种数
cnt:     1 2 3 3 3 4 4 4

我们就可以直接找到这些cnt增加的断点
从这些断点转移即可

tip

这就是faebdc出的题？！
总的来说不是很困难，没有很毒瘤的题，但是都用到了一些很巧妙的优化吧

没有自己写AC代码（懒。。。）

//faebdc标程
//凑活着看吧您内
#include<cstdio>
#include<algorithm>
using namespace std;
const int N = 50500;
int n,m;
int kind[N];
int last[N],ever[N];
int dp[N];
int xu[N];
int main()
{
    freopen("order.in","r",stdin);
    freopen("order.out","w",stdout);
    int i,j;
    scanf("%d%*d",&n);
    for(m=1;m*m<=n;m++);
    m--;
    for(i=1;i<=n;i++)
        scanf("%d",&kind[i]);
    for(i=1;i<=n;i++)
    {
        last[i]=ever[kind[i]];
        ever[kind[i]]=i;
    }
    for(i=1;i<=n;i++)
    {
        if(last[i])
        {
            for(j=1;j<=xu[0];j++)
            {
                if(xu[j]==last[i])
                    break;
            }
            while(j<xu[0])
            {
                xu[j]=xu[j+1];
                j++;
            }
            xu[0]--;
        }
        xu[++xu[0]]=i;
        dp[i]=xu[0]*xu[0];
        int has=0;
        for(j=xu[0];j>1;j--)
        {
            has++;
            if(has>m)
                break;
            dp[i]=min(dp[i],dp[xu[j-1]]+has*has);
        }
    }
    printf("%d\n",dp[n]);
    return 0;
}