[AHOI2009] BZOJ2431 逆序对数列-动态规划-前缀和优化

121 篇文章 0 订阅
40 篇文章 0 订阅

传送门

题目大意:求逆序对数为k的长为n的排列有多少?n,k<=1000

题解:这显然是个DP(废话)

而且dp[n][k]表示长为n逆序对数为k的方案数。

而且状态转移显然要从dp[n-1][...]推过来。

考虑如果把n插入到1~(n-1)的某个逆序对数为x的排列的第p个数的后面,

那么会新产生(n-1)-(p+1)+1=n-p-1对新的逆序对,其中0<=p<=n-1是显然的。

因此如果当前状态dp[n][k],那么上一个状态就是dp[n-1][x=k-(n-p-1)=p+k+1-n]。

那么显然要求p+k+1-n>=0,因此p>=n-k-1,因此max(0,n-k-1)<=p<=n-1。

状态转移方程为:dp[n][k]=sum{dp[n-1][p+k+1-n]},max(0,n-k-1)<=p<=n-1。

状态是O(n^2)的,转移是O(n)的,O(n^3)肯定要T。

优化其实是显然的,因为要求和,没有限制条件,还是连续区间,显然用前缀和优化。

具体来说,将max(0,n-k-1)<=p<=n-1代入到状态转移方程中“来源状态”的第二维,

可以得到dp[n][k]对于dp[n-1][..]的求和区间是[max(0,n-k-1)+k+1-n,k]。

所以再用sumn[n][k]表示sum{dp[n][x]},0<=x<=k。

那么dp[n][k]=sumn[n-1][k]-sumn[n-1][max(0,n-k-1)+k-n]。

这里实现有一个细节,就是第二个sumn的第二维有可能是-1,这个时候sumn=0,要特判一下。

另外一个细节是,对于每一个sumn[n][比较大的...],后面比较大的...也要算全了,即使后面的dp[n][比较大的...]=0(不存在),因为状态转移需要。

这样转移的时间复杂度降到O(1),总复杂度O(N^2)。

我的代码就写到这里了,又去看了看别人的做法,其实可以优化空间复杂度和编码复杂度(一会看看我的代码就知道我的代码有多乱和容易写错了)

空间复杂度O(N^2)看起来很不清真啊,首先由于dp[n]和sumn[n]都是由n-1转移来的,所以可以滚动数组优化到O(N)。

(注意并不是dp[N][2],而是dp[N],可以滚动到真的只需要一维,具体实现留给大家思考)。

但编程还是要考虑上文提到的两个细节。注意到求和区间是非减的,所以可以用一个变量代替sumn数组,实现方式类似于可行性完全背包分组+单调队列的做法,编程复杂度特别低。而且不需要考虑上述两个细节。

代码:

//BZOJ 2431
//AHOI 2009
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define MAXN 1010
#define mod 10000
using namespace std;
int dp[MAXN][MAXN],sumn[MAXN][MAXN];
int main()
{
	int n,k;scanf("%d%d",&n,&k);
	dp[1][0]=sumn[1][0]=sumn[1][1]=1;
	for(int i=2;i<=n;i++)
	{
		dp[i][0]=sumn[i][0]=1;
		for(int j=1;j<=min(k,i*(i-1)/2);j++)
		{
			int a=max(i-j-1,0),b;
			if(j-i+a<0) b=0;
			else b=sumn[i-1][j-i+a];
			dp[i][j]=(sumn[i-1][j]-b)%mod;
			sumn[i][j]=(sumn[i][j-1]+dp[i][j])%mod;
		}
		for(int j=min(k,i*(i-1)/2)+1;j<=min(k,i*(i+1)/2);j++)
			sumn[i][j]=sumn[i][j-1];
	}
	dp[n][k]=(dp[n][k]+mod)%mod;
	printf("%d\n",dp[n][k]);return 0;
}

  • 1
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
很抱歉,我似乎没有找到与彩旗飘飘Python相关的引用内容。但是,如果你想在Python中实现彩旗飘飘的效果,你可以使用Pygame库来创建动画。以下是一个简单的示例代码,可以帮助你开始: ```python import pygame import random # 初始化pygame pygame.init() # 设置屏幕尺寸 screen_width = 800 screen_height = 600 screen = pygame.display.set_mode((screen_width, screen_height)) pygame.display.set_caption("彩旗飘飘") # 加载彩旗图片 flag_image = pygame.image.load("flag.png") # 获取彩旗图片的宽度和高度 flag_width = flag_image.get_width() flag_height = flag_image.get_height() # 设置彩旗的起始位置 flag_x = screen_width // 2 - flag_width // 2 flag_y = screen_height // 2 - flag_height // 2 # 设置彩旗的速度 flag_speed_x = random.randint(-5, 5) flag_speed_y = random.randint(-5, 5) # 游戏主循环 running = True while running: for event in pygame.event.get(): if event.type == pygame.QUIT: running = False # 移动彩旗 flag_x += flag_speed_x flag_y += flag_speed_y # 如果彩旗碰到边界,反转速度方向 if flag_x <= 0 or flag_x + flag_width >= screen_width: flag_speed_x = -flag_speed_x if flag_y <= 0 or flag_y + flag_height >= screen_height: flag_speed_y = -flag_speed_y # 清空屏幕 screen.fill((255, 255, 255)) # 绘制彩旗 screen.blit(flag_image, (flag_x, flag_y)) # 更新屏幕 pygame.display.flip() # 退出游戏 pygame.quit() ``` 请注意,你需要将代码中的"flag.png"替换为你自己的彩旗图片。此外,你还可以根据需要调整彩旗的起始位置、速度等参数。希望这个示例代码能帮助到你!<span class="em">1</span><span class="em">2</span> #### 引用[.reference_title] - *1* [3d max制作彩旗飘飘](https://download.csdn.net/download/m0_71585230/88011500)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_2"}}] [.reference_item style="max-width: 50%"] - *2* [[Rqnoj-371][AHOI1997]彩旗飘飘](https://blog.csdn.net/w745241408/article/details/7176600)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_2"}}] [.reference_item style="max-width: 50%"] [ .reference_list ]

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值