ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

CF1474F 1 2 3 4 ...

2021-08-07 20:34:22  阅读:143  来源: 互联网

标签:... 结尾 int 转折点 ++ CF1474F ls include


一、题目

点此看题

二、解法

理解本题给出序列的方式很重要,我们把它放在坐标轴上,那么我们枚举一个转折点,然后找它后面最高的转折点,求最大差值就可以求出最长上升子序列长度,关键在于求出子序列个数。

不难发现不同的最低点和最高点组合范围一定不交,如果相交可以通过调整获得更优解,所以每一段可以单独处理。

坐标的范围达到了 \(1e9\),肯定往矩阵加速上想,正常递推的顺序是按照 \(x\) 轴从左往右,可以记录以这个点结尾的方案数,但是硬要这么转移你需要记录每个 \(y\) 结尾的方案数,那这做锤子。

问题在于无法记录 \(y\),那么我们按 \(y\) 的顺序递推不就不需要记录它了吗?所以设计状态 \(f_{v,i}\) 表示第 \(i\) 段以 \(y=v\) 结尾有多少种方案,我们先把所有转折点的 \(y\) 离散化,然后对于一小段它们的转移方式是相同的,所以可以矩阵加速第一维,时间复杂度 \(O(n^4\log n)\)

注意每一段要保证左闭右闭,所以端点一定要特别注意,还有如果答案长度为 \(1\) 需要特判。

三、总结

走到绝境的时候应该思考,\(dp\) 状态、转移等是否合理,不要受困于定式思维。

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
#define int long long
const int M = 55;
const int MOD = 998244353;
int read()
{
	int x=0,f=1;char c;
	while((c=getchar())<'0' || c>'9') {if(c=='-') f=-1;}
	while(c>='0' && c<='9') {x=(x<<3)+(x<<1)+(c^48);c=getchar();}
	return x*f;
}
int n,m,k,s,ans,res,a[M],l[M],r[M],d[M*M];
struct mat
{
	int a[M][M];
	void clear() {memset(a,0,sizeof a);}
	mat() {clear();}
	mat operator * (const mat &b) const
	{
		mat r;
		for(int i=0;i<=m;i++)
			for(int j=0;j<=m;j++)
				for(int k=0;k<=m;k++)
					r.a[i][k]=(r.a[i][k]+a[i][j]*b.a[j][k])%MOD;
		return r;
	}
	void print()
	{
		for(int i=0;i<=m;i++,puts(""))
			for(int j=0;j<=m;j++)
				printf("%d ",a[i][j]);
	}
}A,B;
mat qkpow(mat a,int b)
{
	mat r;
	for(int i=0;i<=m;i++) r.a[i][i]=1;
	while(b>0)
	{
		if(b&1) r=r*a;
		a=a*a;
		b>>=1;
	}
	return r;
}
void solve(int L,int R)
{
	m=k=0;
	for(int i=L,t=s;i<=R;t+=a[i],i++)
	{
		l[++m]=(a[i]>0)?t+1:t-1;
		r[m]=t+a[i];
		if(m==1) l[1]=s;
		d[++k]=l[m],d[++k]=l[m]-1,d[++k]=l[m]+1;
		d[++k]=r[m],d[++k]=r[m]-1,d[++k]=r[m]+1;
	}
	sort(d+1,d+1+k);
	A.clear();A.a[0][0]=1;
	k=unique(d+1,d+1+k)-d-1;
	for(int i=1,ls=s-1;i<=k;i++)
	{
		B.clear();
		if(d[i]<=ls || d[i]>s+ans) continue;
		for(int j=1;j<=m;j++)
			if(min(l[j],r[j])<=ls+1 && max(l[j],r[j])>=d[i])
			{
				for(int k=0;k<=j;k++)
					B.a[k][j]=1;
				if(l[j]>=r[j]) B.a[j][j]=0;
			}
		A=A*qkpow(B,d[i]-ls);ls=d[i];
	}
	for(int i=0;i<=m;i++)
		res=(res+A.a[0][i])%MOD;
}
signed main()
{
	n=read();read();
	for(int i=1;i<=n;i++)
	{
		a[i]=read();
		if(!a[i]) {i--;n--;continue;}
	}
	for(int i=1;i<=n;s+=a[i],i++)
		for(int j=i,t=s;j<=n;j++)
		{
			t+=a[j];
			ans=max(ans,t-s);
		}
	if(!ans)
	{
		printf("1 %lld\n",(-s+1)%MOD);
		return 0;
	}
	s=0;
	for(int i=1;i<=n;s+=a[i],i++)
	{
		int r=-1;
		for(int j=i,t=s;j<=n;j++)
		{
			t+=a[j];
			if(t-s==ans) r=j;
		}
		if(r>0) solve(i,r),i=r;
	}
	printf("%lld %lld\n",ans+1,res);
}

标签:...,结尾,int,转折点,++,CF1474F,ls,include
来源: https://www.cnblogs.com/C202044zxy/p/15113113.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有