多校&&HDU 5726 GCD
2016-07-19 20:42
246 查看
【题意】给了你n个数,然后q个查询,每次查询有多少个l,r区间队等于当前这个区间【L,R】的gcd!
【解题方法】区间gcd用线段树维护,RMQ也可!然后map来暴力统计就行了,复杂度nlognlogn!
【AC 代码】
#include <bits/stdc++.h>
using namespace std;
#define ll long long
const int maxn=1000010;
struct node{
int l,r;
ll gcd;
}Tree[maxn<<2];
ll a[maxn];
ll gcd(ll a,ll b)
{
return b==0?a:gcd(b,a%b);
}
void pushup(int rt)
{
Tree[rt].gcd=gcd(Tree[rt*2].gcd,Tree[rt*2+1].gcd);
}
void Build(int l,int r,int rt)
{
Tree[rt].l=l,Tree[rt].r=r;
if(l==r)
{
Tree[rt].gcd=a[l];
return ;
}
int mid=(l+r)/2;
Build(l,mid,rt*2);
Build(mid+1,r,rt*2+1);
pushup(rt);
}
ll queryans(int L,int R,int rt)
{
if(L<=Tree[rt].l&&Tree[rt].r<=R)
{
return Tree[rt].gcd;
}
int mid=(Tree[rt].l+Tree[rt].r)/2;
ll ans=0;
if(L<=mid) ans=gcd(ans,queryans(L,R,rt*2));
if(mid<R) ans=gcd(ans,queryans(L,R,rt*2+1));
return ans;
}
int n;
map<ll,ll>ans;
map<ll,ll>mp1;//mp1代表以x[i]结尾的所有区间的gcd的个数
map<ll,ll>mp2;//临时变量
int main()
{
int T;
scanf("%d",&T);
int cas=1;
while(T--)
{
scanf("%d",&n);
ans.clear();
mp1.clear();
mp2.clear();
for(int i=1; i<=n; i++) scanf("%I64d",&a[i]);
Build(1,n,1);
//puts("success");
mp1[a[1]]++;
ans[a[1]]++;
for(int i=2; i<=n; i++)
{
ll now=a[i];
mp2[now]++;
ans[now]++;
for(auto it=mp1.begin(); it!=mp1.end(); it++)
{
int nex=gcd(now,it->first);
ans[nex]+=it->second;
mp2[nex]+=it->second;
}
mp1.clear();
for(auto it=mp2.begin(); it!=mp2.end(); it++)
{
mp1[it->first]=it->second;
}
mp2.clear();
}
int q;
printf("Case #%d:\n",cas++);
scanf("%d",&q);
while(q--)
{
int l,r;
scanf("%d%d",&l,&r);
ll temp=queryans(l,r,1);
printf("%I64d %I64d\n",temp,ans[temp]);
}
}
}
【解题方法】区间gcd用线段树维护,RMQ也可!然后map来暴力统计就行了,复杂度nlognlogn!
【AC 代码】
#include <bits/stdc++.h>
using namespace std;
#define ll long long
const int maxn=1000010;
struct node{
int l,r;
ll gcd;
}Tree[maxn<<2];
ll a[maxn];
ll gcd(ll a,ll b)
{
return b==0?a:gcd(b,a%b);
}
void pushup(int rt)
{
Tree[rt].gcd=gcd(Tree[rt*2].gcd,Tree[rt*2+1].gcd);
}
void Build(int l,int r,int rt)
{
Tree[rt].l=l,Tree[rt].r=r;
if(l==r)
{
Tree[rt].gcd=a[l];
return ;
}
int mid=(l+r)/2;
Build(l,mid,rt*2);
Build(mid+1,r,rt*2+1);
pushup(rt);
}
ll queryans(int L,int R,int rt)
{
if(L<=Tree[rt].l&&Tree[rt].r<=R)
{
return Tree[rt].gcd;
}
int mid=(Tree[rt].l+Tree[rt].r)/2;
ll ans=0;
if(L<=mid) ans=gcd(ans,queryans(L,R,rt*2));
if(mid<R) ans=gcd(ans,queryans(L,R,rt*2+1));
return ans;
}
int n;
map<ll,ll>ans;
map<ll,ll>mp1;//mp1代表以x[i]结尾的所有区间的gcd的个数
map<ll,ll>mp2;//临时变量
int main()
{
int T;
scanf("%d",&T);
int cas=1;
while(T--)
{
scanf("%d",&n);
ans.clear();
mp1.clear();
mp2.clear();
for(int i=1; i<=n; i++) scanf("%I64d",&a[i]);
Build(1,n,1);
//puts("success");
mp1[a[1]]++;
ans[a[1]]++;
for(int i=2; i<=n; i++)
{
ll now=a[i];
mp2[now]++;
ans[now]++;
for(auto it=mp1.begin(); it!=mp1.end(); it++)
{
int nex=gcd(now,it->first);
ans[nex]+=it->second;
mp2[nex]+=it->second;
}
mp1.clear();
for(auto it=mp2.begin(); it!=mp2.end(); it++)
{
mp1[it->first]=it->second;
}
mp2.clear();
}
int q;
printf("Case #%d:\n",cas++);
scanf("%d",&q);
while(q--)
{
int l,r;
scanf("%d%d",&l,&r);
ll temp=queryans(l,r,1);
printf("%I64d %I64d\n",temp,ans[temp]);
}
}
}
相关文章推荐
- 基本规则
- 怎样在Git码云上上传项目
- 青蛙跳台阶问题
- 冒泡排序代码
- tuxedo linux 安装
- 动态规划--最长公共子串
- hdoj1872稳定排序(sort+结构体)
- 弄提纲【NOIP2016提高A组模拟7.19】
- mysql数据对象
- SpringMVC(二)ContextLoaderListenner
- LIS
- 美团Android自动化之旅—生成渠道包
- 学渣!面对疾风吧!哈撒给
- 【杭电oj】2085 - 核反应堆(数学)
- hdu5726 GCD(gcd +二分+rmq)
- log(三)——MDCAdapter之实现LogbackMDCAdapter
- 38. Count and Say
- KMP
- android.system.ErrnoException: open failed: ENOENT (No such file or directory) 07-19 20:27:45.011 66
- XZ_iOS之KVC