蚯蚓_NOIP2016D2T2_单调队列

48 阅读 0 评论 32 点赞

我是靠谱客的博主活力豌豆，最近开发中收集的这篇文章主要介绍蚯蚓_NOIP2016D2T2_单调队列，觉得挺不错的，现在分享给大家，希望可以做个参考。

概述

Description

本题中，我们将用符号[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个数)。蛐蛐国王当然知道怎么做啦！但是他想考考你……

Input

第一行包含六个整数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。

Output

第一行输出[m/t]个整数，按时间顺序，依次输出第t秒，第2t秒，第3t秒……被切断蚯蚓（在被切断前）的长度。
第二行输出[(n+m)/t]个整数，输出m秒后蚯蚓的长度；需要按从大到小的顺序
依次输出排名第t，第2t，第3t……的长度。
同一行中相邻的两个数之间，恰好用一个空格隔开。即使某一行没有任何数需要输出，你也应输出一个空行。
请阅读样例来更好地理解这个格式。

Analysis

考试的时候是懵逼的，打了q=0的分就走了
经典的合并果子模型。开三个队列分别保存原始蚯蚓、切出来的前半段、切出来的后半段。同时增加可以变成单独减少保持相对的大小关系。每次从队头取三个当中的最长蚯蚓切，再放进队尾
复杂度玄学，有点不明白的

Code

/*
ID:wjp13241
PROG:earthworm
LANG:C++
*/
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <ctime>
#include <algorithm>
#include <iostream>
#include <vector>
#include <cmath>
#include <stack>
#include <queue>
#include <map>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define dfo(i,a,b) for(int i=a;i>=b;i--)
#define fore(i,x,e) for(int i=ls[x];i;i=e[i].next)
#define fil(x,t) memset(x,t,sizeof(x))
#define FILEIN(s) freopen(s,"r",stdin)
#define FILEOUT(s) freopen(s,"w",stdout)
#define STP system("pause")
#define min(x,y) x<y?x:y
#define max(x,y) x>y?x:y
#define MP(x,y) make_pair(x,y)
#define PuB(v,x) v.push_back(x)
#define PoB(v) v.pop_back()
#define ld long double
#define ll long long
#define db double
#define INF 0x3f3f3f3f
#define LIM 100000000
#define EPS 1e-4
#define N 200201
#define E N*20+1
#define L 16
using namespace std;
int a1[N*70],a2[N*70],a3[N*70];
int h1=1,t1=0,h2=1,t2=0,h3=1,t3=0;
__attribute__((optimize("O2")))
inline int read(){
char ch=getchar();int x=0,v=1;
while (ch>'9'||ch<'0'){if (ch=='-')v=-1;ch=getchar();}
while (ch<='9'&&ch>='0'){x=x*10+ch-'0';ch=getchar();}
return x*v;
}
__attribute__((optimize("O2")))
inline void write(int x){
char ch[L]={};int i=0;
if (x<0) putchar('-');
do{ch[++i]='0'+x%10;}while (x/=10);
while (i)putchar(ch[i--]);
}
__attribute__((optimize("O2")))
inline int cmp(int x,int y){return x>y;}
__attribute__((optimize("O2")))
inline int swap(int &x,int &y){x^=y;y^=x;x^=y;}
__attribute__((optimize("O2")))
inline int insert(int x,db p,int q,int i)
{
a2[++t2]=((int)floor(x*p))-i*q;
a3[++t3]=(x-(int)floor(x*p))-i*q;
}
__attribute__((optimize("O2")))
int main()
{
// FILEIN("earthworm.in");
// FILEOUT("earthworm.out");
int n=read(),m=read(),q=read(),u=read(),v=read(),t=read();
db p=u/(v*1.0);
fo(i,1,n)
a1[++t1]=read();
sort(a1+1,a1+t1+1,cmp);
fo(i,1,m)
{
int a=a1[h1],b=a2[h2],c=a3[h3];
if (h1>t1) a=-INF;
if (h2>t2) b=-INF;
if (h3>t3) c=-INF;
int now;
if (a>=b&&a>=c)
{
h1++;
now=a;
}
else if (b>=a&&b>=c)
{
h2++;
now=b;
}
else if (c>=a&&c>=b)
{
h3++;
now=c;
}
insert(now+i*q-q,p,q,i);
if (i%t==0)
{
write(now+i*q-q);
if (i/t<m/t)
putchar(' ');
}
}
putchar('n');
fo(i,1,n+m)
{
int a=a1[h1],b=a2[h2],c=a3[h3];
if (h1>t1) a=-INF;
if (h2>t2) b=-INF;
if (h3>t3) c=-INF;
int now;
if (a>=b&&a>=c)
{
h1++;
now=a;
}
else if (b>=a&&b>=c)
{
h2++;
now=b;
}
else if (c>=a&&c>=b)
{
h3++;
now=c;
}
if (i%t==0)
{
write(now+m*q);
if (i!=n+m)
putchar(' ');
}
}
putchar('n');
return 0;
}