首页 > 代码库 > 【bzoj4721】[Noip2016]蚯蚓
【bzoj4721】[Noip2016]蚯蚓
题目描述
本题中,我们将用符号[c]表示对c向下取整,例如:[3.0」= [3.1」=[3.9」=3。蛐蛐国最近蚯蚓成灾了!隔壁跳蚤国的跳蚤也拿蚯蚓们没办法,蛐蛐国王只好去请神刀手来帮他们消灭蚯蚓。蛐蛐国里现在共有n只蚯蚓(n为正整数)。每只蚯蚓拥有长度,我们设第i只蚯蚓的长度为a_i(i=1,2,...,n),并保证所有的长度都是非负整数(即:可能存在长度为0的蚯蚓)。每一秒,神刀手会在所有的蚯蚓中,准确地找到最长的那一只(如有多个则任选一个)将其切成两半。神刀手切开蚯蚓的位置由常数p(是满足0<p<1的有理数)决定,设这只蚯蚓长度为x,神刀手会将其切成两只长度分别为[px]和x-[px]的蚯蚓。特殊地,如果这两个数的其中一个等于0,则这个长度为0的蚯蚓也会被保留。此外,除了刚刚产生的两只新蚯蚓,其余蚯蚓的长度都会增加q(是一个非负整常数)。蛐蛐国王知道这样不是长久之计,因为蚯蚓不仅会越来越多,还会越来越长。蛐蛐国王决定求助于一位有着洪荒之力的神秘人物,但是救兵还需要m秒才能到来......(m为非负整数)蛐蛐国王希望知道这m秒内的战况。具体来说,他希望知道:?m秒内,每一秒被切断的蚯蚓被切断前的长度(有m个数)?m秒后,所有蚯蚓的长度(有n+m个数)。蛐蛐国王当然知道怎么做啦!但是他想考考你......
输入
第一行包含六个整数n,m,q,u,v,t,其中:n,m,q的意义见问题描述;
u,v,t均为正整数;你需要自己计算p=u/v(保证0<u<v)t是输出参数,其含义将会在输出格式中解释。
第二行包含n个非负整数,为ai,a2,...,an,即初始时n只蚯蚓的长度。
同一行中相邻的两个数之间,恰好用一个空格隔开。
保证1<=n<=10^5,0<m<7*10^6,0<u<v<10^9,0<=q<=200,1<t<71,0<ai<10^8。
输出
第一行输出[m/t]个整数,按时间顺序,依次输出第t秒,第2t秒,第3t秒……被切断蚯蚓(在被切断前)的长度。
第二行输出[(n+m)/t]个整数,输出m秒后蚯蚓的长度;需要按从大到小的顺序
依次输出排名第t,第2t,第3t……的长度。
同一行中相邻的两个数之间,恰好用一个空格隔开。即使某一行没有任何数需要 输出,你也应输出一个空行。
请阅读样例来更好地理解这个格式。
样例输入
3 7 1 1 3 1
3 3 2
样例输出
3 4 4 4 5 5 6
6 6 6 5 5 4 4 3 2 2
题解
队列题
由于先切的两段一定大于等于后切的对应的两段,
于是可以用两个队列分别存储[px]段和x-[px]段,还有一个存储原始长度。
每次把所有没切的长度+q,相当于把切了的长度-q,然后查找时再一起加上。
注意除法需要long long。
还有,在bzoj上,可能会PE无数次,只能按照样例格式输出,不能多任何的空格和换行。
1 #include <cstdio> 2 #include <algorithm> 3 #define inf 0x7fffffff 4 using namespace std; 5 int q[3][8000001] , l[3] , r[3] = {-1 , -1 , -1} , a[100001]; 6 bool cmp(int a , int b) 7 { 8 return a > b; 9 } 10 int getmax() 11 { 12 int k = -inf , ans = 0; 13 if(l[0] <= r[0] && q[0][l[0]] > k) k = q[0][l[0]] , ans = 0; 14 if(l[1] <= r[1] && q[1][l[1]] > k) k = q[1][l[1]] , ans = 1; 15 if(l[2] <= r[2] && q[2][l[2]] > k) k = q[2][l[2]] , ans = 2; 16 return ans; 17 } 18 int main() 19 { 20 int n , m , g , u , v , t , i , p , c; 21 scanf("%d%d%d%d%d%d" , &n , &m , &g , &u , &v , &t); 22 for(i = 0 ; i < n ; i ++ ) 23 scanf("%d" , &a[i]); 24 sort(a , a + n , cmp); 25 for(i = 0 ; i < n ; i ++ ) 26 q[0][++r[0]] = a[i]; 27 for(i = 1 ; i <= m ; i ++ ) 28 { 29 p = getmax(); 30 c = q[p][l[p]++] + g * (i - 1); 31 if(i % t == 0) 32 { 33 if(i != t) 34 printf(" "); 35 printf("%d" , c); 36 } 37 q[1][++r[1]] = (int)((long long)c * u / v) - g * i; 38 q[2][++r[2]] = c - (int)((long long)c * u / v) - g * i; 39 } 40 printf("\n"); 41 for(i = 1 ; i <= m + n ; i ++ ) 42 { 43 p = getmax(); 44 c = q[p][l[p]++] + g * m; 45 if(i % t == 0) 46 { 47 if(i != t) 48 printf(" "); 49 printf("%d" , c); 50 } 51 } 52 printf("\n"); 53 return 0; 54 }
【bzoj4721】[Noip2016]蚯蚓
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。