用户工具

站点工具


2020-2021:teams:legal_string:jxm2001:contest:牛客练习赛83

牛客练习赛83

D-数列递推

题意

给定 $n,f_0$,求 $f_{1\sim n}$,其中

$$ f_i=\sum_{j=1}^if_{i\mod j} $$

题解

设 $i=kb+r$,考虑整数分块枚举 $k$,设 $t=i\mod k$,此时有 $r=t,t+k,t+2k,\cdots i-k\ast\text{lef}$。

当 $k\lt \sqrt n$ 时,如果能提前维护 $f_t,f_{t+2k},f_{t+3k}\cdots$ 的前缀和,就可以 $O(1)$ 计算贡献。

当 $k\ge \sqrt n$ 时,显然 $b,r$ 唯一,直接计算贡献即可。于是整数分块部分的时间复杂度为 $O(\sqrt i)$。

对每个 $i$,枚举 $k\lt \sqrt n$,更新 $f_t,f_{t+2k},f_{t+3k}\cdots$ 的前缀和的时间复杂度为 $O(\sqrt n)$,于是总时间复杂度 $O(n\sqrt n)$。

查看代码

查看代码

const int MAXN=1e5+5,MAXM=405,Mod=998244353;
int f[MAXN];
int s[MAXM][MAXN];
int main()
{
	int n=read_int(),v0=read_int(),m=sqrt(n)+1;
	f[0]=1;
	_for(i,0,m)
	s[i][0]=1;
	_rep(i,1,n){
		int lef=1,rig=0;
		while(lef<=i){
			rig=i/(i/lef);
			int k=i/lef;
			if(k<m)
			f[i]=(f[i]+s[k][i-lef*k])%Mod;
			else
			f[i]=(f[i]+f[i-lef*k])%Mod;
			lef=rig+1;
		}
		_for(k,1,m){
			if(i>=k)
			s[k][i]=(s[k][i-k]+f[i])%Mod;
			else
			s[k][i]=f[i];
		}
	}
	_rep(i,1,n){
		f[i]=(1LL*f[i]*v0%Mod+Mod)%Mod;
		space(f[i]);
	}
	return 0;
}

E-小L的疑惑

题意

给定互素的 $a,b$,求 $ax+by(x,y\ge 0)$ 不能表示的数中第 $k$ 大的数。

题解

首先给定结论当 $a,b$ 互素时,$ax+by(x,y\ge 0)$ 不能表示的正数等价于所有形如 $ab-na-mb(n,m\ge 1)$ 的正数,具体见 证明

接下来问题转化为求所有形如 $ab-na-mb(n,m\ge 1)$ 的正数中第 $k$ 大的元素。

考虑维护两个队列,队列一维护所有 $ab-a-mb(m=1)$,且保证递减。

队列二维护所有 $ab-na-mb(n\ge 1,m\ge 2)$,且保证递减。每次选择当前两个队列的队首的较大者进行 $\text{pop}$,连续 $k$ 次即可得到第 $k$ 大。

队列一的维护仅需要每次 $\text{pop}$ 后 $m$ 加一即可。队列二初始时为空,每次 $\text{pop}$ 元素 $t$ 后将 $t-b$ 加入队列。

现证明这样得到的队列二一定满足递减性质。假设先前 $\text{pop}$ 的所有元素确实为前 $k$ 大,记为 $a_1,a_2\cdots a_k$。

记本次 $\text{pop}$ 的元素为 $a_{k+1}$。则队列二的当前所有元素一定是由 $a_t-b$ 得到的,由于 $a_t\gt a_{k+1}$,所以 $a_t-b\gt a_{k+1}-b$,满足单调性。

于是总时间复杂度 $O(k)$。另外假定 $a\gt b$,可以考虑二分答案,然后枚举 $n$,然后统计 $m$,时间复杂度 $O(\sqrt k\log (ab))$。

查看代码

查看代码

int main()
{
	LL a=read_int(),b=read_int(),k=read_int(),pos=1;
	queue<LL> q;
	while(true){
		if(q.empty()||q.front()<a*b-pos*a-b){
			k--;
			if(k==0){
				enter(a*b-pos*a-b);
				break;
			}
			q.push(a*b-pos*a-2*b);
			pos++;
		}
		else{
			k--;
			if(k==0){
				enter(q.front());
				break;
			}
			q.push(q.front()-b);
			q.pop();
		}
	}
	return 0;
}
2020-2021/teams/legal_string/jxm2001/contest/牛客练习赛83.txt · 最后更改: 2021/05/23 10:01 由 jxm2001