Highways

http://acm.hust.edu.cn/vjudge/contest/view.action?cid=32290#problem/C

// File Name: uvala3720.cpp
// Author: bo_jwolf
// Created Time: 2013年09月28日 星期六 18:58:58

#include<vector>
#include<list>
#include<map>
#include<set>
#include<deque>
#include<stack>
#include<bitset>
#include<algorithm>
#include<functional>
#include<numeric>
#include<utility>
#include<sstream>
#include<iostream>
#include<iomanip>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<ctime>

using namespace std;
const int maxn = 305 ;
int gcd[ maxn ][ maxn ] ;

int Gcd( int a , int b )
{
	return b == 0 ? a : Gcd( b , a % b ) ;
}

int main()
{
	for( int i = 1 ; i < maxn ; ++i )
	{
		for( int j = i ; j < maxn ; ++j )
		{
			gcd[ i ][ j ] = Gcd( i , j ) ; 
		}
	}
	int n  , m , temp; 
	long long int ans ;
	while( scanf( "%d%d" , &n , &m ) != EOF )
	{
		if( !n && !m )
		{
			break ; 
		}
		ans = 0  ; 
		if( n > m )
			swap( n , m ) ;
		for( int i = 1 ; i < n ;++i )
		{
			for( int j = i ; j < m ; ++j )
			{
				if( gcd[ i ][ j ] <= 2 )
				{
					if( i == j )
					{
						temp = ( n - i  ) * ( m - j  ) ; 
						if( gcd[ i ][ j ] == 1 )
						{
							ans += temp ; 
						}
						else
						{
							ans -= temp ; 
						}
					}
					else
					{
						temp = ( n - i  ) * ( m - j  ) ;
						if( i < m && j < n )
						{
							temp += ( m - i  ) * ( n - j  ) ;
						}
						if( gcd[ i ][ j ] == 1 )
						{
							ans += temp ;
						}
						else
						{
							ans -= temp ;
						}
					}
				}
			}
		}
		ans *= 2 ;
		printf( "%lld\n" , ans ) ;
	}			
return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值