静态区间第K大值。主席树和划分树都可解。

 /* 3727 */
 #include <iostream>
 #include <sstream>
 #include <string>
 #include <map>
 #include <queue>
 #include <set>
 #include <stack>
 #include <vector>
 #include <deque>
 #include <algorithm>
 #include <cstdio>
 #include <cmath>
 #include <ctime>
 #include <cstring>
 #include <climits>
 #include <cctype>
 #include <cassert>
 #include <functional>
 #include <iterator>
 #include <iomanip>
 using namespace std;
 //#pragma comment(linker,"/STACK:102400000,1024000")

 #define sti                set<int>
 #define stpii            set<pair<int, int> >
 #define mpii            map<int,int>
 #define vi                vector<int>
 #define pii                pair<int,int>
 #define vpii            vector<pair<int,int> >
 #define rep(i, a, n)     for (int i=a;i<n;++i)
 #define per(i, a, n)     for (int i=n-1;i>=a;--i)
 #define clr                clear
 #define pb                 push_back
 #define mp                 make_pair
 #define fir                first
 #define sec                second
 #define all(x)             (x).begin(),(x).end()
 #define SZ(x)             ((int)(x).size())
 // #define lson            l, mid, rt<<1
 // #define rson            mid+1, r, rt<<1|1

 typedef struct {
     int q;
     int l, r, k;
 } Ques_t;

 ;
 ;
 ;
 int a[maxn], m;
 Ques_t Q[maxq];
 int T[maxn];
 int lson[maxm], rson[maxm], c[maxm];
 int tot, cur, q;

 void init() {
     tot = cur = m = ;
 }

 int Build(int l, int r) {
     int rt = tot++;

     c[rt] = ;
     if (l == r)
         return rt;

     ;

     lson[rt] = Build(l, mid);
     rson[rt] = Build(mid+, r);

     return rt;
 }

 int Insert(int rt, int x, int delta) {
     int nrt = tot++, ret = nrt;
     , r = m - , mid;

     c[nrt] = c[rt] + delta;
     while (l < r) {
         mid = (l + r) >> ;
         if (x <= mid) {
             lson[nrt] = tot++;
             rson[nrt] = rson[rt];
             nrt = lson[nrt];
             rt = lson[rt];
             r = mid;
         } else {
             lson[nrt] = lson[rt];
             rson[nrt] = tot++;
             nrt = rson[nrt];
             rt = rson[rt];
             l = mid + ;
         }
         c[nrt] = c[rt] + delta;
     }

     return ret;
 }

 int Query_kth(int lrt, int rrt, int k) {
     , r = m - , mid;
     int tmp;

     while (l < r) {
         mid = (l + r) >> ;
         tmp = c[lson[rrt]] - c[lson[lrt]];
         if (tmp >= k) {
             lrt = lson[lrt];
             rrt = lson[rrt];
             r = mid;
         } else {
             k -= tmp;
             lrt = rson[lrt];
             rrt = rson[rrt];
             l = mid + ;
         }
     }

     return l;
 }

 int Query_Rank(int rt, int L, int R, int l, int r) {
     if (L==l && R==r)
         return c[rt];

     ;

     if (R <= mid)
         return Query_Rank(lson[rt], L, R, l, mid);
     else if (L > mid)
         , r);
     else
         , R, mid+, r);
 }

 void solve() {
     __int64 ans1 = , ans2 = , ans3 = ;
     int tmp, id;

     sort(a, a+m);
     m = unique(a, a+m) - a;
     T[cur] = Build(, m-);

     rep(i, , q) {
         ) {
             id = lower_bound(a, a+m, Q[i].k) - a;
             T[cur+] = Insert(T[cur], id, );
             ++cur;
             continue;
         }

         ) {
             id = Query_kth(T[Q[i].l-], T[Q[i].r], Q[i].k);
             tmp = a[id];
             ans1 += tmp;
         } ) {
             id = lower_bound(a, a+m, Q[i].k) - a;
             tmp = Query_Rank(T[cur], , id, , m-);
             ans2 += tmp;
         } else {
             id = Query_kth(T[], T[cur], Q[i].k);
             tmp = a[id];
             ans3 += tmp;
         }

         #ifndef ONLINE_JUDGE
             // printf("tmp = %d\n", tmp);
         #endif
     }

     printf("%I64d\n%I64d\n%I64d\n", ans1, ans2, ans3);
 }

 int main() {
     ios::sync_with_stdio(false);
     #ifndef ONLINE_JUDGE
         freopen("data.in", "r", stdin);
         freopen("data.out", "w", stdout);
     #endif

     ;
     ];

     while (scanf("%d", &q)!=EOF) {
         init();
         rep(i, , q) {
             scanf("%s", op);
             ] == 'I') {
                 scanf("%d", &Q[i].k);
                 a[m++] = Q[i].k;
                 Q[i].q = ;
             } ] == ') {
                 scanf("%d %d %d", &Q[i].l, &Q[i].r, &Q[i].k);
                 Q[i].q = ;
             } ] == ') {
                 scanf("%d", &Q[i].k);
                 Q[i].q = ;
             } else {
                 scanf("%d", &Q[i].k);
                 Q[i].q = ;
             }
         }
         printf("Case %d:\n", ++t);
         solve();
     }

     #ifndef ONLINE_JUDGE
         printf("time = %d.\n", (int)clock());
     #endif

     ;
 }

数据发生器。

 from copy import deepcopy
 from random import randint, shuffle
 import shutil
 import string

 def GenDataIn():
     with open("data.in", "w") as fout:
         t = 10
         bound = 10**5
         for tt in xrange(t):
             n = randint(200, 300)
             fout.write("%d\n" % (n))
             L = range(1, n+1)
             shuffle(L)
             cur = 0
             has = []
             for i in xrange(n):
                 op = randint(0, 3)
                 if cur==0:
                     op=0
                 if op==0:
                     idx = randint(0, len(L)-1)
                     x = L[idx]
                     L.remove(x)
                     has.append(x)
                     cur += 1
                     fout.write("Insert %d\n" % x)
                 elif op==1:
                     l = randint(1, cur)
                     r = randint(l, cur)
                     k = randint(1, r-l+1)
                     fout.write("Query_1 %d %d %d\n" % (l, r, k))
                 elif op==2:
                     idx = randint(0, len(has)-1)
                     x = has[idx]
                     fout.write("Query_2 %d\n" % x)
                 else:
                     x = randint(1, cur)
                     fout.write("Query_3 %d\n" % x)

 def MovDataIn():
     desFileName = "F:\eclipse_prj\workspace\hdoj\data.in"
     shutil.copyfile("data.in", desFileName)

 if __name__ == "__main__":
     GenDataIn()
     MovDataIn()

【HDOJ】3727 Jewel的更多相关文章

  1. 【HDOJ】4729 An Easy Problem for Elfness

    其实是求树上的路径间的数据第K大的题目.果断主席树 + LCA.初始流量是这条路径上的最小值.若a<=b,显然直接为s->t建立pipe可以使流量最优:否则,对[0, 10**4]二分得到 ...

  2. 【HDOJ】【3506】Monkey Party

    DP/四边形不等式 裸题环形石子合并…… 拆环为链即可 //HDOJ 3506 #include<cmath> #include<vector> #include<cst ...

  3. 【HDOJ】【3516】Tree Construction

    DP/四边形不等式 这题跟石子合并有点像…… dp[i][j]为将第 i 个点开始的 j 个点合并的最小代价. 易知有 dp[i][j]=min{dp[i][j] , dp[i][k-i+1]+dp[ ...

  4. 【HDOJ】【3480】Division

    DP/四边形不等式 要求将一个可重集S分成M个子集,求子集的极差的平方和最小是多少…… 首先我们先将这N个数排序,容易想到每个自己都对应着这个有序数组中的一段……而不会是互相穿插着= =因为交换一下明 ...

  5. 【HDOJ】【2829】Lawrence

    DP/四边形不等式 做过POJ 1739 邮局那道题后就很容易写出动规方程: dp[i][j]=min{dp[i-1][k]+w[k+1][j]}(表示前 j 个点分成 i 块的最小代价) $w(l, ...

  6. 【HDOJ】【3415】Max Sum of Max-K-sub-sequence

    DP/单调队列优化 呃……环形链求最大k子段和. 首先拆环为链求前缀和…… 然后单调队列吧<_<,裸题没啥好说的…… WA:为毛手写队列就会挂,必须用STL的deque?(写挂自己弱……s ...

  7. 【HDOJ】【3530】Subsequence

    DP/单调队列优化 题解:http://www.cnblogs.com/yymore/archive/2011/06/22/2087553.html 引用: 首先我们要明确几件事情 1.假设我们现在知 ...

  8. 【HDOJ】【3068】最长回文

    Manacher算法 Manacher模板题…… //HDOJ 3068 #include<cstdio> #include<cstring> #include<cstd ...

  9. 【HDOJ】【1512】Monkey King

    数据结构/可并堆 啊……换换脑子就看了看数据结构……看了一下左偏树和斜堆,鉴于左偏树不像斜堆可能退化就写了个左偏树. 左偏树介绍:http://www.cnblogs.com/crazyac/arti ...

随机推荐

  1. JDK环境变量配置说明

    摘要:被人问到,为什么要配置Path/ClassPath JAVA_HOME,突然感觉有些讲不清楚."新人有资格问一个怪问题,但是老鸟不能给一个烂回答".所以今天为了让自己进一步学 ...

  2. c#延迟加载

    public class BlogUser { public int Id { get; private set; } public Lazy<List<Article>> A ...

  3. C#实现DNS解析服务和智能DNS服务

    C#实现DNS解析服务有一个开源项目ARSoft.Tools.Net, ARSoft.Tools.Net是一个非常强大的开源DNS控件库,包含.Net SPF validation, SenderID ...

  4. ACM:SCU 4437 Carries - 水题

    SCU 4437  Carries Time Limit:0MS     Memory Limit:0KB     64bit IO Format:%lld & %llu  Practice  ...

  5. SQLite的原子提交原理

    本文描述了sqlite为保证数据库文件不被损坏而采取的种种手段.. 以下是原译者的摘要:http://www.kuqin.com/shuoit/20150618/346693.html 摘要: 本文源 ...

  6. table清除样式大全

    table{width:100%;text-align:center;border-collapse:collapse;border-spacing:1;border-spacing:0; }tabl ...

  7. Smarty缓存的5个知识点

    (1)页面缓存:整个页面全局的缓存 需要4个步骤: ①开启缓存  $smarty->caching = true; ②设置缓存的生命周期  $smarty->cache_lifetime ...

  8. java文档

    http://www.boyunjian.com/javadoc/com.dyuproject.protostuff/protostuff-me/1.0.5/_/com/dyuproject/prot ...

  9. nodejs入门(一)

    1.nodejs简介  Nodejs不是一个js应用.而是一个js运行平台. Node.js 使用事件驱动, 非阻塞I/O 模型而得以轻量和高效 2. 1).Nodejs内置了一个HTTP模块 var ...

  10. 【Spring】Spring的定时任务

    > 参考的优秀文章 Task Execution and Scheduling > 版本说明 <dependencies> <dependency> <gro ...