当前位置: 代码迷 >> 综合 >> 【NOIP2016】洛谷2827 蚯蚓
  详细解决方案

【NOIP2016】洛谷2827 蚯蚓

热度:13   发布时间:2024-01-13 11:16:36.0

题目描述

本题中,我们将用符号[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……的长度。

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

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

后切的一定没有先切的长,所以可以开三个队列分别存原来的和切出来的两个,每次从三个队头取出最大值。
最开始要排序,之后每次取出和放回都是O(1)的,所以复杂度O(nlogn+m)。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int rd()
{int x=0;char c=getchar();while (c<'0'||c>'9') c=getchar();while (c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}return x;
}
int buf[15];
void wt(int x)
{int l=0;if (!x) putchar('0');while (x){buf[++l]=x%10;x/=10;}while (l--)putchar(buf[l+1]+'0');
}
bool cmp(int x,int y)
{return x>y;
}
int m,n,q,u,v,t,a1[8000010],a2[8000010],a3[8000010];
int main()
{int i,j,k,x,y,z,K,h1,t1,h2,t2,h3,t3;bool flag;scanf("%d%d%d%d%d%d",&n,&m,&q,&u,&v,&t);for (i=1;i<=n;i++)a1[i]=rd();sort(a1+1,a1+n+1,cmp);h1=1;t1=n;h2=h3=1;t2=t3=0;flag=0;for (K=1;K<=m;K++){if (h1<=t1&&(h2>t2||a1[h1]>a2[h2])&&(h3>t3||a1[h1]>a3[h3]))x=a1[h1++];else{if (h2<=t2&&(h3>t3||a2[h2]>a3[h3]))x=a2[h2++];elsex=a3[h3++];}x+=(K-1)*q;if (K%t==0){if (!flag) flag=1;else putchar(' ');wt(x);}y=(double)x*u/v;z=x-y;a2[++t2]=y-K*q;a3[++t3]=z-K*q;}putchar('\n');flag=0;for (K=1;K<=m+n;K++){if (h1<=t1&&(h2>t2||a1[h1]>a2[h2])&&(h3>t3||a1[h1]>a3[h3]))x=a1[h1++];else{if (h2<=t2&&(h3>t3||a2[h2]>a3[h3]))x=a2[h2++];elsex=a3[h3++];}if (K%t==0){if (!flag) flag=1;else putchar(' ');wt(x+m*q);}}putchar('\n');
}