【HDU5857】Median(方法)

最新推荐文章于 2021-08-08 12:22:12 发布

oranges_c

最新推荐文章于 2021-08-08 12:22:12 发布

阅读量232

点赞数

CC 4.0 BY-SA版权

文章标签：方法

本文链接：https://blog.youkuaiyun.com/oranges_c/article/details/52262849

方法专栏收录该内容

4 篇文章

订阅专栏

本文介绍了一种用于查询两个区间合并后的中位数的算法实现。通过预先读取并存储数据，在给定两个区间的起始和结束位置后，能够高效地计算出合并区间内的中位数。此算法适用于数据规模较大且需要频繁查询的情况。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

记录一个菜逼的成长。。

按区间分类讨论。

#include <cstdio>
#include <iostream>
#include <cstring>
#include <string>
#include <algorithm>
#include <cstdlib>
#include <vector>
#include <set>
#include <map>
#include <queue>
#include <stack>
#include <list>
#include <deque>
#include <cctype>
#include <bitset>
#include <cmath>
using namespace std;
#define ALL(v) (v).begin(),(v).end()
#define cl(a) memset(a,0,sizeof(a))
#define bp __builtin_popcount
#define pb push_back
#define fin freopen("D://in.txt","r",stdin)
#define fout freopen("D://out.txt","w",stdout)
#define lson t<<1,l,mid
#define rson t<<1|1,mid+1,r
#define seglen (node[t].r-node[t].l+1)
typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;
typedef pair<LL,LL> PLL;
typedef vector<PII> VPII;
const int INF = 0x3f3f3f3f;
const int MOD = 1e9 + 7;
template <typename T>
inline void read(T &x){
    T ans=0;
    char last=' ',ch=getchar();
    while(ch<'0' || ch>'9')last=ch,ch=getchar();
    while(ch>='0' && ch<='9')ans=ans*10+ch-'0',ch=getchar();
    if(last=='-')ans=-ans;
    x = ans;
}
inline bool DBread(double &num)
{
    char in;double Dec=0.1;
    bool IsN=false,IsD=false;
    in=getchar();
    if(in==EOF) return false;
    while(in!='-'&&in!='.'&&(in<'0'||in>'9'))
        in=getchar();
    if(in=='-'){IsN=true;num=0;}
    else if(in=='.'){IsD=true;num=0;}
    else num=in-'0';
    if(!IsD){
        while(in=getchar(),in>='0'&&in<='9'){
            num*=10;num+=in-'0';}
    }
    if(in!='.'){
        if(IsN) num=-num;
            return true;
    }else{
        while(in=getchar(),in>='0'&&in<='9'){
                num+=Dec*(in-'0');Dec*=0.1;
        }
    }
    if(IsN) num=-num;
    return true;
}
template <typename T>
inline void write(T a) {
    if(a < 0) { putchar('-'); a = -a; }
    if(a >= 10) write(a / 10);
    putchar(a % 10 + '0');
}
/******************head***********************/
const int maxn = 100000 + 10;
LL a[maxn];
int n,m;
LL query(int l1,int r1,int l2,int r2,int pos)
{
    if(l1 + pos - 1 < l2)return a[l1+pos-1];
    if(pos > (l2 - l1) + 2 * (r1 - l2 + 1)){    //   1 2 3 4    3 4 5 6 7 8 9 10 11
        pos -= (l2 - l1) + 2 * (r1 - l2 + 1);
        return a[r1 + pos];
    }
    pos -= l2 - l1;
    int tmp = pos / 2;
    if(pos & 1)return a[l2+tmp];
    return a[l2+tmp-1];
}
int main()
{
    //fin;
    //fout;
    int T;
    double ans;read(T);
    while(T--){
        read(n),read(m);
        for( int i = 1; i <= n; i++ )
            read(a[i]);
        int l1,r1,l2,r2;
        while(m--){
            scanf("%d%d%d%d",&l1,&r1,&l2,&r2);
            if(l1 > l2)swap(l1,l2);
            if(r1 > r2)swap(r1,r2);
            int len1 = r1 - l1 + 1,len2 = r2 - l2 + 1;
            int sum = len1 + len2,mid = sum >> 1;
            if(r1 < l2){
                if(len1 == len2){
                    ans = (a[r1] + a[l2]) / 2.0;
                }
                else {
                    if(len1 > len2){
                        if(sum & 1) ans = a[l1 + mid];
                        else ans = (a[l1 + mid] + a[l1 + mid - 1]) / 2.0;
                    }
                    else {
                        if(sum & 1)ans = a[r2 - mid];
                        else ans = (a[r2 - mid] + a[r2 - mid + 1]) / 2.0;
                    }
                }
            }
            else {
                if(sum & 1)ans = query(l1,r1,l2,r2,sum/2+1);
                else ans = (query(l1,r1,l2,r2,sum/2) + query(l1,r1,l2,r2,sum/2+1)) / 2.0;
            }
            printf("%.1f\n",ans);
        }
    }
    return 0;
}