矩阵
**
- 题目
- 分析
- 代码
- 题目
- 分析
首先 我们需要理解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