uva 1618（枚举技巧）

最新推荐文章于 2021-04-14 20:23:29 发布

playwfun

最新推荐文章于 2021-04-14 20:23:29 发布

阅读量926

点赞数

分类专栏：高效

本文链接：https://blog.csdn.net/playwfun/article/details/51125144

版权

高效专栏收录该内容

9 篇文章 0 订阅

订阅专栏

本文介绍了一种解决特定条件下四数问题的高效算法。该算法利用离散化技术和范围查询优化（RMQ），实现时间复杂度为O(n^2)的解决方案。文章详细解释了如何通过两次区间查找找到符合条件的数，并通过预处理达到快速查询的目的。

摘要由CSDN通过智能技术生成

直接说解法：

我们只考虑1<=p < q < r < s <= n , Nq > Ns > Np > Nr 的这种情况，先枚举 p , q 的位置，那么先选s的位置，s的位置越靠后越好，这样会留下更大的空间选r

又知道所有数不同，离散化位1 到 n ，那么假设p , q 位置确定，以pos数组记录每个1 -- n 每个离散化后的数对应的数组位置，那么 s的最佳位置为 pos数组上

Np + 1 --> Nq - 1 这个区间的最大值记为 max_pos , 那么下面只需要寻找在数组a上的区间 q + 1 ---> max_pos - 1 是否有值比Np小就可以了。

由于这两次查找区间最值都可以通过RMQ预处理做到O（1）的查询，所以总时间复杂度位 O（n^2）。

#include <cstdio>
#include <algorithm>
#include <cstring>
#include <iostream>
#include <map>
#include <vector>
using namespace std;

typedef long long ll;
#define rep( i , n) for(int i = 0 ; i<(int)n;i++)
#define rep1(i ,x , y) for(int i = (int)x ; i<=(int)y;i++)
const int N = 5050;

int max_[N][20],min_[N][20],lg2[N]; //20不一定是唯一的。需要计算log(N)/log(2)
void ST(int *a,int n)
{
    lg2[0]=-1;
    for(int i=1; i<=n; i++)
        lg2[i]=lg2[i-1]+(i&(i-1)?0:1);
    for(int i=0; i<n; i++) max_[i][0]=a[i]; //a第一个数从零开始
    for(int j=1; j<=lg2[n]; j++)
        for(int i=0; lg2[n-i]>=j; i++)
            max_[i][j]=max(max_[i][j-1],max_[i+(1<<(j-1))][j-1]);
}
int RMQ_Max(int x,int y)
{
    int k=lg2[y-x+1];
    return max(max_[x][k],max_[y-(1<<k)+1][k]);
}
void ST2(int *a,int n)
{
    lg2[0]=-1;
    for(int i=1; i<=n; i++)
        lg2[i]=lg2[i-1]+(i&(i-1)?0:1);
    for(int i=0; i<n; i++) min_[i][0]=a[i]; //a第一个数从零开始
    for(int j=1; j<=lg2[n]; j++)
        for(int i=0; lg2[n-i]>=j; i++)
            min_[i][j]=min(min_[i][j-1],min_[i+(1<<(j-1))][j-1]);
}
int RMQ_Min(int x,int y)
{
    int k=lg2[y-x+1];
    return min(min_[x][k],min_[y-(1<<k)+1][k]);
}

struct San :vector<int>
{
    void prepare()
    {
        sort(begin(),end());
        erase(unique(begin(),end()),end());
    }
    int get(int x)
    {
        return (int)(lower_bound(begin(),end(),x)-begin())+1;
    }
} rank;
int n,a[N],pos[N];
bool cal()
{
    ST2(a , n);
    ST(pos , n);

    int ok = 0;
    for(int i = 1 ; i<=n ; i++)
    {
        for(int j = i + 1 ; j<=n ; j++) if(a[j] > a[i])
            {
                int L = a[i] + 1;
                int R = a[j] - 1;
                if(L > R) continue;
                int max_pos = RMQ_Max(L , R);
                R = (--max_pos);
                L = j + 1;
                if(L > R) continue;
                int min_val = RMQ_Min(L , R);
                if(min_val < a[i]){
                    ok = 1;
                    break;
                }
            }
        if(ok) break;
    }
    return ok;
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d",&n);
        rank.clear();
        rep1(i , 1 , n) scanf("%d",&a[i]),rank.push_back(a[i]);
        rank.prepare();
        for(int i = 1 ; i<=n ; i++)
        {
            a[i] = rank.get(a[i]);
            pos[a[i]] = i;
        }
        if(cal())
        {
            printf("YES\n");
            continue;
        }
        for(int i = 1 ; i<=n ; i++)
        {
            a[i] = n + 1 - a[i];
            pos[a[i]] = i;
        }
        if(cal())
        {
            printf("YES\n");
            continue;
        }
        printf("%s\n","NO");
    }
    return 0;
}