作るRounddogハッピー(2019鮎杭州電動マルチ第十校フィールド1011 + HDU6701 +ヒューリスティックパーティション)

トピックリンク

ポータル

問題の意味

サブインターバルの数を満たすために探している\(a_l、A_ {L + 1}、\ドットは、A_R \)は同じではないa_l、A_ {(\(MAX (R-L - L + 1}、\ドット、A_R) 1)\のLeq K \)

考え

私はヒューリスティックパーティションを聞いた後、次のボードを設定する方法を学ぶために行った、地元のねじれが一晩のゲームを変更した後、ゲームをオーバーサンプリングしすることはできません書きました。

問題のこのパーティションのヒューリスティックアイデアは、プロセスであるに見えます\([L、R]を\ ) 位置の間隔の最大値を見つけるために\(MID \)を、左のセクションの長さより短くまたは短い右部分の半分、そして暴力を見て、短い統計の貢献の一環。

前処理された第一に(私は\)\左点として、右の点は、間隔内の同じ数ではない\(R [I] \)とで\(Iは\)間隔の右エンドポイントの左の点は同じ数ではない\ (L [I] \) 次いで\(ST \)位置の間隔最大値を処理するテーブル。

我々は最大の位置を知っているパーティションの時には、そのときカウント拠出ため(例として、左側の部分を取る)不等式の列挙\(Lの\) そして最近では、タイトルの意味を満たすために適切な場所を転置演算し\ (R&LT \ GEQ [MID]。1 -K-L + \) 次いでため\(Lの\)左右点の範囲内の点の数が同じではない\(R&LT [L] \)最も満足のタイトルを意味します(右端点)、この値を減算する寄与する左点の複雑さである\(O(nlog(N))を\)

コード

#include <set>
#include <map>
#include <deque>
#include <queue>
#include <stack>
#include <cmath>
#include <ctime>
#include <bitset>
#include <cstdio>
#include <string>
#include <vector>
#include <cassert>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;

typedef long long LL;
typedef pair<LL, LL> pLL;
typedef pair<LL, int> pLi;
typedef pair<int, LL> pil;;
typedef pair<int, int> pii;
typedef unsigned long long uLL;

#define lson (rt<<1),L,mid
#define rson (rt<<1|1),mid + 1,R
#define lowbit(x) x&(-x)
#define name2str(name) (#name)
#define bug printf("*********\n")
#define debug(x) cout<<#x"=["<<x<<"]" <<endl
#define FIN freopen("/home/dillonh/CLionProjects/Dillonh/in.txt","r",stdin)
#define IO ios::sync_with_stdio(false),cin.tie(0)

const double eps = 1e-8;
const int mod = 1000000007;
const int maxn = 300000 + 7;
const double pi = acos(-1);
const int inf = 0x3f3f3f3f;
const LL INF = 0x3f3f3f3f3f3f3f3fLL;

LL ans;
int t, n, K;
int a[maxn], L[maxn], R[maxn], vis[maxn], dp[maxn][20], pos[maxn][20];

void init() {
    for(int j = 1; j < 20; ++j) {
        if(1<<(j-1) > n) break;
        for(int i = 1; i + (1<<j) - 1 <= n; ++i) {
            if(dp[i][j-1] >= dp[i+(1<<(j-1))][j-1]) {
                dp[i][j] = dp[i][j-1];
                pos[i][j] = pos[i][j-1];
            } else {
                dp[i][j] = dp[i+(1<<(j-1))][j-1];
                pos[i][j] = pos[i+(1<<j-1)][j-1];
            }
        }
    }
}

int query(int l, int r) {
    int k = log(r - l + 1) / log(2);
    if(dp[l][k] >= dp[r-(1<<k)+1][k]) return pos[l][k];
    else return pos[r-(1<<k)+1][k];
}

void solve(int l, int r) {
    if(l > r) return;
    int mid = query(l, r);
    if(r - mid > mid - l) {
        for(int i = l; i <= mid; ++i) {
            int rs = a[mid] - K + i - 1;
            rs = max(rs, mid);
            int dd = min(r, R[i]);
            if(rs > dd) continue;
            ans += dd - rs + 1;
        }
    } else {
        for(int i = mid; i <= r; ++i) {
            int ls = K - a[mid] + i + 1;
            ls = min(ls, mid);
            int dd = max(l, L[i]);
            if(ls < dd) continue;
            ans += ls - dd + 1;
        }
    }
    solve(l, mid - 1), solve(mid + 1, r);
}

int main() {
#ifndef ONLINE_JUDGE
    FIN;
#endif
    scanf("%d", &t);
    while(t--) {
        scanf("%d%d", &n, &K);
        for(int i = 1; i <= n; ++i) {
            scanf("%d", &a[i]);
            dp[i][0] = a[i];
            pos[i][0] = i;
            vis[i] = 0;
        }
        int r = 2;
        vis[a[1]] = 1;
        for(int i = 1; i <= n; ++i) {
            while(r <= n && !vis[a[r]]) {
                vis[a[r]] = 1;
                ++r;
            }
            vis[a[i]] = 0;
            R[i] = r - 1;
        }
        vis[a[n]] = 1;
        int l = n - 1;
        for(int i = n; i >= 1; --i) {
            while(l >= 1 && !vis[a[l]]) {
                vis[a[l]] = 1;
                --l;
            }
            vis[a[i]] = 0;
            L[i] = l + 1;
        }
        init();
        ans = 0;
        solve(1, n);
        printf("%lld\n", ans);
    }
    return 0;
}

おすすめ

転載: www.cnblogs.com/Dillonh/p/11390927.html
おすすめ