POJ 2886 Who Gets the Most Candies? （三种姿势点草动态求第k大）

最新推荐文章于 2019-05-03 13:52:16 发布
最新推荐文章于 2019-05-03 13:52:16 发布
阅读量683 收藏 1
点赞数
分类专栏： 数据结构 - 线段树 数据结构 文章标签： 动态求第k大
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/lwt36/article/details/46527141
版权

Who Gets the Most Candies?

题意: 

N个孩子顺时针坐成一个圆圈且从1到N编号，每个孩子手中有一张标有非零整数的卡片。第K个孩子先出圈，如果他手中卡片上的数字A大于零，下一个出圈的是他左手边第A个孩子。否则，下一个出圈的是他右手边第(-A)个孩子。第p个出圈的孩子会得到F(p)个糖果，F(p)为p的因子数。求得到糖果数最多的是哪个孩子及得到多少糖果。

思路:

模拟约瑟夫环的感觉 注意处理取模就好了 最好在纸上画一下

坑: 

从1开始的时候  - - 千万千万别忘记了取模成0的时候 其实是模数啊 - - (又让我想起了某一次无缘D题 差1s就是这个原因 刚交瞬间就发现了）

负数加模数国际惯例都懂的


树状数组(logn)^2查询做法:

附本题一个很棒的case

//
//  Created by TaoSama on 2015-06-06
//  Copyright (c) 2015 TaoSama. All rights reserved.
//
//#pragma comment(linker, "/STACK:1024000000,1024000000")
#include <algorithm>
#include <cctype>
#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iomanip>
#include <iostream>
#include <map>
#include <queue>
#include <string>
#include <set>
#include <vector>

#define debug(x) printf(#x": %d\n", x);

using namespace std;
const int INF = 0x3f3f3f3f;
const int MOD = 1e9 + 7;
const int N = 5e5 + 10;

int n, k, d[N];

void init() {
    for(int i = 1; i <= 5e5; ++i) d[i] = 1;
    for(int i = 2; i <= 5e5; ++i)
        for(int j = i; j <= 5e5; j += i)
            d[j]++;
}

int b[N];

void add(int i, int x) {
    for(; i <= 5e5; i += i & -i) b[i] += x;
}

int sum(int i) {
    int ret = 0;
    for(; i; i -= i & -i) ret += b[i];
    return ret;
}

int binSearch(int k) {
    int l = 1, r = 5e5;
    while(l <= r) {
        int m = l + r >> 1;
        if(sum(m) < k) l = m + 1;
        else r = m - 1;
    }
    return l;
}

int a[N];
char name[N][15];

int main() {
#ifdef LOCAL
    freopen("in.txt", "r", stdin);
//  freopen("out.txt","w",stdout);
#endif
    //ios_base::sync_with_stdio(0);

    init();
    while(scanf("%d%d", &n, &k) == 2) {
        for(int i = 1; i <= n; ++i) {
            scanf("%s%d", name[i], a + i);
            add(i, 1);
        }
        int ans, p = k, Max = -INF;  //p真实位置 k名次
        for(int i = 1; i <= n; ++i) {
            add(p, -1);
            if(d[i] > Max) Max = d[i], ans = p;
            int mod = n - i;
            if(!mod) break;
//            cout << "mod: " << mod << " p: " << p << ' ' << sum(p) << endl;
            if(a[p] > 0) k = (sum(p) + a[p]) % mod;
            else k = ((sum(p) + a[p]) % mod + mod) % mod + 1;
//            cout << ((sum(p) + a[p]) % mod + mod) % mod << endl;
            if(k == 0) k = mod;
            p = binSearch(k);
//            cout << "k: " << k << " p: " << p << endl;
        }
        printf("%s %d\n", name[ans], Max);
    }
    return 0;
}

/*
7 3
a 3
b 2
c -5
d 4
e 8
f 2
g -6

ans: g 4
*/

线段树logn查询做法:

用到了反素数打表 - - 我会说光看懂这个东西就用了好几天

详情见 http://blog.csdn.net/lwt36/article/details/46526987

//
//  Created by TaoSama on 2015-06-11
//  Copyright (c) 2015 TaoSama. All rights reserved.
//
//#pragma comment(linker, "/STACK:1024000000,1024000000")
#include <algorithm>
#include <cctype>
#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iomanip>
#include <iostream>
#include <map>
#include <queue>
#include <string>
#include <set>
#include <vector>

using namespace std;
const int INF = 0x3f3f3f3f;
const int MOD = 1e9 + 7;
const int N = 5e5 + 10;

typedef long long LL;

int p[] = {2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53};
map<LL, LL> mp, rmp; //divisors .. anti-prime
LL a[100], b[100], ct;

void dfs(int k, LL cur, LL cnt, int limit) {
    if(cur > N) return;
    if(!mp.count(cnt)) mp[cnt] = cur;
    else mp[cnt] = min(cur, mp[cnt]);
    for(int i = 1; i <= limit; ++i) {
        if(1.0 * cur * p[k] > N) break;
        dfs(k + 1, cur *= p[k], cnt * (i + 1), i);
    }
}

void init() {
    dfs(0, 1, 1, 60);
    map<LL, LL>::iterator i = mp.begin();
    for(; i != mp.end(); ++i) rmp.insert(make_pair(i->second, i->first));
    i = rmp.begin();
    for(; i != rmp.end(); ++i) {
        if(ct > 0 && i->second < b[ct - 1]) continue;
        a[ct] = i->first;
        b[ct++] = i->second;
    }
}

int n, k, v[N];
char name[N][10];

int sum[N * 3];

#define root 1, n, 1
#define lson l, m, rt << 1
#define rson m + 1, r, rt << 1 | 1

void build(int l, int r, int rt){
    sum[rt] = r - l + 1;
    if(l == r) return;
    int m = l + r >> 1;
    build(lson);
    build(rson);
}

int query(int k, int l, int r, int rt){
    sum[rt]--;
    if(l == r) return l;
    int m = l + r >> 1;
    if(k <= sum[rt << 1]) return query(k, lson);
    else return query(k - sum[rt << 1], rson);
}

int main() {
#ifdef LOCAL
    freopen("in.txt", "r", stdin);
//  freopen("out.txt","w",stdout);
#endif
    ios_base::sync_with_stdio(0);

    init();
    while(scanf("%d%d", &n, &k) == 2) {
        build(root);
        for(int i = 1; i <= n; ++i) scanf("%s%d", name[i], v + i);

        int ans, idx = 0;
        while(a[idx] <= n && idx < ct) idx++;
        int cnt = a[idx - 1], Max = b[idx - 1]; //第cnt个就是最大的

        for(int i = 1; i <= cnt; ++i) {
            int mod = n - i;
            ans = query(k, root);
            if(mod == 0) break;
            if(v[ans] > 0) k = (k - 1 + v[ans]) % mod;
            else k = ((k - 1 + v[ans]) % mod + mod) % mod + 1;
            if(k == 0) k = mod;
        }
        printf("%s %d\n", name[ans], Max);
    }
    return 0;
}

树状数组logn查询做法:

部分内容来自 http://www.cnblogs.com/wuyiqi/archive/2011/12/25/2301071.html


现在假设要求sum[x]的值，一般我们都是从后往前求和，如x=15
15-lowbit(15)=14;
14-lowbit(14)=12;
12-lowbit(12)=8;
8-lowbit(b)=0;
答案就是sum[15]+sum[14]+sum[12]+sum[8];

现在我们可以这样来求，从不超过15的只有一个1的最大二进制数，也可以理解为指数从log(15)取整=3开始

2^3=8, 加上2^2=4, 加上2^1=2, 加上2^0=1 数字依次为8,12,14,15，也就是把普通的求和过程反向。

int find_kth(int k) {
    int m = 31 - __builtin_clz(n);  //求32位log2(n)
    int cnt = 0, ret = 0;
    for(int i = m; i >= 0; --i) {
        ret += 1 << i;
        if(ret >= n || cnt + bit[ret] >= k) ret -= 1 << i;
        //可能会有很多个数都满足cnt+bit[ret]>=k 找的最大的ret满足cnt+bit[ret]<k
        else cnt += bit[ret];
        //cnt累加比当前ans小的总数
    }
    //ret是cnt（即小于等于ret的数的个数）小于k的情况的最大值 ret+1即第k大的数
    return ret + 1;
}



参考code: 

//
//  Created by TaoSama on 2015-06-11
//  Copyright (c) 2015 TaoSama. All rights reserved.
//
//#pragma comment(linker, "/STACK:1024000000,1024000000")
#include <algorithm>
#include <cctype>
#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iomanip>
#include <iostream>
#include <map>
#include <queue>
#include <string>
#include <set>
#include <vector>

using namespace std;
const int INF = 0x3f3f3f3f;
const int MOD = 1e9 + 7;
const int N = 5e5 + 10;

typedef long long LL;

int p[] = {2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53};
map<LL, LL> mp, rmp; //divisors .. anti-prime
LL a[100], b[100], ct; //a: anti-prime b: its divisors

void dfs(int k, LL cur, LL cnt, int limit) {
    if(cur > N) return;
    if(!mp.count(cnt)) mp[cnt] = cur;
    else mp[cnt] = min(cur, mp[cnt]);
    for(int i = 1; i <= limit; ++i) {
        if(1.0 * cur * p[k] > N) break;
        dfs(k + 1, cur *= p[k], cnt * (i + 1), i);
    }
}

//when try to iterate, do not forget to control the bound (< ct)
void init() {
    dfs(0, 1, 1, 60);
    map<LL, LL>::iterator i = mp.begin();
    for(; i != mp.end(); ++i) rmp.insert(make_pair(i->second, i->first));
    i = rmp.begin();
    for(; i != rmp.end(); ++i) {
        if(ct > 0 && i->second < b[ct - 1]) continue;
        a[ct] = i->first;
        b[ct++] = i->second;
    }
}

int n, k, v[N];
char name[N][10];

int bit[N];

void add(int i, int x) {
    for(; i <= n; i += i & -i) bit[i] += x;
}

int find_kth(int k) {
    int m = 31 - __builtin_clz(n);
    int cnt = 0, ret = 0;
    for(int i = m; i >= 0; --i) {
        ret += 1 << i;
        if(ret >= n || cnt + bit[ret] >= k) ret -= 1 << i;
        else cnt += bit[ret];
    }
    return ret + 1;
}

int main() {
#ifdef LOCAL
    freopen("in.txt", "r", stdin);
//  freopen("out.txt","w",stdout);
#endif
    ios_base::sync_with_stdio(0);

    init();
    while(scanf("%d%d", &n, &k) == 2) {
        memset(bit, 0, sizeof bit);
        for(int i = 1; i <= n; ++i) {
            scanf("%s%d", name[i], v + i);
            add(i, 1);
        }

        int ans, idx = 0;
        while(a[idx] <= n && idx < ct) idx++;
        int cnt = a[idx - 1], Max = b[idx - 1]; //第cnt个就是最大的
        for(int i = 1; i <= cnt; ++i) {
            int mod = n - i;
            ans = find_kth(k);
            if(mod == 0) break;
            add(ans, -1);
            if(v[ans] > 0) k = (k - 1 + v[ans]) % mod;
            else k = ((k - 1 + v[ans]) % mod + mod) % mod + 1;
            if(k == 0) k = mod;
        }
        printf("%s %d\n", name[ans], Max);
    }
    return 0;
}


TaoSama
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
矩形面积并、矩形面积交、矩形周长并（线段树、扫描线总结）
TaoSama's Microcosm
10-05 1万+
HDU 1542 [POJ 1151] Atlantis （矩形面积并） 题意: 求N<=100个矩形的面积并求N<=100个矩形的面积并分析: 离散化: 这些技巧都是老生常谈的了, 不然浮点数怎么建树, 离散化xx坐标就可以了 扫描线: 首先把矩形按yy轴分成两条边, 上边和下边, 对xx轴建树, 扫描线可以看成一根平行于xx轴的直线. 从y=0y=0开始往上扫, 下边表示要计
poj2886 Who Gets the Most Candies?
emmmmmm
07-17 2227
题意：输入n，k表示n个人，第k个第一个离开圆环，之后n行，每行一个名字和一个数字m，表示从这个人开始顺时针（数字大于0）或逆时针（数字小于0）第m个人离开圆环。             每个人离开时会有一个得分，得分=F(离开次序i); F(i)=能整除i的数的个数；求这个人的名字和得分。 思路：首先打表枚举[1，500000]的F(i);  然后模拟求这个人的名字，对于[1，n]中
poj 2886 Who Gets the Most Candies? （树状数组+二分+反素数）
舒哥的blog
07-26 887
Who Gets the Most Candies? Time Limit: 5000MS   Memory Limit: 131072K Total Submissions: 11597   Accepted: 3616 Case Time Limit: 2000MS Description N children are si
POJ2886 Who Gets the Most Candies?【线段树 点修改】
逐梦者
05-03 208
Who Gets the Most Candies? http://poj.org/problem?id=2886 Time Limit:5000MS Memory Limit:131072K Total Submissions:16835 Accepted:5298 Case Time Limit:2000MS Description ...
poj2886--Who Gets the Most Candies?(线段树+反素数)
随风的专栏
07-11 991
题目链接点击打开链接 题目大意：给出n个人的姓名和手里的一个号码，n个人排成一圈，号码有正有负，代表着正向还是反向移动k个位置，比赛从第k个人开始，把被选到的人踢出，问按踢出的顺序中因子数最多的是谁？ 建立线段树，把n个人被踢的顺序找到，然后求出n个人中因子数最多的（最小的数）是谁，这里要用到反素数，详看链接点击打开链接 #include #include #include using
POJ2092:计数排序,求第K大的元素
04-15
【标题】"POJ2092:计数排序,求第K大的元素"是一个编程题目，主要涉及计数排序算法以及如何在数组中找出第K大的元素。计数排序是一种非基于比较的排序算法，它适用于整数排序，尤其在数据范围不大的情况下效率极高。...
POJ2388-Who's in the Middle
07-31
《POJ2388-Who's in the Middle》是一个典型的计算机编程竞赛题目，主要考察参赛者的算法设计和实现能力。这个题目源自北京大学的在线评测系统POJ（PKU Online Judge），是编程爱好者和学生提升编程技能的重要平台之...
西工大C语言POJ习题答案.docx
最新发布
11-15
西工大C语言POJ习题答案 本资源提供了一系列C语言POJ习题的答案，涵盖了基本编程、数据类型、运算符、控制结构、函数、数组、指针等多个方面。下面是对每道题目的解释和分析： 1. 第一道题目：两个整数的和 这道...
POJ3273.rar_M?n
09-23
《POJ3273 Monthly Expense 题解——数据结构与动态规划的应用》 POJ3273，这是一道经典的算法竞赛题目，它涉及到的问题是：如何有效地合并连续的数，使得合并后的集合数量不超过M个，并且这些集合的和尽可能小。这...
西工大poj第一季
12-04
poj习题第一题，很有意思，是学习poj的开始
HDU 5493 Queue （线段树） 2015合肥网络赛
TaoSama's Microcosm
09-28 1531
题意: 给定N<105个人的序列,N个(height,k)二元组描述这个序列height:=这个人的身高,k:=这个人的左边或者右边有k个人比他高构造一个字典序最小的序列满足这些条件给定N<10^5个人的序列, N个(height, k)二元组描述这个序列 height:= 这个人的身高, k:= 这个人的左边或者右边有k个人比他高 构造一个字典序最小的序列满足这些条件 分析:
线段树小结
TaoSama's Microcosm
12-10 1215
线段树小结 A Summary for Segment Tree线段树小结\ A\ Summary\ for\ Segment\ Tree0. Anouncement0.\ Anouncement 本文部分图片以及部分内容来自互联网，内容过多就不一一注明出处了，冒犯之处还请海涵。本文部分图片以及部分内容来自互联网，内容过多就不一一注明出处了，冒犯之处还请海涵。 Some of the
Codeforces #321 E. Kefa and Watch （线段树、哈希）
TaoSama's Microcosm
09-23 946
题意: 给定数字序列 2种操作           1. 将[l,r]区间数字变成c           2. 查询[l,r]区间数字串是不是周期为d的串            x is called a period of string s (1 ≤ x ≤ |s|), if si  =  si + x for all i from 1 to |s|  -  x. 分析:
UVALive 3938 Ray, Pass me the dishes （线段树区间合并）
TaoSama's Microcosm
09-16 940
题意: 求q次询问的静态区间连续最大和起始位置和终止位置 输出字典序最小的解 分析: 很显然的区间合并的题目啦   常理就是维护pre, sub, suf 最大的 前缀和, 连续和, 后缀和 然后这个题还要记录解的位置所以还要区间总和sum   然后解的位置preR, subL, subR, sufL 然后更新就好了 自己注意下字典序就好了 建树和查询都类似 总结: 死守着HH线段树风格
ACdream 简单数据结构 专题
TaoSama's Microcosm
09-22 880
A - 风之国 题意: 给出X轴上N个点的坐标   X轴上的点按大小两两互相连通 现在要必须干掉K条路 求最小花费 分析: dp+线段树优化 我表示不会 copy一遍题解          首先题目给出的城市坐标不是按X轴升序的,  那么我们就将每个城市根据坐标从左到右映射到X轴上的1~n          设dp[i]:= 到第 i 个点(包括第i个点)时,  处理掉前面所有必须不连
HDU 4267 A Simple Problem with Integers （线段树）
TaoSama's Microcosm
09-30 700
题意: N<=5×104,每次将区间[a,b]内满足(i−a)%k=0(i∈[a,b])的数Ai增加c,最后单点查询值N <= 5×10^4, 每次将区间[a,b]内满足(i-a)\%k=0(i \in[a,b])的数A_i增加c, 最后单点查询值 分析: 很裸的区间增减,单点查询的模型,但是好久没写BIT了我没有看出来,就想着线段树去了,然后只有单点查询的,我维护了区间增减值的和。。这题
HDU 5497 Inversion （BIT）
TaoSama's Microcosm
10-04 696
题意: N<=105的序列,求删去一个M<N的连续子序列后,剩余序列的最小逆序对数N<=10^5的序列, 求删去一个M<N的连续子序列后, 剩余序列的最小逆序对数 分析: 求逆序对当然要BIT辣,然后维护删去固定大小的序列可以用two pointers求逆序对当然要BIT辣, 然后维护删去固定大小的序列可以用two \ pointers 开两个BIT,b[]维护连续M子序列前面的,c[
Codeforces #310（div2）
TaoSama's Microcosm
06-28 655
A. Case of the Zeros and Ones 题意: 我爱消除 - - 相邻的01就会被消除 问最后能剩下多少~ 思路: 仔细观察我们发现 只要存在不相同的数 那么它们一定会被消除 最后一定为全0或者全1序列 那么答案就是abs(n_'0' - n_'1')了 参考code: // // Created by TaoSama on 2015-06-27 // Cop

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值