【LOJ】#6435. 「PKUSC2018」星際穿越
阿新 • • 發佈:2018-12-15
題解
想出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(); }