高速公路 (Highway,CERC 2006,LA 3720)

感谢 http://www.cnblogs.com/zhexipinnong/archive/2013/05/07/3064893.html

#include <iostream>
#include <string.h>
#include <string>
#include <fstream>
#include <algorithm>
#include <stdio.h>
#include <vector>
#include <queue>
#include <set>
#include <cmath>
using namespace std;
const double eps = 1e-8;
const int INF=0x7fffffff;
unsigned long long uINF = ~0LL;
#define MAXN 10000007
typedef long long LL;
LL gcd(LL a,LL b)
{
    return b==0?a:gcd(b,a%b);
}
int g[301][301];
void init()
{
    for(int i=1;i<=300;i++)
        for(int j=1;j<=300;j++)
        g[i][j]=gcd(i,j);
}
int main()
{
    int n,m;
    init();
    while(scanf("%d%d",&n,&m),n+m!=0)
    {
        LL ans=0;
        n--;m--;
        if(n>m)swap(n,m);

        for(int i=1;i<=n;i++)
            for(int j=i;j<=m;j++)
            if(g[i][j]<=2)
            {
                LL temp=(n-i+1)*(m-j+1);
                if(i==j)
                {
                    if(g[i][j]==1)ans+=temp;
                    else ans-=temp;
                }
                else
                {
                    if(i<=m&&j<=n)temp+=(m-i+1)*(n-j+1);
                    if(g[i][j]==1)ans+=temp;
                    else ans-=temp;
                }
            }
            ans*=2;
        printf("%lld
",ans);
    }

    return 0;
}