NOIP2016D2T2 蚯蚓

Posted

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了NOIP2016D2T2 蚯蚓相关的知识,希望对你有一定的参考价值。

题目描述

本题中,我们将用符号?c?表示对c向下取整,例如:?3.0?=?3.1?=?3.9?=3。

蛐蛐国最近蚯蚓成灾了!隔壁跳蚤国的跳蚤也拿蚯蚓们没办法,蛐蛐国王只好去请神刀手来帮他们消灭蚯蚓。

蛐蛐国里现在共有n只蚯蚓(n为正整数)。每只蚯蚓拥有长度,我们设第i只蚯蚓的长度为a_i(i=1,2,...,n)ai?(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个非负整数,为a_i,a_2,...,a_nai?,a2?,...,an?,即初始时n只蚯蚓的长度。

同一行中相邻的两个数之间,恰好用一个空格隔开。

保证1n10^5,0<m7?10^6,0u<v10^9,0q200,1t71,0<ai10^8。

 

输出格式:

第一行输出?m/t?个整数,按时间顺序,依次输出第t秒,第2t秒,第3t秒……被切断蚯蚓(在被切断前)的长度。

第二行输出?(n+m)/t?个整数,输出m秒后蚯蚓的长度;需要按从大到小的顺序,依次输出排名第t,第2t,第3t……的长度。

同一行中相邻的两个数之间,恰好用一个空格隔开。即使某一行没有任何数需要 输出,你也应输出一个空行。

请阅读样例来更好地理解这个格式。

算法分析

首先想到是优先队列,显然超时了。所以必须使用O(n)的算法,也就意味着要实现一个不需要更新放进去就是有序的队列算法。

所以想到构建3个队列,一个放原本的有序数列,其余两个分别放被切开的较大部分和较小部分。

为什么是有序的呢,第一个显而易见,后两个只要证明a[i]*p+q>(a[i+1]+q)*p,已知a[i]>a[i+1],0<p<1,q>0,所以q>p*q,易得原式。

所以就可以直接实现了。但是我为什么打了这么丑的一个代码呢(。。)

代码

 1 #include<bits/stdc++.h>
 2 using namespace std;
 3 int n,m,d,u,v,t,tail1=0,tail2=0;
 4 long long q[3][7000010],tim[3][7000010],a[100010];
 5 
 6 bool cmp(const long long& a,const long long& b){
 7     return a>b;
 8 }
 9 
10 int main(){
11     scanf("%d%d%d%d%d%d",&n,&m,&d,&u,&v,&t);
12     for (int i=1; i<=n; i++)
13         scanf("%lld",&a[i]);
14     sort(a+1,a+n+1,cmp);
15     for (int i=1; i<=n; i++){
16         q[0][i]=a[i];
17         tim[0][i]=1;
18     }
19     long long len;
20     int i=1,j=1,k=1;
21     for (int T=1; T<=m; T++){
22         long long l0=q[0][i]+d*(T-tim[0][i]),l1=q[1][j]+d*(T-tim[1][j]),l2=q[2][k]+d*(T-tim[2][k]);
23         if (i>n) l0=-1; if (j>tail1) l1=-1; if (k>tail2) l2=-1;
24         if (l0>=l1&&l0>=l2){
25             len=l0; i++;
26         }
27             else if (l1>=l0&&l1>=l2){
28                 len=l1; j++;
29             }
30                 else{
31                     len=l2; k++;
32                 }
33         if (T%t==0) printf("%lld ",len);
34         q[1][++tail1]=(int)(len*(double)u/(double)v);
35         tim[1][tail1]=T+1;
36         q[2][++tail2]=len-q[1][tail1];
37         tim[2][tail2]=T+1;
38     }
39     printf("\n");
40     for (int T=1; T<=n+m; T++){
41         long long l0=q[0][i]+d*(m-tim[0][i]+1),l1=q[1][j]+d*(m-tim[1][j]+1),l2=q[2][k]+d*(m-tim[2][k]+1);
42         if (i>n) l0=-1; if (j>tail1) l1=-1; if (k>tail2) l2=-1;
43         if (l0>=l1&&l0>=l2) len=l0,i++;
44             else if (l1>=l0&&l1>=l2) len=l1,j++;
45                 else len=l2,k++;
46         if (T%t==0) printf("%lld ",len);
47     }
48     printf("\n");
49     return 0;
50 } 

 

以上是关于NOIP2016D2T2 蚯蚓的主要内容,如果未能解决你的问题,请参考以下文章

[Noip2016]蚯蚓 D2 T2 队列

NOIP考前注意事项

NOIP2016D2T1 组合数问题

[noip2013 d2t2]花匠

NOIP2015 D2T2 子串

codevs 4560 NOIP2015 D2T2 子串