2019 HDOJ Multi-University Training Contest Stage 2(杭电多校)

服务器时不时爆炸，有点难受。

题目链接：http://acm.hdu.edu.cn/userloginex.php?cid=849

---

A：

神仙题。不可做题。

B：

dp。

C：

推式子题。

D：

边分治。

E：

可以数学推理的题。但是显然打表更快找出规律。对打出来的结果做两次差分即可。

/* basic header */
#include <bits/stdc++.h>
/* define */
#define ll long long
#define dou double
#define pb emplace_back
#define mp make_pair
#define sot(a,b) sort(a+1,a+1+b)
#define rep1(i,a,b) for(int i=a;i<=b;++i)
#define rep0(i,a,b) for(int i=a;i<b;++i)
#define eps 1e-8
#define int_inf 0x3f3f3f3f
#define ll_inf 0x7f7f7f7f7f7f7f7f
#define lson (curpos<<1)
#define rson (curpos<<1|1)
#define mid (curl+curr>>1)
/* namespace */
using namespace std;
/* header end */

const ll mod = 998244353;

ll qp(ll x, ll b) {
    ll ret = 1, base = x;
    while (b) {
        if (b & 1) ret = ret * base % mod;
        base = base * base % mod;
        b >>= 1;
    }
    return ret;
}

int main() {
    ll n;
    while (~scanf("%lld", &n)) {
        ll k = 0, tmp = 3;
        rep1(i, 2, n) k += tmp, tmp += 2;
        printf("%lld
", k * qp(9, mod - 2) % mod);
    }
    return 0;
}

F：

fwt题。

G：

不均等博弈题。用surrealnumber秒杀。

H：

最小割。

I：

求出本质不同的回文串的数量分布(求每种回文串的个数)，然后对每种check一下叠加答案。manacher或者字符串hash都可以。

J：

签到题。1e6+3的阶乘后面全是0已经暗示一切。

/* basic header */
#include <bits/stdc++.h>
/* define */
#define ll long long
#define dou double
#define pb emplace_back
#define mp make_pair
#define sot(a,b) sort(a+1,a+1+b)
#define rep1(i,a,b) for(int i=a;i<=b;++i)
#define rep0(i,a,b) for(int i=a;i<b;++i)
#define eps 1e-8
#define int_inf 0x3f3f3f3f
#define ll_inf 0x7f7f7f7f7f7f7f7f
#define lson (curpos<<1)
#define rson (curpos<<1|1)
#define mid (curl+curr>>1)
/* namespace */
using namespace std;
/* header end */

const int mod = 1e6 + 3;
ll p[mod + 10];
int n;

int main() {
    p[0] = p[1] = 1;
    rep0(i, 2, mod) p[i] = 1ll * i * p[i - 1] % mod;
    while (~scanf("%d", &n)) {
        if (n >= mod) {
            puts("0");
            continue;
        }
        n %= mod;
        printf("%lld
", p[n]);
    }
    return 0;
}

K：

对于某个区间，首先考虑区间第1、2、3大的边能否构成三角形，然后考虑区间第2、3 、4大的边能否构成三角形，以此类推。

只要考虑区间前44大即可，因为区间最坏情况是斐波那契数列，然而第45项就爆了1e9的范围。时间复杂度O(44*nlogn)。

/* basic header */
#include <iostream>
#include <cstring>
/* define */
#define ll long long
#define dou double
#define pb emplace_back
#define mp make_pair
#define sot(a,b) sort(a+1,a+1+b)
#define rep1(i,a,b) for(int i=a;i<=b;++i)
#define rep0(i,a,b) for(int i=a;i<b;++i)
#define eps 1e-8
#define int_inf 0x3f3f3f3f
#define ll_inf 0x7f7f7f7f7f7f7f7f
#define lson (curpos<<1)
#define rson (curpos<<1|1)
/* namespace */
using namespace std;
/* header end */

const int maxn = 1e5 + 10;
int a[maxn];

struct Node {
    int l, r, val[50];
};

Node segt[maxn << 2];

void maintain(Node &fa, Node &ls, Node &rs) {
    int i = 0, j = 0;
    rep0(c, 0, 50) {
        if (i >= 50)
            fa.val[c] = rs.val[j++];
        else if (j >= 50)
            fa.val[c] = ls.val[i++];
        else if (ls.val[i] < rs.val[j])
            fa.val[c] = rs.val[j++];
        else fa.val[c] = ls.val[i++];
    }
}

void build(int curpos, int curl, int curr) {
    segt[curpos].l = curl, segt[curpos].r = curr;
    // memset(segt[curpos].val, 0, sizeof(segt[curpos].val));
    if (curl < curr) { // if is not leaf node
        int mid = curl + curr >> 1;
        build(lson, curl, mid); build(rson, mid + 1, curr);
        maintain(segt[curpos], segt[lson], segt[rson]);
    } else
        segt[curpos].val[0] = a[curl]; // if is leaf node
}

Node query(int curpos, int curl, int curr) {
    if (segt[curpos].l == curl && segt[curpos].r == curr) return segt[curpos];
    else {
        int mid = segt[curpos].l + segt[curpos].r >> 1;
        if (curr <= mid) return query(lson, curl, curr);
        else if (curl > mid) return query(rson, curl, curr);
        else {
            Node lnode = query(lson, curl, mid), rnode = query(rson, mid + 1, curr);
            Node ret; ret.l = curl, ret.r = curr;
            maintain(ret, lnode, rnode);
            return ret;
        }
    }
}

int main() {
    int n, q;
    while (~scanf("%d%d", &n, &q)) {
        rep1(i, 1, n) scanf("%d", &a[i]);
        build(1, 1, n);
        while (q--) {
            int l, r; scanf("%d%d", &l, &r);
            Node cur = query(1, l, r);
            ll ans = -1;
            rep0(i, 0, 48) {
                if ((ll)cur.val[i] < (ll)cur.val[i + 1] + (ll)cur.val[i + 2]) {
                    ans = (ll)cur.val[i] + (ll)cur.val[i + 1] + (ll)cur.val[i + 2];
                    break;
                }
            }
            printf("%lld
", ans);
        }
    }
    return 0;
}

L：

又是线段树题，背锅。