KSGT

来源:互联网 发布:id是什么软件 编辑:程序博客网 时间:2024/05/21 10:57

上午,强开女神,瘫软在地

下午,怒写体积交

晚上,试水匈牙利


区间合并的问题,只是要维护两组的变量,这两组的变量有优先级的差异,才用的是两颗线段树各维护一组变量的技巧,减少了代码规模,其实想写一个基类,两个子类,继承重载的方法

注意!
1.在区间合并中,用query把长度转换成具体的[l, r]区间,可以大大较少update里面的思维复杂度
2.push_down 与 push_up只有全空与全满的情况,因为在几种update中只有全空、全满的更新,没有其他的
(正是因为陷入了半满的怪圈,才导致强上女神而不得)

#define _CRT_SECURE_NO_WARNINGS#include <iostream>#include <cstdio>#include <cstring>#include <algorithm>#include <vector>using namespace std;const int maxn = 100000 + 10;struct SegTree{    int l, r;    int left_time, max_time, right_time;    void init(int l, int r)    {        this->l = l, this->r = r;        left_time = max_time = right_time = r - l + 1;    }}goddess_segTree[maxn << 2], friend_segTree[maxn << 2];int t, n;void build(int k, int l, int r){    goddess_segTree[k].init(l, r);    friend_segTree[k].init(l, r);    if (l == r) return;    int m = (l + r) >> 1;    build(k << 1, l, m);    build(k << 1 | 1, m + 1, r); }int query(int k, SegTree * segTree, int len){    if (segTree[k].left_time >= len) return segTree[k].l;    if (segTree[k << 1].max_time >= len) return query(k << 1, segTree, len);    if (segTree[k << 1].right_time + segTree[k << 1 | 1].left_time >= len)        return segTree[k << 1].r - segTree[k << 1].right_time + 1;    return query(k << 1 | 1, segTree, len);}void push_down(int k){    if (friend_segTree[k].l == friend_segTree[k].r) return;    if (friend_segTree[k].max_time == 0)    {        friend_segTree[k << 1].left_time = friend_segTree[k << 1].max_time            = friend_segTree[k << 1].right_time = 0;        friend_segTree[k << 1 | 1].left_time = friend_segTree[k << 1 | 1].max_time            = friend_segTree[k << 1 | 1].right_time = 0;    }    if (friend_segTree[k].max_time == friend_segTree[k].r - friend_segTree[k].l + 1)    {        friend_segTree[k << 1].left_time = friend_segTree[k << 1].max_time = friend_segTree[k << 1].right_time            = friend_segTree[k << 1].r - friend_segTree[k << 1].l + 1;        friend_segTree[k << 1 | 1].left_time = friend_segTree[k << 1 | 1].max_time = friend_segTree[k << 1 | 1].right_time            = friend_segTree[k << 1 | 1].r - friend_segTree[k << 1 | 1].l + 1;    }    //_________________________________________________________________________    if (goddess_segTree[k].max_time == 0)    {        goddess_segTree[k << 1].left_time = goddess_segTree[k << 1].max_time            = goddess_segTree[k << 1].right_time = 0;        goddess_segTree[k << 1 | 1].left_time = goddess_segTree[k << 1 | 1].max_time            = goddess_segTree[k << 1 | 1].right_time = 0;    }    if (goddess_segTree[k].max_time == goddess_segTree[k].r - goddess_segTree[k].l + 1)    {        goddess_segTree[k << 1].left_time = goddess_segTree[k << 1].max_time = goddess_segTree[k << 1].right_time            = goddess_segTree[k << 1].r - goddess_segTree[k << 1].l + 1;        goddess_segTree[k << 1 | 1].left_time = goddess_segTree[k << 1 | 1].max_time = goddess_segTree[k << 1 | 1].right_time            = goddess_segTree[k << 1 | 1].r - goddess_segTree[k << 1 | 1].l + 1;    }}void push_up(int k){    if (friend_segTree[k].l == friend_segTree[k].r) return;    friend_segTree[k].left_time = friend_segTree[k << 1].left_time;    if (friend_segTree[k << 1].l + friend_segTree[k << 1].left_time - 1 == friend_segTree[k << 1].r)        friend_segTree[k].left_time += friend_segTree[k << 1 | 1].left_time;    friend_segTree[k].right_time = friend_segTree[k << 1 | 1].right_time;    if (friend_segTree[k << 1 | 1].r - friend_segTree[k << 1 | 1].right_time + 1 == friend_segTree[k << 1 | 1].l)        friend_segTree[k].right_time += friend_segTree[k << 1].right_time;    friend_segTree[k].max_time = max(friend_segTree[k << 1].max_time, friend_segTree[k << 1 | 1].max_time);    friend_segTree[k].max_time = max(friend_segTree[k].max_time,         friend_segTree[k << 1].right_time + friend_segTree[k << 1 | 1].left_time);    //_________________________________________________________________________    goddess_segTree[k].left_time = goddess_segTree[k << 1].left_time;    if (goddess_segTree[k << 1].l + goddess_segTree[k << 1].left_time - 1 == goddess_segTree[k << 1].r)        goddess_segTree[k].left_time += goddess_segTree[k << 1 | 1].left_time;    goddess_segTree[k].right_time = goddess_segTree[k << 1 | 1].right_time;    if (goddess_segTree[k << 1 | 1].r - goddess_segTree[k << 1 | 1].right_time + 1 == goddess_segTree[k << 1 | 1].l)        goddess_segTree[k].right_time += goddess_segTree[k << 1].right_time;    goddess_segTree[k].max_time = max(goddess_segTree[k << 1].max_time, goddess_segTree[k << 1 | 1].max_time);    goddess_segTree[k].max_time = max(goddess_segTree[k].max_time,        goddess_segTree[k << 1].right_time + goddess_segTree[k << 1 | 1].left_time);}void update(int k, int l, int r, int flag){    if (l > friend_segTree[k].r || r < friend_segTree[k].l) return;    if (l <= friend_segTree[k].l && friend_segTree[k].r <= r)    {        friend_segTree[k].left_time = friend_segTree[k].max_time            = friend_segTree[k].right_time = 0;        if (flag)        {            goddess_segTree[k].left_time = goddess_segTree[k].max_time                = goddess_segTree[k].right_time = 0;        }        return;    }    push_down(k);    update(k << 1, l, r, flag);    update(k << 1 | 1, l, r, flag);    push_up(k);}void solve_ds(int len){    if (friend_segTree[1].max_time >= len)    {        int beg = query(1, friend_segTree, len);        printf("%d,let's fly\n", beg);        update(1, beg, beg + len - 1, 0);        return;    }    puts("fly with yourself");}void solve_ns(int len){    if (friend_segTree[1].max_time >= len)    {        int beg = query(1, friend_segTree, len);        printf("%d,don't put my gezi\n", beg);        update(1, beg, beg + len - 1, 1);        return;    }    if (goddess_segTree[1].max_time >= len)    {        int beg = query(1, goddess_segTree, len);        printf("%d,don't put my gezi\n", beg);        update(1, beg, beg + len - 1, 1);        return;    }    puts("wait for me");}void update_study(int k, int l, int r){    if (l > friend_segTree[k].r || r < friend_segTree[k].l) return;    if (l <= friend_segTree[k].l && friend_segTree[k].r <= r)    {        friend_segTree[k].left_time = friend_segTree[k].max_time = friend_segTree[k].right_time            = friend_segTree[k].r - friend_segTree[k].l + 1;        goddess_segTree[k].left_time = goddess_segTree[k].max_time = goddess_segTree[k].right_time            = goddess_segTree[k].r - goddess_segTree[k].l + 1;        return;    }    push_down(k);    update_study(k << 1, l, r);    update_study(k << 1 | 1, l, r);    push_up(k);}void solve_study(int l, int r){    puts("I am the hope of chinese chengxuyuan!!");    update_study(1, l, r);}int main(){    int ttt; scanf("%d", &ttt);    int kase = 0;    while (ttt--)    {        printf("Case %d:\n", ++kase);        scanf("%d%d", &t, &n);        build(1, 1, t);        for (int i = 0; i < n; i++)        {            char cmd[10];            int a, b;            scanf("%s%d", cmd, &a);            if (cmd[0] == 'D') solve_ds(a);            else if (cmd[0] == 'N') solve_ns(a);            else if (cmd[0] == 'S')            {                scanf("%d", &b);                solve_study(a, b);            }        }    }    return 0;}

这种代码仿佛是在写项目….


周长交 -> 面积交 -> 体积交
一维 二维 三维

一维的周长中,差值来计算

二维的面积中,pre_len * deta_h

三维的体积中,以体积层为枚举单位,利用当前层内的面积 * deta_z

主要难点在于枚举体积层,其他的,比如重合三次,因为刚刚完成了面积交,所以很顺畅

#define _CRT_SECURE_NO_WARNINGS#include <iostream>#include <cstdio>#include <cstring>#include <vector>#include <algorithm>using namespace std;typedef long long ll;const int maxn = 3000 + 10;struct Line{    ll l, r, h;    ll z1, z2;    ll flag;    Line(ll l = 0, ll r = 0, ll h = 0, ll z1 = 0, ll z2 = 0, ll flag = 0)    {        this->l = l, this->r = r, this->h = h, this->z1 = z1, this->z2 = z2, this->flag = flag;    }    bool operator < (const Line & rhs) const    {        return h < rhs.h;    }}line[maxn];struct SegTree{    ll l, r;    ll cnt;    ll sum1, sum2, sum3;    void init(ll l, ll r)    {        this->l = l, this->r = r;        sum1 = sum2 = sum3 = 0;        cnt = 0;    }}segTree[maxn << 2];ll new_x[maxn], new_z[maxn];ll cur1;void build(int k, ll l, ll r){    segTree[k].init(l, r);    if (l + 1 >= r) return;    ll m = (l + r) >> 1;    build(k << 1, l, m);    build(k << 1 | 1, m, r);}ll depress_x(){    vector<ll> pool;    for (int i = 0; i < cur1; i++)    {        pool.push_back(line[i].l);        pool.push_back(line[i].r);    }    sort(pool.begin(), pool.end());    pool.erase(unique(pool.begin(), pool.end()), pool.end());    for (int i = 0; i < cur1; i++)    {        ll cur;        cur = find(pool.begin(), pool.end(), line[i].l) - pool.begin();        new_x[cur] = line[i].l; line[i].l = cur;        cur = find(pool.begin(), pool.end(), line[i].r) - pool.begin();        new_x[cur] = line[i].r; line[i].r = cur;    }    return pool.size() - 1;}ll depress_z(){    vector<ll> pool;    for (int i = 0; i < cur1; i++)    {        pool.push_back(line[i].z1);        pool.push_back(line[i].z2);    }    sort(pool.begin(), pool.end());    pool.erase(unique(pool.begin(), pool.end()), pool.end());    for (int i = 0; i < cur1; i++)    {        ll cur;        cur = find(pool.begin(), pool.end(), line[i].z1) - pool.begin();        new_z[cur] = line[i].z1;        cur = find(pool.begin(), pool.end(), line[i].z2) - pool.begin();        new_z[cur] = line[i].z2;    }    return pool.size() - 1;}void chan_sum(int k){    if (segTree[k].cnt >= 3)    {        segTree[k].sum3 = new_x[segTree[k].r] - new_x[segTree[k].l];        segTree[k].sum2 = new_x[segTree[k].r] - new_x[segTree[k].l];        segTree[k].sum1 = new_x[segTree[k].r] - new_x[segTree[k].l];    }    else if (segTree[k].cnt == 2)    {        segTree[k].sum3 = segTree[k << 1].sum1 + segTree[k << 1 | 1].sum1;        segTree[k].sum2 = new_x[segTree[k].r] - new_x[segTree[k].l];        segTree[k].sum1 = new_x[segTree[k].r] - new_x[segTree[k].l];    }    else if (segTree[k].cnt == 1)    {        segTree[k].sum3 = segTree[k << 1].sum2 + segTree[k << 1 | 1].sum2;        segTree[k].sum2 = segTree[k << 1].sum1 + segTree[k << 1 | 1].sum1;        segTree[k].sum1 = new_x[segTree[k].r] - new_x[segTree[k].l];    }    else if (segTree[k].cnt == 0)    {        segTree[k].sum3 = segTree[k << 1].sum3 + segTree[k << 1 | 1].sum3;        segTree[k].sum2 = segTree[k << 1].sum2 + segTree[k << 1 | 1].sum2;        segTree[k].sum1 = segTree[k << 1].sum1 + segTree[k << 1 | 1].sum1;    }}void update(int k, ll l, ll r, ll c){    if (l > segTree[k].r || r < segTree[k].l) return;    if (l <= segTree[k].l && segTree[k].r <= r)    {        segTree[k].cnt += c;        chan_sum(k);        return;    }    if (segTree[k].l + 1 == segTree[k].r)    {        if (segTree[k].cnt >= 3)        {            segTree[k].sum3 = new_x[segTree[k].r] - new_x[segTree[k].l];            segTree[k].sum2 = new_x[segTree[k].r] - new_x[segTree[k].l];            segTree[k].sum1 = new_x[segTree[k].r] - new_x[segTree[k].l];        }        else if (segTree[k].cnt == 2)        {            segTree[k].sum3 = 0;            segTree[k].sum2 = new_x[segTree[k].r] - new_x[segTree[k].l];            segTree[k].sum1 = new_x[segTree[k].r] - new_x[segTree[k].l];        }        else if (segTree[k].cnt == 1)        {            segTree[k].sum3 = 0;            segTree[k].sum2 = 0;            segTree[k].sum1 = new_x[segTree[k].r] - new_x[segTree[k].l];        }        else if (segTree[k].cnt == 0)        {            segTree[k].sum3 = 0;            segTree[k].sum2 = 0;            segTree[k].sum1 = 0;        }        return;    }    update(k << 1, l, r, c);    update(k << 1 | 1, l, r, c);    chan_sum(k);}ll square(ll down_lim, ll up_lim){    ll pre_h = 0, pre_len = 0;    ll res = 0;    for (int i = 0; i < cur1; i++)    {        if (line[i].z1 <= down_lim && up_lim <= line[i].z2)        {            update(1, line[i].l, line[i].r, line[i].flag);            res += (line[i].h - pre_h) * pre_len;            pre_h = line[i].h;            pre_len = segTree[1].sum3;        }    }    return res;}ll solve(){    ll w_x = depress_x();    ll w_z = depress_z();    ll res = 0;    for (int i = 1; i <= w_z; i++)    {        build(1, 0, w_x);        ll down_lim = new_z[i - 1];        ll up_lim = new_z[i];        res += (up_lim - down_lim) * square(down_lim, up_lim);    }    return res;}int main(){    int t; scanf("%d", &t);    int kase = 0;    while (t--)    {        cur1 = 0;        int n; scanf("%d", &n);        for (int i = 0; i < n; i++)        {            ll xx1, yy1, zz1, xx2, yy2, zz2;            scanf("%lld%lld%lld%lld%lld%lld", &xx1, &yy1, &zz1, &xx2, &yy2, &zz2);            line[cur1++] = Line(xx1, xx2, yy1, zz1, zz2, 1);            line[cur1++] = Line(xx1, xx2, yy2, zz1, zz2, -1);        }        sort(line, line + cur1);        printf("Case %d: %lld\n", ++kase, solve());    }    return 0;}

复习了匈牙利算法,难为我还能记得,然而建模上面缺乏经验,被队友随意摩擦

原创粉丝点击