矩阵 | 线段树&区间gcd

矩阵

**

• 题目
• 分析
• 代码

- 题目

- 分析

首先 我们需要理解x2∑i=x1 gcd（Ciy1,…,Ciy2)的含义：
我们将其拆分为：
gcd（Cx1y1,…,Cx1y2）+…
.+gcd（Cx2y1,…,Cx2y2）;
由于上述式子中每一个gcd里都包含一个公因子x，我们可以将这个x提取出来，那么这个式子就会转化为：(x1+…+x2)*gcd(y1,…,y2)

经过转化问题已经非常明显了:
题目中给出的a【】数组代表各个x，b【】数组代表各个y
我们需要做的就是通过前缀和求出x1到x2的和，再利用线段树维护b数组的区间gcd，两者相乘即为最终答案。

-代码

#include<cstdio>//矩阵
#include<iostream>
#include<algorithm>
#include<string.h>
using namespace std;
const int maxn=1e5+7;
long long a[maxn];//极端情况下前缀和为1e11会爆int，记得开long long
int b[maxn];
int q,n;
struct Node//线段树结构体
{
    int l,r;
    int gcd;//维护区间gcd
}tr[maxn*4];//切记开4倍空间

这一题的线段树是真的基础，板子题看看就好

int gcd(int a,int b)
{
    if(b==0)
    {
        return a;
    }
    return gcd(b,a%b);
}

gcd模板，很短，敲两遍就记住了

void build(int d,int l,int r)
{
    tr[d].l=l;
    tr[d].r=r;
    if(l==r)
    {
        tr[d].gcd=b[l];//叶子gcd等于自身不用多说吧
        return ;
    }
    int mid=(l+r)/2,lc=d*2,rc=d*2+1;
    build(lc,l,mid);
    build(rc,mid+1,r);
    tr[d].gcd=gcd(tr[lc].gcd,tr[rc].gcd);//回溯更新区间gcd
}

标准的建树模板

int query(int d,int l,int r)
{
    if(tr[d].l==l&&tr[d].r==r)//找到区间直接返回该区间的gcd
    {
        return tr[d].gcd;
    }
    int mid=(tr[d].l+tr[d].r)/2,lc=d*2,rc=d*2+1;
    if(l>mid)
    {
        return query(rc,l,r);
    }
    else if(r<=mid)
    {
        return query(lc,l,r);
    }
    else
    {
        return gcd(query(lc,l,mid),query(rc,mid+1,r));
    }
}

区间查找板子（这一题都不用写修改来着）

int main()
{
    int x1,x2,y1,y2;
    scanf("%d",&n);
    a[0]=0;//求前缀和必要的前缀0
    for(int i=1;i<=n;i++)//预处理前缀和
    {
        scanf("%lld",&a[i]);
        a[i]+=a[i-1];
    }
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&b[i]);
    }
    build(1,1,n);//建树
    scanf("%d",&q);
    while(q--)
    {
        long long ans;
        scanf("%d %d %d %d",&x1,&y1,&x2,&y2);
        ans=query(1,y1,y2)*(a[x2]-a[x1-1]);//查找y1到y2区间的gcd与x1到x2区间和相乘
        printf("%lld\n",ans);//直接输出答案即可
    }
    return 0;
}

That’s All