NoiminのNoise

競技プログラミング (多め) とWeb (たまに) ,自然言語処理 (ブログではまだ)。数式の書き方を一気に KaTeX に変えようとして記事を全削除してインポートし直すなどしたので,過去にブックマークされた記事は URL が変わってしまっている可能性があります…….

AtCoder Beginner Contest 128 F - Frog Jump

問題: F - Frog Jump

公式解説: https://img.atcoder.jp/abc128/editorial.pdf

問題概要

座標 {\displaystyle 0, 1,  2, \dots, N-2, N-1 (3 \leq N \leq 10^5) } 上を移動する.移動は座標 0 から始めて,{\displaystyle A (\gt 0) } だけ正の方向に進むことと,{\displaystyle B (\gt 0) } だけ負の方向に進むことを交互に繰り返す. {\displaystyle N-1 } にたどり着いた地点で終了する.同じ座標を2度訪れてはいけない.

座標 0 と 座標 {\displaystyle N-1 } 以外の各座標{\displaystyle i }は,訪れると点数{\displaystyle s_i }がもらえる.負の点数もありうる.

{\displaystyle A }{\displaystyle B } を適切に定め,もらえる点数の合計を最大化せよ.

解法概要

座標 {\displaystyle N } 以降に移動することは考えなくて良いので,{\displaystyle N-1 } にたどり着くことができるとしたら正の方向に進んだときだけである.また,{\displaystyle A \leq B } であるような場合も考えなくて良い. したがって,座標の遷移は次のような形になるはずである:

{\displaystyle 0, A, A-B, 2A-B, 2A-2B, \dots, xA-(x-1)B, xA-xB,  (x+1)A-xB=N-1 (2x+1 }回移動するとする) .

これの奇数 / 偶数番目を抜き出した数列を考えると,隣接する要素の差が{\displaystyle A-B } の等差数列となっている:

奇数番目の要素を抜き出した数列: {\displaystyle 0, A-B, 2A-2B, \dots, xA-xB }

偶数番目の要素を抜き出した数列{\displaystyle A, 2A-B, \dots, xA-(x-1)B, (x+1)A-xB=N-1 (2x+1 }

ここで,{\displaystyle A-B=C} とおくと

奇数番目の要素を抜き出した数列: {\displaystyle 0, C, 2C, \dots, xC }

偶数番目の要素を抜き出した数列{\displaystyle A, A+C, \dots, A+(x-1)C, A+xC=N-1}

だいぶスッキリして見える (主観) .偶数番目の要素について,(0 との差ではなく) N-1 との差に注目して式変形をすると:

奇数番目の要素を抜き出した数列: {\displaystyle 0, C, 2C, \dots, xC }

偶数番目の要素を抜き出した数列{\displaystyle N-1-xC, N-1-(x-1)C, \dots, N-1-C, N-1}

したがって,題意の通りに移動したときのスコアは {\displaystyle C,x } の関数{\displaystyle f(C,x) } として表現することができる.

{\displaystyle f(C,x) = f(C, x-1) + s_{xC} + s_{N-1-xC} } と表現できることから,{\displaystyle x, C } が決まればスコアは程数時間で求めることができる.

また,{\displaystyle xC \lt N } でなければならないことから,全体の時間計算量は {\displaystyle O(N \log N) } である.*1

ただし,{\displaystyle C }{\displaystyle N-1 } の約数であるような場合のみ注意が必要である.このような場合は,{\displaystyle xC \leq N-1-xC } となったとき,同じ座標のスコアを 2 回足してしまうことになる.同じ座標を 2 度訪れてはいけないから,{\displaystyle C }{\displaystyle N-1 } の約数である場合には {\displaystyle xC \leq N-1-xC } となった地点でループを抜けなければならない.

ソースコード

#include <iostream>
#include <vector>
#include <algorithm>

using namespace std;

typedef long long ll;

int main() {
    int n;
    cin >> n;
    vector<ll> s(n);
    for(int i=0;i<n;++i){
        cin >> s[i];
    }

    ll score = 0LL;
    for(int c=1;c<n/2;++c) {
        ll tmp_score = 0LL;
        if((n-1) % c) {
            for(int x=0;n-1-c*x>=c;++x) {
                tmp_score += s[c*x] + s[n-1-c*x];
                score = max(score, tmp_score);
            }
        } else {
            for(int x=0;c*x<n-1-c*x;++x) {
                tmp_score += s[c*x] + s[n-1-c*x];
                score = max(score, tmp_score);
            }
        }
        
    }

    cout << score << endl;
}

所感

考察難易度 >> 実装難易度な問題は不得意.

*1:{\displaystyle \log N = \sum_{i=1}^{\infty} \frac{1}{i} } より