1. 程式人生 > >NOIP2016 DAY2 T2蚯蚓

NOIP2016 DAY2 T2蚯蚓

fine 普通 day2 esp n) 一位 之間 整數 while

傳送門

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……的長度。 同一行中相鄰的兩個數之間,恰好用一個空格隔開。即使某一行沒有任何數需要 輸出,你也應輸出一個空行。 請閱讀樣例來更好地理解這個格式。

Sample Input

3 7 1 1 3 1
3 3 2

Sample Output

3 4 4 4 5 5 6
6 6 6 5 5 4 4 3 2 2
感覺挺莫名奇妙的..... 在洛谷上AC了這題,然後去下載官方數據, 居然超時?? 也許是我這臺機子太慢了..................... 正解: 第一遍看完題目時就覺得用優先隊列,但想想.... CCF會這麽容易給你滿分嗎?? 做夢吧。。。 優先隊列65分。 正解是用三個普通的隊列做的。 一開始將所有蚯蚓讀入,放進一個隊列並且排序,從大到小。 然後拿出最長的一條,進行切割,將 x*q 放在一個隊列,將x - x*q放在另一個隊列。 經過證明,這兩個隊列都是單調的。
利用這一點,我們只需每次從這三個隊列裏面挑出一個最大的蚯蚓進行切割然後放進那兩個隊列就行了。 下面貼代碼,有問題留言。
#include<cstdio>
#include<algorithm>
#define N 20000900
using namespace std;

int que[3][N],head[3],tail[3];

int cmp(const int &a,const int &b){
    return a > b;
}

int read(){
    int x = 0;
    char ch = getchar();
    while(ch < 0 || ch > 9)ch = getchar();
    while(ch >= 0 && ch <= 9){
        x = x * 10 + ch -0;
        ch = getchar();
    }
    return x;
}

int main(){
    int n = read(),m = read(),q = read(),u = read(),v = read(),t = read();
    double p = 1.0*u/v;
    for(int i = 1; i <= n; i++)que[0][i] = read();
    sort(que[0]+1,que[0]+n+1,cmp);
    int first = 1,len = 0;
    head[0] = head[1] = head[2] = 1;
    tail[0] = n,tail[1] = tail[2] = 0;
    for(int i = 1; i <= m; i++){
        long long mx = -99999999999,pos = 0;
        for(int j = 0; j <= 2; j++)
          if(head[j] <= tail[j] && que[j][head[j]] > mx){
                  mx = que[j][head[j]];
                  pos = j;
          }
        head[pos]++;
        mx += len;
        if(i % t == 0){
            if(first)first = 0;
            else printf(" ");
            printf("%d",mx);
        }
        int x1 = (int)(mx*u/v);
        int x2 = mx - x1;
        que[1][++tail[1]] = x1 - len-q;
        que[2][++tail[2]] = x2 - len-q;
        len += q;
        
    }
    printf("\n");
    first = 1;
    for(int i = 1; i <= n + m; i++){
        long long mx = -9999999999999,pos = 0;     //mx要設小一點,不然會爆
        for(int j = 0;  j <= 2; j++)
          if(head[j] <= tail[j] && que[j][head[j]] > mx){
               mx = que[j][head[j]];
               pos = j;
          }
        head[pos]++;
        mx += len;
        if(i % t == 0){
            if(first)first = 0;
            else printf(" ");
            printf("%lld",mx);
        }
    }
    printf("\n");
    return 0;
}

NOIP2016 DAY2 T2蚯蚓