P1890 gcd区间
我一开始80分
暴力，模拟
100做法
dp
O(n^2+m)
f[i][j]表示i到j的 gcd
初始化
f[i][i]=i;
f[i][j]=gcd(f[i][j-1],a[j]);
这样查询就到了O(1)
80代码
#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#define p(a) putchar(a)
#define g() getchar()
#define inf 2147483647
#define For(i,a,b) for(register int i=a;i<=b;i++)
using namespace std; void in(int &x)
{
int y=;
char c=g();x=;
while(c<''||c>'')
{
if(c=='-')
y=-;
c=g();
}
while(c<=''&&c>='')x=x*+c-'',c=g();
x*=y;
} void o(int x)
{
if(x<)
{
p('-');
x=-x;
}
if(x>)o(x/);
p(x%+'');
} int gcd(int a,int b)
{
return (b==?a:gcd(b,a%b));
}
int n,m;
int a[];
int l,r,g;
int main()
{
in(n),in(m);
For(i,,n)
in(a[i]);
For(i,,m)
{
in(l),in(r);
g=a[l];
For(j,l+,r)
g=gcd(g,a[j]);
o(g),p('\n');
}
return ;
}
100分
#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#define p(a) putchar(a)
#define g() getchar()
#define inf 2147483647
#define For(i,a,b) for(register int i=a;i<=b;i++)
using namespace std; void in(int &x)
{
int y=;
char c=g();x=;
while(c<''||c>'')
{
if(c=='-')
y=-;
c=g();
}
while(c<=''&&c>='')x=x*+c-'',c=g();
x*=y;
} void o(int x)
{
if(x<)
{
p('-');
x=-x;
}
if(x>)o(x/);
p(x%+'');
} int gcd(int a,int b)
{
return (b==?a:gcd(b,a%b));
}
int n,m;
int a[];
int l,r,g; int f[][]; int main()
{
in(n),in(m);
For(i,,n)
in(a[i]);
For(i,,n)
f[i][i]=a[i]; For(i,,n)
For(j,i+,n)
f[i][j]=gcd(f[i][j-],a[j]);
For(i,,m)
{
in(l),in(r);
o(f[l][r]),p('\n');
}
return ;
}