P1890 gcd区间

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 ;
}