題目

某國為了防御敵國的導(dǎo)彈襲擊，發(fā)展出一種導(dǎo)彈攔截系統(tǒng)。但是這種導(dǎo)彈攔截系統(tǒng)有一個缺陷：雖然它的第一發(fā)炮彈能夠到達任意的高度、并且能夠攔截任意速度的導(dǎo)彈，但是以后每一發(fā)炮彈都不能高于前一發(fā)的高度，其攔截的導(dǎo)彈的飛行速度也不能大于前一發(fā)。某天，雷達捕捉到敵國的導(dǎo)彈來襲。由于該系統(tǒng)還在試用階段，所以只有一套系統(tǒng)，因此有可能不能攔截所有的導(dǎo)彈。
在不能攔截所有的導(dǎo)彈的情況下，我們當然要選擇使國家損失最小、也就是攔截導(dǎo)彈的數(shù)量最多的方案。但是攔截導(dǎo)彈數(shù)量的最多的方案有可能有多個，如果有多個最優(yōu)方案，那么我們會隨機選取一個作為最終的攔截導(dǎo)彈行動藍圖。
我方間諜已經(jīng)獲取了所有敵軍導(dǎo)彈的高度和速度，你的任務(wù)是計算出在執(zhí)行上述決策時，每枚導(dǎo)彈被攔截掉的概率。

輸入格式

第一行包含一個正整數(shù)n，表示敵軍導(dǎo)彈數(shù)量；
下面 行按順序給出了敵軍所有導(dǎo)彈信息：
第i+1行包含2個正整數(shù)hi和vi，分別表示第 枚導(dǎo)彈的高度和速度。

輸出格式

輸出包含兩行。
第一行為一個正整數(shù)，表示最多能攔截掉的導(dǎo)彈數(shù)量；
第二行包含n個0到1之間的實數(shù)，第i個數(shù)字表示第i枚導(dǎo)彈被攔截掉的概率（你可以保留任意多位有效數(shù)字）。

輸入樣例

4

3 30

4 40

6 60

3 30

輸出樣例

2

0.33333 0.33333 0.33333 1.00000

提示

對于100%的數(shù)據(jù)，1≤n≤5*104， 1≤hi ，vi≤109；

均勻分布著約30%的數(shù)據(jù)，所有vi均相等。

均勻分布著約50%的數(shù)據(jù)，滿足1≤hi ，vi≤1000。

題解

二維LIS
考慮cdq分治，套上樹狀數(shù)組可以得到答案

但是要算概率就有些麻煩了
先要算出總方案數(shù)，計算過程中記錄，最后最大值的地方方案數(shù)之和就是總方案數(shù)

至于每個點有多少方案經(jīng)過
我們反著再做一次LIS，這樣一個點往前往后之和 - 1如果等于答案，那么就將往前往后方案數(shù)乘起來就是總的方案數(shù)

碼著真tm累

#include<iostream> #include<cstdio> #include<cmath> #include<map> #include<cstring> #include<algorithm> #define LL long long int #define Redge(u) for (int k = h[u],to; k; k = ed[k].nxt) #define REP(i,n) for (int i = 1; i <= (n); i++) #define BUG(s,n) for (int i = 1; i <= (n); i++) cout<<s[i]<<' '; puts(""); #define lbt(x) (x & -x) #define mp(a,b) make_pair<double,double>(a,b) #define cp pair<double,double> using namespace std; const int maxn = 50005,maxm = 100005,INF = 1000000000; inline int read(){int out = 0,flag = 1; char c = getchar();while (c < 48 || c > 57){if (c == '-') flag = -1; c = getchar();}while (c >= 48 && c <= 57){out = (out << 3) + (out << 1) + c - 48; c = getchar();}return out * flag; } struct node{int t,x,y; double f[2],g[2];}e[maxn],t[maxn]; inline bool operator <(const node& a,const node& b){return a.x == b.x ? a.y < b.y : a.x < b.x; } inline bool cmp(const node& a,const node& b){return a.t < b.t; } int b[maxn],tot,c[maxn],tt; int getx(int x){return lower_bound(b + 1,b + 1 + tot,x) - b;} int gety(int x){return lower_bound(c + 1,c + 1 + tt,x) - c;} int n; double mx[maxn],sum[maxn]; void upd(int u,double v,double s){while (u <= tt){if (v > mx[u]){mx[u] = v;sum[u] = s;}else if (v == mx[u])sum[u] += s;u += lbt(u);} } cp query(int u){cp re = mp(0,0);while (u){if (mx[u] > re.first){re.first = mx[u];re.second = sum[u];}else if (mx[u] == re.first)re.second += sum[u];u -= lbt(u);}return re; } void cls(int u){while (u <= tt){sum[u] = mx[u] = 0;u += lbt(u);} } void init(){n = read();for (int i = 1; i <= n; i++){e[i].t = i;e[i].x = b[i] = read();e[i].y = c[i] = read();}sort(b + 1,b + 1 + n);sort(c + 1,c + 1 + n);tot = tt = 1;for (int i = 2; i <= n; i++) if (b[i] != b[tot]) b[++tot] = b[i];for (int i = 2; i <= n; i++) if (c[i] != c[tt]) c[++tt] = c[i];for (int i = 1; i <= n; i++){e[i].x = getx(e[i].x);e[i].y = gety(e[i].y);} } void cdq(int l,int r,int p){if (l == r){if (!e[l].f[p]) e[l].f[p] = e[l].g[p] = 1;return;}int mid = l + r >> 1,li = l,ri = mid + 1;for (int i = l; i <= r; i++){if (e[i].t <= mid) t[li++] = e[i];else t[ri++] = e[i];}for (int i = l; i <= r; i++) e[i] = t[i];cdq(l,mid,p);sort(e + l,e + mid + 1);cp tmp; li = l; ri = mid + 1;while (li <= mid && ri <= r){if (e[li].x <= e[ri].x) upd(e[li].y,e[li].f[p],e[li].g[p]),li++;else {tmp = query(e[ri].y);if (!tmp.first) {ri++; continue;}if (tmp.first + 1 > e[ri].f[p]){e[ri].f[p] = tmp.first + 1;e[ri].g[p] = tmp.second;}else if (tmp.first + 1 == e[ri].f[p]){e[ri].g[p] += tmp.second;}ri++;}}while (ri <= r){tmp = query(e[ri].y);if (!tmp.first) {ri++; continue;}if (tmp.first + 1 > e[ri].f[p]){e[ri].f[p] = tmp.first + 1;e[ri].g[p] = tmp.second;}else if (tmp.first + 1 == e[ri].f[p]){e[ri].g[p] += tmp.second;}ri++;}for (int i = l; i < li; i++) cls(e[i].y);cdq(mid + 1,r,p); } void solve(){for (int i = 1; i <= n; i++){e[i].x = tot - e[i].x + 1;e[i].y = tt - e[i].y + 1;}sort(e + 1,e + 1 + n);cdq(1,n,0);for (int i = 1; i <= n; i++){e[i].x = tot - e[i].x + 1;e[i].y = tt - e[i].y + 1;e[i].t = n - e[i].t + 1;}sort(e + 1,e + 1 + n);cdq(1,n,1);for (int i = 1; i <= n; i++) e[i].t = n - e[i].t + 1;sort(e + 1,e + 1 + n,cmp);double ans = 0,sum = 0;for (int i = 1; i <= n; i++){if (e[i].f[0] > ans){ans = e[i].f[0];sum = e[i].g[0];}else if (e[i].f[0] == ans){sum += e[i].g[0];}}printf("%.lf\n",ans);for (int i = 1; i <= n; i++){if (e[i].f[0] + e[i].f[1] - 1 < ans) printf("0 ");else printf("%.6lf ",e[i].g[0] * e[i].g[1] / sum);} } int main(){init();solve();return 0; }

轉(zhuǎn)載于:https://www.cnblogs.com/Mychael/p/8580956.html

總結(jié)

以上是生活随笔為你收集整理的BZOJ2244 [SDOI2011]拦截导弹 【cdq分治 + 树状数组】的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。