1. 程式人生 > >【LOJ】#6435. 「PKUSC2018」星際穿越

【LOJ】#6435. 「PKUSC2018」星際穿越

題解

想出70的大眾分之後就棄療了,正解有點神仙

就是首先有個比較顯然的結論,就是要麼是一直往左走,要麼是走一步右邊,然後一直往左走

根據這個可以結合RMQ寫個70分的暴力

我們就考慮,最優的話顯然是走一步左邊就到了目標點,第二步才開始有分叉
假如我們先走了一步左邊,然後就變成了,從\(L[x]\)開始走,下一步可以走到\([L[x],N]\)的所有點最小的轉移點之前,之後再把後來走的點代價都加上1即可
這樣的話,不管是一直走左邊,還是走了一步右邊再走了左邊,情況都被包含了

這個時候考慮這個問題就比較簡單了,可以使用倍增
\(f[i][j]\)表示\([i,n]\)內最小的\(l[x]\)

的值
\(s[i][j]\)表示\(i\)走到\(f[i][j]\)內所有點的距離和

轉移就是
\(f[i][j] = f[f[i][j - 1]][j - 1]\)
\(s[i][j] = s[i][j - 1] + s[f[i][j - 1]][j - 1] + 2^{j - 1} * (f[i][j - 1] - f[i][j])\)

查詢兩端字首和,查的時候直接把\(x\)變成\(L[x]\)進行倍增即可

程式碼

#include <bits/stdc++.h>
#define fi first
#define se second
#define pii pair<int,int>
#define pdi pair<db,int>
#define mp make_pair
#define pb push_back
#define enter putchar('\n')
#define space putchar(' ')
#define eps 1e-8
#define mo 974711
#define MAXN 300005
//#define ivorysi
using namespace std;
typedef long long int64;
typedef double db;
template<class T>
void read(T &res) {
    res = 0;char c = getchar();T f = 1;
    while(c < '0' || c > '9') {
    if(c == '-') f = -1;
    c = getchar();
    }
    while(c >= '0' && c <= '9') {
    res = res * 10 + c - '0';
    c = getchar();
    }
    res *= f;
}
template<class T>
void out(T x) {
    if(x < 0) {x = -x;putchar('-');}
    if(x >= 10) {
    out(x / 10);
    }
    putchar('0' + x % 10);
}
int N,L[MAXN];
int f[MAXN][20];
int64 s[MAXN][20];
void Init() {
    read(N);
    for(int i = 2 ; i <= N ; ++i) read(L[i]);
    f[N][0] = L[N];s[N][0] = N - L[N];
    for(int i = N - 1 ; i >= 1 ; --i) {
    f[i][0] = min(f[i + 1][0],L[i]);s[i][0] = i - f[i][0];
    }
    for(int j = 1 ; j <= 19 ; ++j) {
    for(int i = 1 ; i <= N ; ++i) {
        f[i][j] = f[f[i][j - 1]][j - 1];
        s[i][j] = s[i][j - 1] + s[f[i][j - 1]][j - 1] + 1LL * (f[i][j - 1] - f[i][j]) * (1 << j - 1); 
    }
    }
}
int64 gcd(int64 a,int64 b) {
    return b == 0 ? a : gcd(b,a % b);
}
int64 Calc(int tar,int st) {
    if(tar >= L[st]) return st - tar;
    int64 res = st - L[st];st = L[st];
    int64 sum = 1;
    for(int j = 19 ; j >= 0 ; --j) {
    if(f[st][j] >= tar) {
        res += s[st][j];
        res += 1LL * sum * (st - f[st][j]);
        st = f[st][j];
        sum += 1 << j;
    }
    }
    res += 1LL * (sum + 1) * (st - tar);
    return res;
}
void Solve() {
    int Q;int l,r,x;
    read(Q);
    while(Q--) {
    read(l);read(r);read(x);
    int64 u = Calc(l,x) - Calc(r + 1,x),d = r - l + 1,g = gcd(u,d);
    u /= g;d /= g;
    out(u);putchar('/');out(d);enter;
    }
}
int main() {
#ifdef ivorysi
    freopen("f1.in","r",stdin);
#endif
    Init();
    Solve();
}