这道题,很容易想到 dp。
令 f i , j f_{i,j} fi,j 表示第 i i i 行第 j j j 列的最大 gcd,则状态转移方程为 f i , j = max ( gcd ( f i − 1 , j , a i , j ) , ( gcd ( f i , j − 1 , a i , j ) ) f_{i,j}=\max(\gcd(f_{i-1,j},a_{i,j}),(\gcd(f_{i,j-1},a_{i,j})) fi,j=max(gcd(fi−1,j,ai,j),(gcd(fi,j−1,ai,j))。
然而这明显不对。
样例 2 2 2 中,若按这种方法去做,答案为 2 2 2,而实际上,是 3 3 3。
也就是一个状态最优解不一定为下一个状态的最优解!
可怜的蒟蒻我不知如何是好,只能硬着头皮优化。
把二维数组 f f f 升级成三维。我们把每种情况的 gcd 都存进 f i , j f_{i,j} fi,j 里,状态转移方程变为 f i , j . p u s h b a c k ( gcd ( f i , j − 1 , k , a i , j ) f_{i,j}.pushback(\gcd(f_{i,j-1,k},a_{i,j}) fi,j.pushback(gcd(fi,j−1,k,ai,j) 与 f i , j . p u s h b a c k ( gcd ( f i − 1 , j , k , a i , j ) f_{i,j}.pushback(\gcd(f_{i-1,j,k},a_{i,j}) fi,j.pushback(gcd(fi−1,j,k,ai,j)。
酱紫应该能行了吧!我小心翼翼地交上去,结果……
果然不出我所料,MLE ON #12!
我只好用 set 代替 vector。
硬着头皮再交一发!
TLE ON #31。
我只得继续优化。
考虑以下性质: gcd ( n , k ) ≥ gcd ( m , k ) ( n = m q ) ( q ∈ N + ) \gcd(n,k) \ge \gcd(m,k) (n=mq)(q∈{N_+}) gcd(n,k)≥gcd(m,k)(n=mq)(q∈N+)。
于是,我们可以删掉每个 f i , j , k f_{i,j,k} fi,j,k 的真因子,只留下它自己。
先初始化出每一个
k
k
k 的真因子
(
1
≤
k
≤
1
0
6
)
(1 \le k \le 10^6)
(1≤k≤106)
,在每次计算完
f
i
,
j
,
k
f_{i,j,k}
fi,j,k 后删掉真因子。
预处理时间复杂度为 O ( m × m ) O(m \times \sqrt m) O(m×m) 为 1 e 9 1e9 1e9,求值与优化时间复杂度为 O ( n 2 × k ) O(n^2 \times k) O(n2×k),3秒时限绰绰有余。
Accepted!