给定一个n*m的网格，求面积为奇数的正方形有多少个.
首先是n*m个面积为1的，然后剩下的要么是边长为奇数，要么被这样一个奇数边长所包围。
原因如下：
对于一个边长不平行于坐标抽的正方形，其边长一定是某个长方形的对角线，而且长方形长宽a,b一定是一奇数，一偶数，这样area = a^2+b^2才是奇数。
所以可以对任何奇数i <= min(n, m) 求出这样的边长正方形以及被包围的正方形个数。
注意对于一个奇数例如5，被包围的正方形可以是以1和4的对角线，2和3的对角线为边，这样对任何一个奇数i，被包围的正方形有i/2个，根据对称性还应该*2。
w-i+1表示宽方向能够移动的次数，l-i+1表示长度方向能够移动的次数，例如图5，长和宽方向均能移动2次。
ans = n*m（单位正方形） +
{2*（w-i+1)*(l-i+1)*(i/2) + (w-i+1)*(l-i+1)} i为 3<=i<=min(n, m)的奇数。
后来在最后完成这篇博客的时候想到能不能把式子简化一下呢，然后化简就得到这个：ans = ∑{ (w*l+w+l+1)*i – (w+l+2)*i*i + i*i*i} i为1到min（l,w)之间的奇数。
甚至可以推广到面积为偶数的情况：
公式：ans = ∑{ (w*l+w+l+1)*i – (w+l+2)*i*i + i*i*i} i为1到min（l,w)之间的偶数。
如果是求所有的正方形个数的话，只要把i从1取到min（w,l)就行了。
注意程序中预处理一下i，i*i，i*i*i的和。
下面只上化简后公式实现的代码好了。
time：29ms
//Date: 20140211
#include <iostream>
#include <sstream>
#include <cstdio>
#include <climits>
#include <ctime>
#include <cstring>
#include <cstdlib>
#include <string>
#include <stack>
#include <map>
#include <cmath>
#include <vector>
#include <queue>
#include <algorithm>
#define esp 1e-3
#define pi acos(-1.0)
#define inf 0x0f0f0f0f
#define pb push_back
#define lson l, m, rt<<1
#define rson m+1, r, rt<<1|1
#define mp(a, b) make_pair((a), (b))
#define in freopen("test_in.txt", "r", stdin);
#define out freopen("test_out.txt", "w", stdout);
#define bug puts("********))))))");
#define inout in out
#define stop system("pause");
#define PRD(a) printf("%d\n",(a))
#define PRLD(a) printf("%lld\n", (a))
#define PRID(a) printf("%I64d\n", (a))
#define PRU(a) printf("%u\n", (a))
#define PRLU(a) printf("%llu\n", (a))
#define PRIU(a) printf("%I64u\n", (a))
#define SET(a, v) memset(a, (v), sizeof(a))
#define READ(a, n) {REP(i, n) cin>>a[i];}
#define REP(i, n) for(int i = 0; i < (n); i++)
#define Rep(i, base, n) for(int i = base; i < n; i++)
#define REPS(s) for(int i = 0; s[i]; i++)
#define pf(x) ((x)*(x))
#define Log(a, b) (log((double)b)/log((double)a))
#define Srand() srand((int)time(0))
#define random(number) (rand()%number)
#define random_range(a, b) (int)(((double)rand()/RAND_MAX)*(b-a) + a) /*
1. 点分治，atan和atan2返回弧度值，atan值域为（-pi/2, pi/2), atan2值域为(-pi, pi) ；
2. exp(x)用于计算e^x ；
3. log2和log10分别用来计算对数, log默认以e为底 ；
4. sort中比较函数的原型:bool cmp(const Type &a, const Type &b);
5. lower_bound和upper_bound返回的是相应下标对应的指针 ；
6. dfs时注意利用强剪枝和避免重复状态进行优化 ；
7. 尽量减少不必要的状态表示的维度 ；
8. greater<T> () less<T> () ；
9. 尽量少用strlen,尤其是在递归深度较大，字符串较长的时候，容易超时；少用memset ；
10.不要在函数里面开数组，易暴栈 ；
*/ using namespace std;
typedef long long LL;
typedef unsigned long long ULL;
typedef vector<int> VI;
typedef pair<int,int> pii;
typedef vector<pii> VII;
typedef vector<pii, int> VIII;
typedef VI:: iterator IT;
typedef map<string, int> Mps;
typedef map<int, int> Mpi;
typedef map<int, pii> Mpii;
typedef map<pii, int> Mpiii; template<class T> inline const T& Max(const T& a, const T& b)
{
return a < b ? b : a;
}
template<class T> inline const T& Min(const T& a, const T& b)
{
return a < b ? a : b;
}
template<class T> inline void checkMax(T& a, const T& b)
{
if(a < b) a = b;
}
template<class T> inline void checkMin(T& a, const T& b)
{
if(a > b) a = b;
} const int maxn = + ;
LL f[maxn], g[maxn], h[maxn];
void pre()
{
f[] = g[] = h[] = ;
for(int i = ; i < maxn; i += ) {
f[i] = f[i-] + i;
g[i] = g[i-] + (LL)i*i;
h[i] = h[i-] + (LL)i*i*i;
}
}
int main()
{ LL n, m;
pre();
while(scanf("%lld%lld", &n, &m), n||m) {
LL k = min(n, m);
if((k&) == )
k--;
LL ans = (n*m+n+m+)*f[k]-(n+m+)*g[k]+h[k];
printf("%lld\n", ans);
}
return ;
}