原来真的是按题意模拟啊,还以为有高能的算法可以直接算每个$t$的值。
考虑到先切的蚯蚓一定比后切的蚯蚓长,于是可以弄三个队列分别存放原来的序列和两个切开后的序列,每次取出三个队头的最大值进行扩展。
考虑到每秒钟除了取出来的队头其他的长度都会增加$q$,那么我们可以写一个全局变量$tag$标记现在进行了几轮,然后每一次进队的时候反向减去这一个$tag$就好了。
时间复杂度$O(nlogn)$。
$stl$的$queue$开了$O2$之后就很快了
Code:
#include <cstdio> #include <cstring> #include <queue> #include <algorithm> using namespace std; typedef double db; const int N = 7e6 + 5; const int inf = 1 << 30; int n, m, q, t, a[N * 3]; db p; queue <int> Q[3]; inline void read(int &X) { X = 0; char ch = 0; int op = 1; for(; ch > '9'|| ch < '0'; ch = getchar()) if(ch == '-') op = -1; for(; ch >= '0' && ch <= '9'; ch = getchar()) X = (X << 3) + (X << 1) + ch - 48; X *= op; } bool cmp(const int x, const int y) { return x > y; } inline int bet(int x, int y, int rx, int ry) { return rx > ry ? x : y; } inline int getQ(int now) { if(Q[now].empty()) return -inf; else return Q[now].front(); } inline int getMax() { int res[3]; for(int i = 0; i < 3; i++) res[i] = getQ(i); return bet(bet(0, 1, res[0], res[1]), 2, res[bet(0, 1, res[0], res[1])], res[2]); } int main() { int u, v; read(n), read(m), read(q), read(u), read(v), read(t); p = (db)u / v; for(int i = 1; i <= n; i++) read(a[i]); sort(a + 1, a + 1 + n, cmp); for(int i = 1; i <= n; i++) Q[0].push(a[i]); int tag = 0; for(int i = 1; i <= m; i++) { int now = getMax(); int z = Q[now].front() + tag; Q[now].pop(); if(i % t == 0) printf("%d ", z); int x = (int)z * p, y = z - x; tag += q; x -= tag, y -= tag; Q[1].push(x), Q[2].push(y); } printf("\n"); int len = 0; for(int i = 0; i < 3; i++) for(;!Q[i].empty(); Q[i].pop()) a[++len] = Q[i].front(); sort(a + 1, a + 1 + len, cmp); /* for(int i = 1; i <= len; i++) printf("%d ", a[i]); printf("\n"); */ int cnt = (n + m) / t; for(int i = 1; i <= cnt; i++) printf("%d ", a[i * t] + tag); printf("\n"); return 0; }