Luogu P2827 蚯蚓

看到題目就能夠想到直接開的堆模擬的過程了吧,這個仍是很naive的git

注意在用堆作的時候也是要明智一點的,對於蚯蚓長度的相加確定不能直接遍歷並加上,仍是能夠差分一下的spa

其實說白了就是把集體加->單體減的一個小技巧,仍是挺經常使用的。code

而後看這個數據範圍猜測應該是有什麼\(O(n)\)的作法的,而後這就要發現題目中隱含的單調性隊列

咱們考慮講全部的蚯蚓分個類,全部初始時沒切割過的蚯蚓分爲一類,每次切割產生的較長的蚯蚓分爲一類,而產生的較短的蚯蚓分爲一類it

而後咱們推到一下就能夠發現,對於後面的兩類蚯蚓,它們知足單調性io

由於咱們根據切割的過程能夠發現:class

  • 先切割的蚯蚓長度必定比後切割的蚯蚓長度長
  • 同一種切法,後切割的必定比先切割的短

因此咱們再對初始的蚯蚓長度拍個序,就能夠獲得三個單調的隊列(注意不是單調隊列),每一次比較時咱們取出隊首並切割最長的一隻再丟進第二隊,第三隊便可遍歷

差分的思想仍是要的,不過這裏直接記錄了每一隻蚯蚓進隊的時間,每次注意加上增加的長度技巧

細節比較多,最後對這三隊蚯蚓作一次歸併便可(由於都是有序的)數據

所以複雜度爲\(O(n\ logn+m)\)

CODE

#include<cstdio>
#include<cctype>
#include<algorithm>
using namespace std;
const int N=100005,M=7000005;
int n,m,q,u,v,t,a[N],que[3][M],num[3][M],H[3],T[3],cut[M],ans[N+M],temp[N+M],cnt,tot;
inline char tc(void)
{
    static char fl[100000],*A=fl,*B=fl;
    return A==B&&(B=(A=fl)+fread(fl,1,100000,stdin),A==B)?EOF:*A++;
}
inline void read(int &x)
{
    x=0; char ch; while (!isdigit(ch=tc()));
    while (x=(x<<3)+(x<<1)+ch-'0',isdigit(ch=tc()));
}
inline void write(int x)
{
    if (x>9) write(x/10);
    putchar(x%10+'0');
}
inline bool cmp(int x,int y)
{
    return x>y;
}
inline void swap(int &a,int &b)
{
    int t=a; a=b; b=t;
}
inline void merge(void)
{
    register int i=H[0],j=H[1];
    while (i<=T[0]&&j<=T[1])
    if (que[0][i]+(m-num[0][i])*q>que[1][j]+(m-num[1][j])*q) temp[++cnt]=que[0][i]+(m-num[0][i])*q,++i;
    else temp[++cnt]=que[1][j]+(m-num[1][j])*q,++j;
    for (;i<=T[0];++i) temp[++cnt]=que[0][i]+(m-num[0][i])*q;
    for (;j<=T[1];++j) temp[++cnt]=que[1][j]+(m-num[1][j])*q;
    i=1; j=H[2];
    while (i<=cnt&&j<=T[2])
    if (temp[i]>que[2][j]+(m-num[2][j])*q) ans[++tot]=temp[i],++i; else ans[++tot]=que[2][j]+(m-num[2][j])*q,++j;
    for (;i<=cnt;++i) ans[++tot]=temp[i];
    for (;j<=T[2];++j) ans[++tot]=que[2][j]+(m-num[2][j])*q;
}
inline void print(void)
{
    for (register int i=t;i<=m;i+=t)
    write(cut[i]),putchar(' '); putchar('\n');
    for (register int i=t;i<=n+m;i+=t)
    write(ans[i]),putchar(' ');
}
int main()
{
    //freopen("CODE.in","r",stdin); freopen("CODE.out","w",stdout);
    register int i,j; read(n); read(m); read(q); read(u); read(v); read(t);
    for (i=1;i<=n;++i) read(a[i]); sort(a+1,a+n+1,cmp);
    for (i=1;i<=n;++i) que[0][++T[0]]=a[i]; H[0]=H[1]=H[2]=1;
    for (i=1;i<=m;++i)
    {
        int len=-1,id;
        for (j=0;j<3;++j) 
        if (H[j]<=T[j]) if (que[j][H[j]]+(i-num[j][H[j]]-1)*q>len) len=que[j][H[j]]+(i-num[j][H[j]]-1)*q,id=j;
        cut[i]=len; ++H[id]; int x=1LL*len*u/v,y=len-x; if (x<y) swap(x,y);
        que[1][++T[1]]=x; num[1][T[1]]=i; que[2][++T[2]]=y; num[2][T[2]]=i;
    }
    merge(); print(); return 0;
}
相關文章
相關標籤/搜索