BZOJ3672: [Noi2014]購票(dp 斜率優化 點分治 二分 凸包)

来源:https://www.cnblogs.com/zwfymqz/archive/2019/01/02/10211469.html
-Advertisement-
Play Games

題意 "題目鏈接" Sol 介紹一種神奇的點分治的做法 啥?這都有根樹了怎麼點分治?? 嘿嘿,這道題的點分治不同於一般的點分治。正常的點分治思路大概是先統計過重心的,再遞歸下去 實際上一般的點分治與統計順序關係不大,也就是說我可以先統計再遞歸,或者先遞歸再統計。 但是這題不單單是統計,它是dp,存在 ...


題意

題目鏈接

Sol

介紹一種神奇的點分治的做法

啥?這都有根樹了怎麼點分治??

嘿嘿,這道題的點分治不同於一般的點分治。正常的點分治思路大概是先統計過重心的,再遞歸下去

實際上一般的點分治與統計順序關係不大,也就是說我可以先統計再遞歸,或者先遞歸再統計。

但是這題不單單是統計,它是dp,存在決策順序問題,我們就需要換一種思路了。

首先我們可以這樣考慮:對於每個點\(x\),找出子樹重心\(root\)對除去重心外的部分遞歸執行該操作,那麼回溯回來的時候,我們預設除重心的子樹外答案都已經更新好了。

接下來考慮重心子樹內的點的轉移,我們只需要考慮從\(root\)\(x\)的路徑,顯然排序之後雙指針可以做到\(nlogn\)的複雜度。(對轉移位置按深度排序,對要更新的點按深度 - 限制長度排序,雙指針的時候維護一下凸包,因為\(p\)不單調所以需要在凸包上二分)

複雜度不太會嚴格的證明,但是跑的飛快。因為雖然我們的分治結構變了,但每次還是找重心更新答案,所以複雜度是有保證的。

#include<bits/stdc++.h> 
#define int long long 
#define LL long long 
using namespace std;
const int MAXN = 1e6 + 10, mod = 1e9 + 7, INF = 3e18 + 10;
const double eps = 1e-9;
template <typename A, typename B> inline bool chmin(A &a, B b){if(a > b) {a = b; return 1;} return 0;}
template <typename A, typename B> inline bool chmax(A &a, B b){if(a < b) {a = b; return 1;} return 0;}
template <typename A> inline void debug(A a){cout << a << '\n';}
inline int read() {
    char c = getchar(); int x = 0, f = 1;
    while(c < '0' || c > '9') {if(c == '-') f = -1; c = getchar();}
    while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
    return x * f;
}
int N, T, fa[MAXN], d[MAXN], pp[MAXN], qq[MAXN], lim[MAXN], siz[MAXN], Sum, vis[MAXN], Mx, root, f[MAXN], st[MAXN], top, cnt;
vector<int> v[MAXN];
void dfs(int x) {
    d[x] += d[fa[x]]; siz[x] = 1;
    for(int i = 0; i < v[x].size(); i++) {
        int to = v[x][i]; if(to == fa[x]) continue;
        dfs(to); siz[x] += siz[to];
    }
}
void Find(int x) {
    //printf("%d\n", x);
    int mx = 0; siz[x] = 1;
    for(int i = 0; i < v[x].size(); i++) {
        int to = v[x][i]; if(to == fa[x] || vis[to]) continue;
        Find(to); siz[x] += siz[to];
        chmax(mx, siz[to]); 
    }
    chmax(mx, Sum - siz[x]);
    if(mx < Mx && siz[x] > 1) Mx = mx, root = x;
}
struct Node {
    int id, dis;
    bool operator < (const Node &rhs) const {
        return dis > rhs.dis;
    }
}a[MAXN];
void dfs2(int x) {
    a[++cnt].id = x;
    a[cnt] = (Node) {x, d[x] - lim[x]};
    for(int i = 0; i < v[x].size(); i++) {
        int to = v[x][i]; if(to == fa[x] || vis[to]) continue;
        dfs2(to);
    }
}
int Y(int x) {
    return f[x];
}
int X(int x) {
    return d[x];
}
double slope(int x, int y) {
    return (double) (Y(y) - Y(x)) / (X(y) - X(x));
} 
void insert(int x) {
    while(top > 1 && slope(st[top], x) > slope(st[top - 1], st[top])) top--;
    st[++top] = x;
}
int Search(double x, int id) {
    if(!top) return INF;
    int l = 1, r = top - 1, ans = 1;
    while(l <= r) {
        int mid = l + r >> 1;
        if((slope(st[mid], st[mid + 1]) >= x)) ans = mid + 1, l = mid + 1;
        else r = mid - 1;
    }
    return f[st[ans]] - d[st[ans]] * pp[id];
}
void solve(int x, int tot, int up) {
    vector<int> pot;
    for(int t = x; t != up; t = fa[t]) 
        pot.push_back(t);
    cnt = 0; 
    for(int i = 0; i < v[x].size(); i++) {
        int to = v[x][i]; if(to == fa[x]) continue;
        dfs2(to);//dep´Ó´óµ½´ïС£¬dis´Ó´óµ½Ð¡
    }
    sort(a + 1, a + cnt + 1);
    int now = 0; top = 0;
    for(int i = 1; i <= cnt; i++) {
        int cur = a[i].id;
        while(now <= pot.size() - 1 && d[pot[now]] >= a[i].dis) insert(pot[now++]);
        chmin(f[cur], Search(pp[cur], cur) + qq[cur] + d[cur] * pp[cur]);
    }
}
void work(int x, int tot) {
    if(tot == 1) return ;
    root = x; Sum = tot; Mx = Sum; Find(root);
    int rt = root;
    for(int i = 0; i < v[rt].size(); i++) {
        int to = v[rt][i]; if(to == fa[rt]) continue;
        vis[to] = 1;
    }
    work(x, tot - siz[root] + 1);
    solve(rt, tot, fa[x]);
    for(int i = 0; i < v[rt].size(); i++) {
        int to = v[rt][i]; if(to == fa[rt]) continue;
        work(to, siz[to]);
    }   
}
signed main() {
    //freopen("a.in", "r", stdin);
    N = read(); T = read();
    for(int i = 2; i <= N; i++) {
        fa[i] = read();
        v[fa[i]].push_back(i); v[i].push_back(fa[i]);
        d[i] = read(); pp[i] = read(); qq[i] = read();  lim[i] = read();
    }
    dfs(1);
    memset(f, 0x7f7f7f, sizeof(f)); f[1] = 0;
    work(1, N);
    for(int i = 2; i <= N; i++) cout << f[i] << '\n';
    return 0;
}
/*
7 3
1 2 20 0 3
1 5 10 100 5
2 4 10 10 10
2 9 1 100 10
3 5 20 100 10
4 4 20 0 10
*/

您的分享是我們最大的動力!

-Advertisement-
Play Games
更多相關文章
  • 可迭代對象 以直接作用於 for 迴圈的數據類型有以下幾種: 一類是集合數據類型,如 list 、 tuple 、 dict 、 set 、 str 等; 一類是 generator ,包括生成器和帶 yield 的generator function。 這些可以直接作用於 for 迴圈的對象統稱為 ...
  • 運行結果: 結論:非守護線程結束之後jvm退出,守護線程沒有要守護的線程,也進行退出 守護線程和用戶線程唯一的區別在於: 補充說明: 定義:守護線程--也稱“服務線程”,在沒有用戶線程可服務時會自動離開。 優先順序:守護線程的優先順序比較低,用於為系統中的其它對象和線程提供服務。 設置:通過setDae ...
  • 由來: Java是由Sun Microsystems公司於1995年5月推出的Java面向對象程式設計語言和Java平臺的總稱。由James Gosling和同事們共同研發,併在1995年正式推出。 特性: 1.簡單性 :Java吸收了C++語言的各種優點,使得大多數程式員很容易學習和使用。另一方面 ...
  • 1.什麼是進程、線程 進程可以簡單的理解為應用程式。多線程,相當於多個人共同完成一件事情,每個線程,獨立完成一件任務。因此,一個進程至少有一個線程,一個線程不能獨立存在,它必須是進程的一部分。 進程:每個獨立運行著的程式。具有獨立的記憶體空間和系統資源。(建大廈) 線程:是一個進程內部的一條執行路徑。 ...
  • 在軟體開發領域中,人們經常會用到這一個概念——“設計模式”(design pattern),它是一種針對軟體設計的共性問題而提出的解決方案。在一本聖經級的書籍《設計模式:可復用面向對象軟體的基礎》(1991年,Design Patterns - Elements of Reusable Object ...
  • `synchronized java 同步`的關鍵字,其典型的作用域如下所示. 1 對象鎖 1.1 代碼塊修飾(對象) 此時, 用於保證 函數中的被 大括弧包裹的代碼同步執行. 作用的對象為 的對象實例,例如 函數中的 以及 . Tips: 1. 若在多個線程中被調用,其輸出順序將保證同步,按照 的 ...
  • Django 為大家提供了一個完善的後臺管理系統—admin,但是這個後臺管理系統總體來說不太適合國人的習慣,所以有大神就使用 bootstrap 和 jQuery,為我們開發了一個第三 方的 Django 後臺管理系統—xadmin,這個界面的友好度比較高,目前企業的使用頻率比較高,我們來學習一些 ...
  • 標準輸入對話框(QInputDialog)系列: 主要模塊 效果如下: 完整代碼: 打開文件(QFileDialog)、顏色(QColorDialog)、字體(QFontDialog)對話框 效果如下: 參考: https://zhuanlan.zhihu.com/p/29321561 完整代碼: ...
一周排行
    -Advertisement-
    Play Games
  • 移動開發(一):使用.NET MAUI開發第一個安卓APP 對於工作多年的C#程式員來說,近來想嘗試開發一款安卓APP,考慮了很久最終選擇使用.NET MAUI這個微軟官方的框架來嘗試體驗開發安卓APP,畢竟是使用Visual Studio開發工具,使用起來也比較的順手,結合微軟官方的教程進行了安卓 ...
  • 前言 QuestPDF 是一個開源 .NET 庫,用於生成 PDF 文檔。使用了C# Fluent API方式可簡化開發、減少錯誤並提高工作效率。利用它可以輕鬆生成 PDF 報告、發票、導出文件等。 項目介紹 QuestPDF 是一個革命性的開源 .NET 庫,它徹底改變了我們生成 PDF 文檔的方 ...
  • 項目地址 項目後端地址: https://github.com/ZyPLJ/ZYTteeHole 項目前端頁面地址: ZyPLJ/TreeHoleVue (github.com) https://github.com/ZyPLJ/TreeHoleVue 目前項目測試訪問地址: http://tree ...
  • 話不多說,直接開乾 一.下載 1.官方鏈接下載: https://www.microsoft.com/zh-cn/sql-server/sql-server-downloads 2.在下載目錄中找到下麵這個小的安裝包 SQL2022-SSEI-Dev.exe,運行開始下載SQL server; 二. ...
  • 前言 隨著物聯網(IoT)技術的迅猛發展,MQTT(消息隊列遙測傳輸)協議憑藉其輕量級和高效性,已成為眾多物聯網應用的首選通信標準。 MQTTnet 作為一個高性能的 .NET 開源庫,為 .NET 平臺上的 MQTT 客戶端與伺服器開發提供了強大的支持。 本文將全面介紹 MQTTnet 的核心功能 ...
  • Serilog支持多種接收器用於日誌存儲,增強器用於添加屬性,LogContext管理動態屬性,支持多種輸出格式包括純文本、JSON及ExpressionTemplate。還提供了自定義格式化選項,適用於不同需求。 ...
  • 目錄簡介獲取 HTML 文檔解析 HTML 文檔測試參考文章 簡介 動態內容網站使用 JavaScript 腳本動態檢索和渲染數據,爬取信息時需要模擬瀏覽器行為,否則獲取到的源碼基本是空的。 本文使用的爬取步驟如下: 使用 Selenium 獲取渲染後的 HTML 文檔 使用 HtmlAgility ...
  • 1.前言 什麼是熱更新 游戲或者軟體更新時,無需重新下載客戶端進行安裝,而是在應用程式啟動的情況下,在內部進行資源或者代碼更新 Unity目前常用熱更新解決方案 HybridCLR,Xlua,ILRuntime等 Unity目前常用資源管理解決方案 AssetBundles,Addressable, ...
  • 本文章主要是在C# ASP.NET Core Web API框架實現向手機發送驗證碼簡訊功能。這裡我選擇是一個互億無線簡訊驗證碼平臺,其實像阿裡雲,騰訊雲上面也可以。 首先我們先去 互億無線 https://www.ihuyi.com/api/sms.html 去註冊一個賬號 註冊完成賬號後,它會送 ...
  • 通過以下方式可以高效,並保證數據同步的可靠性 1.API設計 使用RESTful設計,確保API端點明確,並使用適當的HTTP方法(如POST用於創建,PUT用於更新)。 設計清晰的請求和響應模型,以確保客戶端能夠理解預期格式。 2.數據驗證 在伺服器端進行嚴格的數據驗證,確保接收到的數據符合預期格 ...