NOIP2016全国信息学分区联赛提高组第二试 蚯蚓 ssl 2599 队列

题目描述

本题中，我们将用符号[c]表示对c向下取整，例如：[3.0」= [3.1」=[3.9」=3。

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

蛐蛐国里现在共有n只蚯蚓（n为正整数)。每只蚯蚓拥有长度，我们设第i只蚯蚓的长度为a_i(i=1,2,…,n)，并保证所有的长度都是非负整数（即:可能存在长度为0的蚯蚓）。

每一秒，神刀手会在所有的蚯蚓中，准确地找到最长的那一只（如有多个则任选一个）将其切成两半。神刀手切开蚯蚓的位置由常数p(是满足0

输入格式：

第一行包含六个整数n,m,q,u,v,t，其中：n,m,q的意义见【问题描述】；u,v,t均为正整数；你需要自己计算p=u/v(保证0

输出格式：

第一行输出[m/t]个整数，按时间顺序，依次输出第t秒，第2t秒，第3t秒……被切断蚯蚓（在被切断前）的长度。

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

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

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

分析

这题打暴力，用堆来做，理论上是可以拿80分，所以比赛还是打暴力好。
但是，正解又短又机智。
我们假设x1>x2。
那么x1肯定会在x2之前就被剪断。
设x1被剪成p1和q1，p1=⌊px1⌋,q1=x1−⌊px1⌋
设x2被剪成p2和q2，p2=⌊px2⌋,q2=x2−⌊px2⌋
我们可以得到一个很显然的结论，p1>p2，q1>q2，也就是说p1一定会在p2之前被剪断，q1一定会在q2之前被剪断。
我们来简单的证明一下：
假设x1剪断之后到剪断x2经过了i的时间，
那么此时x2=x2+i∗q,p1=⌊px1⌋+i∗q，q1=x1−⌊px1⌋+i∗q
剪断x2之后p2=⌊p(x2+i∗q)⌋,q2=x2+i∗q−⌊p(x2+i∗q)⌋
我们先比较一下p1和p2的大小关系：
我们如果不考虑取整的话，p1−p2=p(x1−x2)+q∗(i−i∗p)，因为p是小于1的，所式子一定是大于0的。
p2和q2的大小关系同理。
得证。
那么知道了这个结论之后要怎么去用呢。
我们设三个队列，一个a[0],a[1],a[2]。
a[0]预先存原序列的从大到小的排序。
a[1]存px的从大到小排序
a[2]存x-px的从大到小排序
那么每次取出最大的队首，然后拆掉存到a[1]和a[2]中，根据结论这样一定合法。

来自Facico

ps：我们题库能对，但是洛谷会超时两个点（很迷）

code

#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<string>
#include<algorithm>

using namespace std;

bool cmp(int x,int y)
{
    return x>y;
}

long long a[3][710];
int  n,m,q,u,v,t;
int head[3],tail[3];
int main()
{
    scanf("%d%d%d%d%d%d",&n,&m,&q,&u,&v,&t);
    for (int i=1;i<=n+m;i++)
        a[0][i]=a[1][i]=a[2][i]=-2000000000;
    for (int i=1;i<=n;i++)
        scanf("%lld",&a[0][i]);
    sort(a[0]+1,a[0]+1+n,cmp);
    tail[0]=n;
    tail[1]=tail[2]=0;
    head[1]=head[2]=head[0]=1;
    for (int i=1;i<=m;i++)
    {
        long long max=-2000000000;
        int k=0;
        for (int j=0;j<=2;j++)
            if (a[j][head[j]]>max)
            {
                max=a[j][head[j]];
                k=j;
            }
        long long tt=0;
        tt=max+(i-1)*q;
        if (i%t==0) printf("%lld ",tt);
        head[k]++;
        long long x,y;
        x=tt*u/v;
        y=tt-x;
        a[1][++tail[1]]=x-i*q;
        a[2][++tail[2]]=y-i*q;
    }
    printf("\n");
    for (int i=1;i<=m+n;i++)
    {
        long long max=-2000000000;
        int k=0;
        for (int j=0;j<=2;j++)
            if (a[j][head[j]]>max)
            {
                max=a[j][head[j]];
                k=j;
            }
        long long tt=0;
        tt=max+m*q;
        if(i%t==0) printf("%lld ",tt);
        head[k]++;
    }

}