蚯蚓

最新推荐文章于 2020-09-23 21:10:42 发布

原创最新推荐文章于 2020-09-23 21:10:42 发布 · 315 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#数据结构

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

蛐蛐国里现在共有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*106，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……的长度。

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

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

样例输入 [复制]
输入样例#1：
3 7 1 1 3 1
3 3 2

输入样例#2：
3 7 1 1 3 2
3 3 2

输入样例#3：
3 7 1 1 3 9
3 3 2
样例输出 [复制]
输出样例#1：
3 4 4 4 5 5 6
6 6 6 5 5 4 4 3 2 2
输出样例#2：
4 4 5
6 5 4 3 2
输出样例#3：
2

最开始直接用了小根堆。然后 t 了两个点。最后发现是单调的,但是就算是写成线性的。在学校的评测机还是卡了。luogu能过

关键点: 发现此题中隐含的单调性.

发现先被切掉的蚯蚓分成的蚯蚓一定比后切掉的蚯蚓分成的蚯蚓大. 　　假设这两只蚯蚓分别为a,ba,b,其中a>ba>b.那么它被切成a_1,a_2a
　　
　　那么就可以将这两堆依次存储, 加上还没被切过的蚯蚓.每次要切时在这三堆里面选择最大的, 切完再依次放回去. 　　所以这么做时间复杂度为O(m)O(m).再优化一下细节基本上就没问题了.

75—80分

#include<iostream>
#include<cstdio>
#include<queue>

using namespace std;
priority_queue<int>q;
int main()
{
    int n,m,qq,u,v,t;
//n:蚯蚓数 ; m : m秒内的状况 ; q ：其余蚯蚓增加的长度  
    int qwq; 
    scanf("%d%d%d%d%d%d",&n,&m,&qq,&u,&v,&t);
    double p;
    p = u * 1.0 / v;
    for(int i = 1;i <= n;++ i)
    {
        scanf("%d",&qwq);
        q.push(qwq);
    }
    int jss = 0;
    while(m --) 
    {
        int awa = q.top() + jss * qq;
        jss ++;
        if(jss % t == 0)
            printf("%d ",awa);
        q.pop();
        int X1 = awa * p;
        int aa = X1 - jss * qq;
        int bb = (awa - X1) - jss * qq;
        q.push(aa);q.push(bb);
    }
    printf("\n");
    int js = 0;
    while(!q.empty())
    {
        js ++;
        if(js % t == 0)
            printf("%d ",q.top() + jss * qq);
        q.pop();
    }
    printf("\n");
    return 0;
}

100分

在这里插入代码片#include<bits/stdc++.h>
using namespace std;
#define N 7000005
inline char nc(){
    static char buf[100000],*p1=buf,*p2=buf;
    return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline int read(){
    char ch=nc();int sum=0;
    while(!(ch>='0'&&ch<='9'))ch=nc();
    while(ch>='0'&&ch<='9')sum=sum*10+ch-48,ch=nc();
    return sum;
}
int n,m,q,u,v,t;
int cut1[N],cut2[N],h1,h2,h,num[N],t0,t1,t2,sigma=0,top=0;
double p;
priority_queue<int>ans;
bool cmp(int x,int y){
	return x>y;
}
signed main(){
    n=read(),m=read();q=read();u=read();v=read();t=read();
	p=(double)u/v;
	for(int i=1;i<=n;i++){
		num[i]=read();
	}
	h1=h2=h=1;t0=n;
	sort(num+1,num+n+1,cmp);
	for(int i=1;i<=m;i++){
		if(h>t0){
			if(cut1[h1]>cut2[h2]) top=cut1[h1++];
			else top=cut2[h2++]; 
		}
		else if(num[h]>=cut1[h1]&&num[h]>=cut2[h2]){
			top=num[h++];
		}
        else if(cut1[h1]>=cut2[h2]&&num[h]<=cut1[h1]){
		   top=cut1[h1++];
	    }
	    else (top=cut2[h2++]);
	    top+=sigma;
	    int a1=(int)(p*(double)top),a2=top-a1;
	    sigma+=q;
	    a1-=sigma;a2-=sigma;
	    cut1[++t1]=a1,cut2[++t2]=a2;
	    if(i%t==0) printf("%d ",top);
    }
    puts("");
    for(int i=h;i<=t0;i++) ans.push(num[i]);
    for(int i=h1;i<=t1;i++) ans.push(cut1[i]);
    for(int i=h2;i<=t2;i++) ans.push(cut2[i]);
    for(int i=1;ans.size();i++){
    	if(i%t==0) printf("%d ",ans.top()+sigma);
    	ans.pop();
	}
	return 0;
}