Codeforces Round #729 div.2 A-E题解

最新推荐文章于 2022-07-17 06:21:18 发布

欣君

最新推荐文章于 2022-07-17 06:21:18 发布

阅读量816

点赞数 13

分类专栏： ----Codeforces 文章标签：算法

本文链接：https://blog.csdn.net/xin_jun/article/details/118476221

版权

----Codeforces 专栏收录该内容

22 篇文章 8 订阅

订阅专栏

视频讲解：BV1hB4y1N7jF

A. Odd Set

题目大意

给定 $\leq n \leq 100)$ 个可能重复的整数 $a_i(0 \leq a_i \leq 100)$ ，问是否能将其分为 $n$ 对，使得每对的两个整数之和为奇数。

题解

只有当每对内包含一个奇数加一个偶数时，才能使其总和为奇数。因此充要条件为奇数个数与偶数个数相等。

参考代码

#include <bits/stdc++.h>
typedef long long ll;
using namespace std;

int main()
{
	int T,n,x,i,sum0,sum1;
	scanf("%d",&T);
	while(T--)
	{
		scanf("%d",&n);
		sum0=sum1=0;
		for(i=1;i<=2*n;i++)
		{
			scanf("%d",&x);
			if(x&1)
				sum1++;
			else
				sum0++;
		}
		if(sum1==sum0)
			printf("Yes\n");
		else
			printf("No\n");
	}
}

B. Plus and Multiply

题目大意

有一个无限大的集合，定义如下：

$1$ 在集合中；
如果 $x$ 在集合中，则 $\cdot a$ 和 $x + b$ 也在集合中；

给定 $\leq n,a,b \leq 10^9)$ ，判断 $n$ 是否在集合中。

题解

如果 $n$ 在集合中，那么 $n$ 必定可以表达为：
$n=a^x+b \cdot y$

其中 $\leq x,y$ 。因此只要枚举 $x$ ，判断是否存在合法的 $y$ 即可。
注意， $a$ 可能等于 $1$ ，因此需要特判，或限定枚举次数。

参考代码

#include <bits/stdc++.h>
typedef long long ll;
using namespace std;

int main()
{
	int T,i;
	ll n,a,b,base;
	bool flag;
	scanf("%lld",&T);
	while(T--)
	{
		scanf("%lld%lld%lld",&n,&a,&b);
		base=1;
		flag=false;
		for(i=1;i<=30&&base<=n;i++,base*=a)
		{
			if((n-base)%b==0)
			{
				flag=true;
				break;
			}
		}
		if(flag)
			printf("Yes\n");
		else
			printf("No\n");
	}
}

C. Strange Function

题目大意

设 $f (i)$ 表示最小的不是 $i$ 的因数的整数。
求 $\sum_{i=1}^n{f(i)}$ 对 $10^9+7$ 取模的结果。
其中 $\leq n \leq 10^{16}$ 。

题解

首先会发现， $f (i)$ 应该是一个相比 $i$ 而言，很小的数。
如果 $f (i) = y$ ，那么 $i$ 需要是 $L C M (1, 2, . . ., (y - 1))$ 的倍数，这个类似于阶乘的条件十分苛刻。
设 $a_i=LCM(1,2,3,...,i)$ ，可以通过递推打表的方式求出所有小于 $10^{16}$ 的 $a_i$ 。
若 $x$ 是 $a_i$ 的倍数，则有 $\geq i+1$ 。
因此可以通过枚举 $i$ ，统计有多少个不超过 $n$ 的 $a_i$ 的倍数，从而累加到答案上。
具体而言，
$ans=(n+\sum_{i=1}^{a_i \leq n}{\lfloor \frac{n}{a_i} \rfloor}) \% \;mod$

时间复杂度 $\cdot \log(n))$ 。

参考代码

#include <bits/stdc++.h>
typedef long long ll;
using namespace std;

const ll MAXN=1e16+100;
const ll mod=1e9+7;
ll a[50];

int main()
{
	int T,cnt,i;
	ll n,ans;
	a[0]=a[1]=1;
	for(cnt=1;a[cnt-1]<MAXN;cnt++)
	{
		for(i=1;i<=cnt;i++)
		{
			if(a[cnt-1]*i%cnt==0)
			{
				a[cnt]=a[cnt-1]*i;
				break;
			}
		}
	}
	scanf("%d",&T);
	while(T--)
	{
		scanf("%lld",&n);
		ans=0;
		for(i=0;i<cnt&&a[i]<=n;i++)
			ans=(ans+n/a[i])%mod;
		printf("%lld\n",ans);
	}
}

D. Priority Queue

题目大意

对于由 $+ x$ 和 $-$ 构成的操作序列 $S$ ，定义函数 $f (S)$ 由以下方式计算得到：

从左到右遍历 $S$ ，并维护一个多重集合 $T$ ；
对于 $S$ 中的每个元素，若其为 $+ x$ ，则添加 $x$ 到 $T$ 中；否则删除 $T$ 中最小的元素（若 $T$ 为空，则不进行任何操作）；
遍历完 $S$ 中的所有元素后，求 $T$ 中所有元素的总和，就是 $f (S)$ ;

给定一个长度为 $\leq n \leq 500)$ 的操作序列 $A$ ，求 $A$ 的所有子序列 $B$ 的 $f (B)$ 之和。

题解

对于每个 $a_x$ 的操作（假设是第 $x$ 次操作），考虑有多少种 $B$ 在遍历完后依旧保留 $a_x$ 。这样的子序列必须满足以下条件：

必须包含第 $x$ 次操作；
设 $s$ 为 $T$ 中小于 $a_x$ 的元素个数。对于第 $x$ 次之后的每一次 $-$ 操作前， $s$ 必须大于 $0$ ；

设 $dp_{i,j}$ 表示有多少种满足上述条件的子序列，并且：

最后一次操作不晚于操作 $i$ （若 $i < x$ 则不要求满足第一个条件）；
$T$ 有 $j$ 个小于 $a_x$ 的元素；

$dp_{i,j}$ 可以通过以下转移式得到：

若 $\neq x$ ，则可以不选择操作 $i$ ，因此 $dp_{i,j}+=dp_{i-1,j}$ ;
若 $i < x$ 且操作 $i$ 是 $-$ 操作，则可以选择操作 $i$ ，因此
- 若 $j = 0$ ，则上一步可以有 $1$ 个小于 $a_x$ 的元素，也可以没有小于 $a_x$ 的元素，即 $dp_{i,0}+=dp_{i-1,0}+dp_{i-1,1}$ ；
- 若 $j > 0$ ，则 $dp_{i,j}+=dp_{i-1,j+1}$
若 $i > x$ 且操作 $i$ 是 $-$ 操作，则上一步有 $j + 1$ 个小于 $a_x$ 的元素，即 $dp_{i,j}+=dp_{i-1,j+1}$ ；
若操作 $i$ 是 $a_i$ 操作，则判断 $a_i$ 和 $a_x$ 大小：
- 若 $a_i < a_x$ ，则 $dp_{i,j}+=dp_{i-1,j-1}$ ；
- 若 $a_i > a_x$ ，则 $dp_{i,j}+=dp_{i-1,j}$ ；
- 注意，若 $a_i=a_x$ ，还需判断 $i$ 和 $j$ 的大小，避免重复统计。

最终答案为
$ans=\sum_{op_x='+'}{\sum_{j}{dp_{n,j}\cdot a_x}}$

参考代码

#include <bits/stdc++.h>
typedef long long ll;
using namespace std;

const int MAXN=550;
const int mod=998244353;
char op[MAXN];
int a[MAXN],sum[MAXN];
ll dp[MAXN][MAXN];

int main()
{
	int n,i,j,x;
	ll ans;
	scanf("%d",&n);
	for(i=1;i<=n;i++)
	{
		scanf(" %c",&op[i]);
		sum[i]=sum[i-1]+(op[i]=='+');
		if(op[i]=='+')
			scanf("%d",&a[i]);
	}
	ans=0;
	for(x=1;x<=n;x++)
	{
		if(op[x]=='-')
			continue;
		memset(dp,0,sizeof(dp));
		dp[0][0]=1;
		for(i=1;i<=n;i++)
		{
			for(j=0;j<=sum[i];j++)
			{
				if(x!=i)
					dp[i][j]=dp[i-1][j];
				if(op[i]=='-')
				{
					if(i<x)
					{
						if(j==0)
							dp[i][j]=(dp[i][j]+dp[i-1][j]+dp[i-1][j+1])%mod;
						else
							dp[i][j]=(dp[i][j]+dp[i-1][j+1])%mod;
					}
					else
						dp[i][j]=(dp[i][j]+dp[i-1][j+1])%mod;
				}
				else
				{
					if(a[i]<a[x]||a[i]==a[x]&&i<x)
						dp[i][j]=(dp[i][j]+dp[i-1][j-1])%mod;
					else
						dp[i][j]=(dp[i][j]+dp[i-1][j])%mod;
				}
			}
		}
		for(j=0;j<=sum[n];j++)
			ans=(ans+dp[n][j]*a[x])%mod;
	}
	printf("%lld\n",ans);
}

E1+E2. Abnormal Permutation Pairs

题目大意

设 $p$ 和 $q$ 为 $1$ 到 $n$ 的排列，求有多少对 $(p, q)$ 满足以下条件：

$p$ 的字典序比 $q$ 小；
$p$ 中的逆序数对比 $q$ 中的多；

答案对 $\leq mod \leq 10^9)$ 取模，注意 $m o d$ 不一定为素数。

对于 Easy version， $\leq n \leq 50$ ；
对于 Hard version， $\leq n \leq 500$ ；

题解

若 $p$ 的字典序比 $q$ 小，则对于他们第一个不同的位置 $x$ ，必有 $p_x < q_x$ 。
考虑以 $x$ 为分割点， $p$ 和 $q$ 之间的逆序数对区别。
为便于表述，设 $P [l, r]$ 表示序列 $P$ 的子序列 ${P_l,P_{l+1},...,P_{r}\}$ ， $i n v (P)$ 表示序列 $P$ 的逆序数对数量。

根据之前定义，

$p [1, x - 1] = q [1, x - 1]$ ，即 $i n v (p [1, x - 1]) = i n v (q [1, x - 1])$ 。
若 $\leq i \leq x-1$ ， $\leq j \leq n$ ，则在序列 $p$ 中 $(i, j)$ 是逆序数对的数量，也等于在序列 $q$ 中 $(i, j)$ 是逆序数对的数量。

因此 $i n v (p) > i n v (q)$ ，等价于 $i n v (p [x, n]) > i n v (q [x, n])$ 。
由于逆序数对仅和大小有关，因此 $p [x, n]$ 和 $q [x, n]$ 可以简化为 $1$ 到 $n - x + 1$ 的排列。

Easy Version

设 $f (i, j)$ 表示长度为 $i$ 且有 $j$ 个逆序数对的排列数量。
考虑 $f (i, j)$ 的动态转移方程式。假设该排序首位元素的大小是 $k$ ，那么与 $k$ 相关的逆序数对，共有 $k - 1$ 个。如果删掉 $k$ ，并将所有大于 $k$ 的元素都减 $1$ 则状态变为 $f (i - 1, j - (k - 1))$ 。因此可以得到转移式：
$f(i,j)=\sum_{k=1}^{min(i,j+1)}{f(i-1,j-k+1)}$

整理后可以表示为：
$f(i,j)=\sum_{k=max(j-i+1,0)}^{j}{f(i-1,k)}$

利用前缀和，每次转移可以优化到 $O (1)$ ，即可以在 $O(n^3)$ 复杂度内求得所有 $f (i, j)$ 。

回到 $p < q, i n v (p) > i n v (q)$ 的条件。设第一个不同的位置为 $x$ ，则其等价于 $p_x < q_x,inv(p[x,n])>inv(q[x,n])$ 。
设 $p_x$ 在 $p [x, n]$ 中是第 $k_1$ 小的数， $q_x$ 在 $q [x, n]$ 中是第 $k_2$ 小的数，则在 $p [x, n]$ 中与 $p_x$ 相关的逆序数对有 $k_1-1$ 个，在 $q [x, n]$ 中 $q_x$ 相关的逆序数对有 $k_2-1$ 个。因此条件可以转化为：
$\iff p_x < q_x \iff k_1 < k_2$

$\iff inv(p[x,n]) >inv(q[x,n])$

$\iff k_1+inv(p[x+1,n])>k_2+inv(q[x+1,n])$

其中 $i n v (p [x + 1, n])$ 与 $i n v (q [x + 1, n])$ 可以通过 $f (i, j)$ 求得。

设 $ans_i$ 表示 $n = i$ 时的答案。
计算 $ans_i$ 时，若 $p_1=q_1$ ，那么 $p_1$ 共有 $i$ 种选择，剩余 $i - 1$ 个数共有 $ans_{i-1}$ 种方案。
若 $p_1 \neq q_1$ ，即第一个不同的位置 $x = 1$ ，则需枚举 $p_x,q_x,inv(p[2,n]),inv(q[2,n])$ 四项：

$ans_i=i\cdot ans_{i-1}+\sum_{k_1=1}^{n-x+1}{\sum_{k_2=k_1+1}^{n-x+1}{\sum_{j_1=0}^{\frac{(i-1)(i-2)}{2}}f(i-1,j_1)\cdot \sum_{j_2=0}^{j_1+k_1-k_2-1}{f(i-1,j_2)}}}$

其中 $j_2$ 循环可以用前缀和优化，最终复杂度为 $O(n^5)$ ，如果采用 FFT 加速，可以优化到 $O(n^3 \log(n))$ ，或许也可以通过 Hard Version。

Hard Version

由于我们计算时，实际需要满足的不等式为 $k_1+inv(p[x+1,n])>k_2+inv(q[x+1,n])$ ，其在 $k_1,k_2$ 固定的情况下，成立条件只和 $i n v (p [x + 1, n])$ 与 $i n v (p [x + 1, n])$ 之差有关，考虑以此优化 $d p$ 。

设 $g (i, j)$ 表示有多少对长度为 $i$ 的排列 $p$ 与排列 $q$ ，满足 $i n v (p) - i n v (q) = k$ 。

考虑其转移式，假设长度为 $i - 1$ 的一对排列 $(p^{'}, q^{'})$ ，其逆序数对数量分别为 $j_1$ 和 $j_2$ 。
现在分别向其中插入 $i$ ，使其变为长度为 $i$ 的一堆排列 $(p, q)$ ，假设插入的位置分别在第 $k_1(0 \leq k_1 \leq i-1)$ 个元素与第 $k_2(0 \leq k_2 \leq i-1)$ 个元素之后，则 $inv(p)=j_1+i-1-k_1$ ， $inv(q)=j_2+i-1-k_2$ 。
因此 $inv(p')-inv(q')=j_1-j_2$ 变为了 $inv(p)-inv(q)=j_1-j_2+k_2-k_1$ 。
设 $w=k_2-k_1$ ，当 $w$ 固定时，这样的 $k_2-k_1$ 共有 $n - ∣ w ∣$ 对。
枚举 $w$ 可以得到一个 $O(n^4)$ 的转移式：

$g(n,j)=\sum_{w=-n+1}^{n-1}{(n-|w|) \cdot g(n-1,j+w)}$

$O(n^4)$ 的复杂度依旧超标，再次对其优化：

$\begin{aligned} g(n,j)-g(n,j-1) &=\sum_{w=-n+1}^{n-1}{(n-|w|) \cdot g(n-1,j+w)}-\sum_{w=-n+1}^{n-1}{(n-|w|) \cdot g(n-1,j-1+w)} \\ &=\sum_{w=0}^{n-1}{g(n-1,j+w)}-\sum_{w=-n+1}^{0}{g(n-1,j-1+w))} \end{aligned}$

设前缀和 $s(n,j)=\sum_{k < j}{g(n,k)}$ ，上式可以化简为
$\begin{aligned} g(n,j)&=g(n,j-1) \\&+(s(n-1,j+n-1)-s(n-1,j-1)) \\&-(s(n-1,j-1)-s(n-1,j-n-1)) \end{aligned}$

由于 $j$ 可能是负数，所以还要再加一个偏移量 $B$ 。
这样复杂度就优化到 $O(n^3)$ ，最终答案 $ans_n$ 可以通过以下递推式得到：
$ans_i=i \cdot ans_{i-1}+\sum_{j=1}^i{(i-j)(s(i-1,\frac{i(i+1)}{2})-s(i-1,j))}$

参考代码

#include <bits/stdc++.h>
typedef long long ll;
using namespace std;

const int MAXN=550;
const int B=MAXN*MAXN/2;
ll g[2][B<<1],s[2][B<<1],ans[MAXN];

int main()
{
	ll n,mod;
	int i,j,now,bef,mx;
	scanf("%lld%lld",&n,&mod);
	g[0][B]=s[0][B]=1;
	for(i=B;i<2*B;i++)
		s[0][i]=1;
	for(i=1;i<=n;i++)
	{
		now=i&1;
		bef=now^1;
		memset(g[now],0,sizeof(g[now]));
		memset(s[now],0,sizeof(s[now]));
		mx=i*(i-1)/2;
		g[now][-mx+B]=1;
		for(j=-mx+B+1;j<=mx+B;j++)
		{
			g[now][j]=g[now][j-1]+(s[bef][j+i-1]-s[bef][j-1])-(s[bef][j-1]-s[bef][j-i-1]);
			g[now][j]=(g[now][j]%mod+mod)%mod;
		}
		mx=(i+1)*(i+2)/2;
		for(j=-mx+B;j<=mx+B;j++)
			s[now][j]=(s[now][j-1]+g[now][j])%mod;
		for(j=1;j<i;j++)
			ans[i]=(ans[i]+(s[bef][(i+1)*i/2+B]-s[bef][j+B]+mod)%mod*(i-j))%mod;
		ans[i]=(ans[i]+i*ans[i-1])%mod;
	}
	printf("%lld\n",ans[n]);
}

欣君

关注

13
点赞
踩
2

收藏

觉得还不错? 一键收藏
6
评论
Codeforces Round #729 div.2 A-E题解

视频讲解：TBDA. Odd Set题目大意给定 2n(1≤n≤100)2n(1 \leq n \leq 100)2n(1≤n≤100) 个可能重复的整数 ai(0≤ai≤100)a_i(0 \leq a_i \leq 100)ai(0≤ai≤100) ，问是否能将其分为 nnn 对，使得每对的两个整数之和为奇数。题解只有当每对内包含一个奇数加一个偶数时，才能使其总和为奇数。因此充要条件为奇数个数与偶数个数相等。参考代码#include <bits/stdc++.h>type
复制链接

扫一扫