Uva 10177 (2/3/4)-D Sqr/Rects/Cubes/Boxes?

给出一个N*N的正方形方格，问能找到多少个正方形？矩形呢（不包括正方形）？

题目扩展到了3维和4维，三维的如下图，对应的是正方体和立方体（不包括正方体）

用ans1[k]表示k维下的正方形个数，ans2[k]表示k维下矩形个数，推下公式：

先考虑二维的情况，然后推广。

边长为 i 的二维正方形个数：[a_2[i]=(N-i+1)^2]

总的正方形个数：
[ecause sum_{i=1}^{N}(N-i+1)^2=sum_{i=1}^{N}i^2]
[ herefore sum_{i=1}^{N}a[i]=sum_{i=1}^{N}i^2=frac{N(N+1)(2N+1)}{6}]
长i宽j的矩形（包括正方形）个数：

[c_2[i,j]=(N-i+1)^2(N-j+1)^2]

总的矩形个数（包括正方形）：
[sum_{i=1}^{N}sum_{j=1}^{N}c[i,j]=sum_{i=1}^{N}sum_{j=1}^{N}(i imes j)=[frac{N(N+1)}{2}]^2]

推广到k维。

边长为i的k维正方形个数：

[a_k[i]=(N-i+1)^k]

总的k维正方形个数：
[sum_{i+1}^{N}a_k[i]=sum_{i+1}^{N}i^k]

k维矩形（第j个边长分别为(i_j)包括k维正方形）个数：
[c_k[i]=prod_{j=1}^{k}(N-i_j+1),i_j=1,2,cdots,N]

k维矩形（包括k维正方形）总数：
[sum_{j=1}^{N}sum_{i_j=1}^{N}c_k[i]=sum_{j=1}^{N}sum_{i_j=1}^{N}prod_{j=1}^{k}i_j]
[ecause sum_{i=1}^{N}sum_{j=1}^{M}(i imes j)=1 imessum_{j=1}^{M}+2 imessum_{j=1}^{M}+cdots+N imessum_{j=1}^{M}=(sum_{i=1}^{N}i) imes(sum_{j=1}^{M}j)]
[ herefore sum_{j=1}^{N}sum_{i_j=1}^{N}prod_{j=1}^{k}i_j=prod_{j=1}^{k}(sum_{i_j=1}^{N}i_j)=prod_{j=1}^{k}frac{N(N+1)}{2}=[frac{N(N+1)}{2}]^k]
[ herefore sum_{j=1}^{N}sum_{i_j=1}^{N}c_k[i]=[frac{N(N+1)}{2}]^k]

结果：
[ ext{ans1[k]}=sum_{i=1}^{N}i^k]
[ ext{ans2[k]}=[frac{N(N+1)}{2}]^k- ext{ans1[k]}]

复杂度O(n)。

1A

# include <stdio.h>

int main()
{
	int n;
	while (scanf("%d", &n) != EOF) {
		long long ans1[3], ans2[3];
		ans1[0] = ans1[1] = ans1[2] = 0;
		for (int i = 1; i <= n; ++i) {
			ans1[0] += i*i;
			ans1[1] += i*i*i;
			ans1[2] += i*i*i*i;
		}
		long long int t = n*(n+1)/2;
		ans2[0] = t*t-ans1[0];
		ans2[1] = t*t*t-ans1[1];
		ans2[2] = t*t*t*t-ans1[2];
		printf("%lld %lld %lld %lld %lld %lld
", ans1[0], ans2[0], ans1[1], ans2[1], ans1[2], ans2[2]);
	}
	
	return 0;
}