CF1174E (Ehab and the Expected GCD Problem) 题解——数论与dp的完美结合

Description

Solution

Lemma 1

$g_{i+1}=t{g}_i(t\in N^{\ast })$，且$t\le 3$。

Certification 1

根据$\gcd$的定义，不难发现$g_{i+1}$一定是$g_i$的约数。

例如，一个合法的$g$序列是 { 12 , 6 , 6 , 6 , 6 , 3 , 2 , 2 , 1 , 1 , 1 , 1 , 1 , 1 , 1 , 1 } \{ 12,6,6,6,6,3,2,2,1,1,1,1,1,1,1,1\} {12,6,6,6,6,3,2,2,1,1,1,1,1,1,1,1}。如果我们将这个序列给缩一下，那么会变为 g ′ = { 12 , 6 , 3 , 2 , 1 } g'=\{12,6,3,2,1\} g′={12,6,3,2,1}。

此时$f(p)$就是$g^{\prime }$的长度。当$f(p)$最大时，我们肯定要贪心地让$g_{i+1}^{\prime }$比$g_i^{\prime }$小的倍数越少越好，否则$g^{\prime }$的长度必然被削减。

令这个小的倍数为$t^{\prime }$，不难发现$t^{\prime }>2$。考虑当$t^{\prime }=4$时，我们显然可以在其中加一项——例如 g ′ = { 24 , 12 , 3 , 1 } g'=\{24,12,3,1\} g′={24,12,3,1}可以优化为 g ′ = { 24 , 12 , 6 , 3 , 1 } g'=\{24,12,6,3,1\} g′={24,12,6,3,1}，我们在其中添了一项$6$，这样不仅没有影响$g^{\prime }$的合法性也使得$g^{\prime }$的长度变大。

扩展得，当$4\le t^{\prime }$时$g^{\prime }$一定不是最优。所以Lemma 1中$t\le 3$。

证毕。

Lemma 2

$g_{i+1}=t{g}_i(t\in N^{\ast })$，且$t\le 3$；更进一步的，$t=3$的次数最多只有1次。

Certification 1

考虑到$2^3\le 3^2$，所以连续两次$t=3$肯定差于三次$t=2$，所以$t=3$的次数最多只有$1$次。

---

考虑$dp$。

状态设计: $d{p}_{i,j,k}$: 看到第$i$个数，目前$\gcd$为$2^j{3}^k$。

状态转移:
①$\gcd$不变。$$d{p}_{i,j,k}=(\lfloor \frac{n}{p}\rfloor -i+1)d{p}_{i-1,j,k}$$
②$\gcd$除以$2$($j+1\le \lfloor \log n\rfloor$): $$d{p}_{i,j,k}=(\lfloor \frac{n}{2^j{3}^k}\rfloor -\lfloor \frac{n}{2^{j+1}{3}^k}\rfloor )d{p}_{i,j+1,k}$$
③$\gcd$除以$3$($k=0$): $$d{p}_{i,j,k}=(\lfloor \frac{n}{2^j{3}^k}\rfloor -\lfloor \frac{n}{2^j{3}^{k+1}}\rfloor )d{p}_{i,j,k+1}$$

为什么③的合法转移要求$k=0$呢？因为，对于一次③的转移相当于使$t$成为了一次$3$。根据之前的引理，$t=3$的次数只能有$1$次。

边界:
①$d{p}_{1,\lfloor \log n\rfloor ,0}=1$
②$d{p}_{1,\lfloor \log n\rfloor -1,1}=1(2^{\lfloor \log n\rfloor -1}\times 3\le n)$

答案: $d{p}_{n,0,0}$。

总时间复杂度$O(n{\log }_{2}n)$。

---

Summary

本题考查了性质挖掘能力与基本的$dp$能力，可以说是一道不可多得的好题。

我在做这题的时候，性质全部挖掘了出来，可是$dp$的状态设计却并不是最优，这导致我列出了一个关于下降幂的状态转移，从而时间复杂度多出了一个$n$。所以说，我是不是很菜？

这还用问……

Code

①在实现的时候，你必须采用滚动数组，否则空间开不下。
②一些卡常技巧: 少取模，开$O2$。

#include <bits/stdc++.h>
#define int long long
using namespace std;
const int g1=25,g2=25,mod=1e9+7;

int n,L2,L3;
int dp[2][g1][2],bin[g1],sea[g2];

int get_log_two(int tmp){
	int p=1,cnt=0;
	while (p<=tmp)  p*=2,cnt++;
	return cnt-1;
}

int get_log_three(int tmp){
	int p=1,cnt=0;
	while (p<=tmp)  p*=3,cnt++;
	return cnt-1;
}

int add(int tmpx,int tmpy){return (tmpx+(max(tmpy,0ll)))%mod;}
void init(){
	bin[0]=1,sea[0]=1;
	for (int i=1;i<=20;i++)  bin[i]=(bin[i-1]*2)%mod;
	for (int i=1;i<=13;i++)  sea[i]=(sea[i-1]*3)%mod;
}

signed main(){
	cin>>n;init();
	L2=get_log_two(n),L3=get_log_three(n);
	dp[0][L2][0]=1;
	if (((1ll<<L2)/2)*3<=n)  dp[0][L2-1][1]=1;
	for (int i=2;i<=n;i++){
		for (int j=0;j<=L2;j++){
			for (int k=0;k<2;k++){
				int p=bin[j]*sea[k];
				if (p>n)  continue;
				
				dp[1][j][k]=add(dp[1][j][k],((n/p)-i+1)*dp[0][j][k]);
				if (j<L2)  dp[1][j][k]=add(dp[1][j][k],((n/p)-(n/(2*p)))*dp[0][j+1][k]);
				if (k==0)  dp[1][j][k]=add(dp[1][j][k],((n/p)-(n/(3*p)))*dp[0][j][k+1]);
			}
		}
		for (int j=0;j<=L2;j++){
			for (int k=0;k<2;k++)  dp[0][j][k]=dp[1][j][k],dp[1][j][k]=0;
		}
	}
	cout<<dp[0][0][0]<<endl;
	
	return 0;
}