多校联赛之二分加rmq

最新推荐文章于 2019-06-06 13:50:57 发布

原创最新推荐文章于 2019-06-06 13:50:57 发布 · 407 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#多校联赛之二分加rmq

多校联赛之二分加rmq 专栏收录该内容

1 篇文章

订阅专栏

本文回顾了作者参加多校联赛首场比赛的经历，虽然未能解题成功，但深入剖析了两道典型题目：最小生成树问题及RMQ问题。通过他人代码学习到树形DP与RMQ的有效解决方法，分享了实现细节与算法优化技巧。

多校联赛开始了，第一场没有做出一道题，感觉的话挺对不起自己的，你的不努力换来的就是这个样子，永远都是，你的思路总是停留在浅薄的平面上，第一题是最小生成树，但是看出这个有什么用呢，关键在于求那个期望，那个期望值才是关键，果然，关键的东西还是掌握在少数人的手里，用的是一个深搜，也是很巧妙，据说是树形dp,我并没有学习过，但是我大致看明白代码了，还有一道是一个rmq的题目。其实刚开始我没有想到rmq，还是后来在别人的提醒下想到了rmq，但是这并没有什么卵用，关键在于后面的计数，后面的计数直接暴力的话是超时的，然后你就要观察这个数列的性质，最后发现他是递减的，然后二分，时间复杂度也降了下来，最后就完美的解决了，可是这还是掌握在少数人的手里面，很是对不起自己啊。

题目链接：点击打开链接

代码：直接贴别人代码了，自己原来是年少轻狂，写了一个n^2的代码，然后最后连数组的大小都判断不出来，最后还是超时，并没有什么卵用。

#include<stdio.h>
#include<algorithm>
using namespace std;
#include<bits/stdc++.h>
#include<map>
map<int ,long long >q;
int  treemaxzuida[100010][25];
int q1[100000],q2[100000];
int zuida(int a,int b)
{
    if(a<b)
        swap(a,b);
    return b==0?a:zuida(b,a%b);
}
void  rmq_zuida(int num)
{
    for(int j=1; (1<<j)<=num; ++j)
        for(int i=1; i+(1<<j)-1<=num; i++)
            treemaxzuida[i][j]=zuida(treemaxzuida[i][j-1],treemaxzuida[i+(1<<(j-1))][j-1]);

}
//int get_max(int a,int b)
//{
//    int k=int(log(b-a+1)/log(2));
//    return  zuida(treemaxzuida[a][k],treemaxzuida[b-(1<<k)+1][k]);
//}
int get_max(int l,int r)
{
    int k=0;
    while((1<<(k+1))<=r-l+1) k++;
    return zuida(treemaxzuida[l][k],treemaxzuida[r-(1<<k)+1][k]);
}
int main()
{
    int a,a1;
    scanf("%d",&a1);
    int cas=1;
    while(a1--)
    {
        scanf("%d",&a);
        memset(treemaxzuida,0,sizeof(treemaxzuida));
        for(int i=1; i<=a; i++)
            scanf("%d",&treemaxzuida[i][0]);
        rmq_zuida(a);
        q.clear();
        int nex,l,r,mid,d2,d1,ans;
        for(int i=1; i<=a; i++)
        {
            nex=i;
            while(nex<=a)
            {
                d1=get_max(i,nex);
                l=nex;
                r=a;
                while(l<=r)
                {
                    mid=(l+r)>>1;
                    d2=get_max(i,mid);
                    if(d2>=d1) l=mid+1,ans=mid;
                    else r=mid-1;
                }

                q[d1]+=(ans-nex)+1;
                nex=r+1;
            }
        }
//        for(int i=1; i<=a; i++)
//            for(int j=i; j<=a; j++)
//                q[get_max(i,j)]++;
        int b;
        memset(q1,0,sizeof(q1));
        memset(q2,0,sizeof(q2));
        scanf("%d",&b);
        for(int i=0; i<b; i++)
        {
            scanf("%d%d",&q1[i],&q2[i]);

        }
        printf("Case #%d:\n",cas);
        for(int i=0; i<b; i++)
            printf("%d %I64d\n",get_max(q1[i],q2[i]),q[get_max(q1[i],q2[i])]);
        cas++;

    }
    return 0;
}

自己写的代码，各种小心，还有一处爆int，我还wrong了好几发，，，，，，咳