2023“钉耙编程”联赛 Day 3 E 题 - Out of Control 题解

原题描述

题面由我校翻译组激情制造。

有一个云服务 $\text{API}$ 来帮助用户存储历史时间戳。

每个用户的数据结构最初是一个空堆栈。

每次向 $\text{API}$ 发送带有整数 $x$ (表示当前时间戳)的请求时，服务器都会将 $x$ 附加到堆栈的末尾。

若 $x$ 小于堆栈中存储的前一个时间戳时，服务器将认为输入是错误的，将加入前一个时间戳入堆栈而不是 $x$。

您已经发布了 $n$ 次 $\text{API}$，您的请求时间戳在第 $i$ 次调用中是 $x_i$。

然而，网络已经失控。

服务器可能以任意顺序接收您的请求，甚至可能错过某些请求。

知道了这个问题，您要求随叫随到的工程师查看数据库中的堆栈。

假设服务器正好收到 $k$ 个请求，那么它可能有多少不同的堆栈?

输入格式

第一行包含单个整数 $T$ $(1\le T\le 100)$，测试用例的数量。对于每个测试用例:

输入的第一行包含单个整数 $n$ $(1\le n\le 3000)$，表示请求总数。

第二行包含 $n$ 个整数 $x_1,x_2,\dots ,x_n$ $(1\le x_i\le 10^9)$，表示每个请求的时间戳。

可以保证所有 $n$ 的和不超过 $30000$。

输出格式

对于每个测试用例，输出 $n$ 行，其中第 $i$ $(1\le i\le n)$ 行包含一个整数，表示 $k=i$ 时不同堆栈的数量。

请注意，答案可能非常大，所以请打印它模 $10^9+7$ 代替。

样例

样例输入

2
3
1 2 3
3
2 3 3

样例输出

3
5
5
2
2
2

样例解释

在第一个示例中：

当 $k=1$ 时，栈可以是 $[1],[2],[3]$。
当 $k=2$ 时，栈可以是 $[1,2],[1,3],[2,2],[2,3],[3,3]$。
当 $k=3$ 时，栈可以是 $[1,2,3],[1,3,3],[2,2,3],[2,3,3],[3,3,3]$。

在第二个示例中：
当 $k=1$ 时，栈可以是 $[2],[3]$。
当 $k=2$ 时，栈可以是 $[2,3],[3,3]$。
当 $k=3$ 时，栈可以是 $[2,3,3],[3,3,3]$。

题目分析

首先发现我们不关心值的大小，只关心排名，而且貌似需要嗯造到动态规划状态里，所以需要离散化。

我们考虑转换问题，发现类似于填充一个不下降子序列，但是有点限制。

接下来考虑动态规划，我们设 $f_{i,j}$ 表示填到第 $i$ 位，末尾为 $j$。

边界比较清楚，$f_{1,j}=1$。

我们发现这个填充并不是没有条件的，事实上，我们发现 $i\le t_j$，$t_j$ 表示有多少数字小于等于 $j$，因为我们在前面填出的多的 $j$ 本质是拿比 $j$ 小的数字代替，因此怎么填也填不到 $t_j$ 后面。

我们不难得到转移方程：

$$f_{i,j}=\sum _{k=1}^j{f}_{i-1,k}(i\le t_j)$$

条件不成立则值为 $0$。

不得不说，这个还是比较好想的，但是我们需要优化，加个前缀和即可。

我们用定义一个 $sum$，有：

$$su{m}_{i,j}=\sum _{k=1}^j{f}_{i,j}$$

以上方程可以写成：

$$f_{i,j}=su{m}_{i-1,j}(i\le t_j)$$

对于第 $i$ 个输出：

$$an{s}_i=\sum _{j=1}^M{f}_{i,j}$$

这里的 $M$ 表示离散化之后的最大值。

实现代码

#include<bits/stdc++.h>
#define LL long long
using namespace std;
const LL N=3e3+5;
const LL mod=1e9+7;
LL T,n,x[N],a[N],t[N],f[N][N],sum[N][N],cnt,ans[N];
int main()
{
	scanf("%lld",&T);
	while(T--)
	{
		scanf("%lld",&n);
		memset(t,0,sizeof(t));
		memset(f,0,sizeof(f));
		memset(sum,0,sizeof(sum));
		memset(ans,0,sizeof(ans));
		for(int i=1;i<=n;i++)
		{
			scanf("%lld",&x[i]);
		}
		sort(x+1,x+n+1);
		cnt=0;
		for(int i=1;i<=n;i++)
		{
			if(x[i]!=x[i-1])cnt++;
			a[i]=cnt;
			t[cnt]++;
		}

		for(int i=1;i<=cnt;i++)t[i]+=t[i-1];
		for(int i=1;i<=n;i++)
		{
			for(int j=1;j<=cnt;j++)
			{
				if(i==1)f[i][j]=1;
				else if(i<=t[j])f[i][j]=sum[i-1][j];
				else f[i][j]=0;
				ans[i]=(ans[i]+f[i][j])%mod;
				sum[i][j]=(sum[i][j-1]+f[i][j])%mod;
			}
		}
		for(int i=1;i<=n;i++)
		{
			printf("%lld\n",ans[i]);
		}
	}
}

如果你看完了整个文章，下面有投票可以玩，再奖励你看一个科比。