bzoj3505: [Cqoi2014]数三角形
2014-09-17 17:31
309 查看
首先共有(n+1)*(m+1)个点,所以先++n,++m
不考虑三点共线的情况,有(n*m)*(n*m-1)*(n*m-2)/6个
三点都在水平和竖直的有m*(m-1)*(m-2)*n/6+n*(n-1)*(n-2)*m/6个
然后考虑斜着的
枚举每个点作为第一个点,使后两个点坐标大于第一个
过第一个点做个水平和竖直的线,右上和左下的矩形可用
可用dp求出一个矩阵中第一个点在角上,三点共线的个数
不考虑三点共线的情况,有(n*m)*(n*m-1)*(n*m-2)/6个
三点都在水平和竖直的有m*(m-1)*(m-2)*n/6+n*(n-1)*(n-2)*m/6个
然后考虑斜着的
枚举每个点作为第一个点,使后两个点坐标大于第一个
过第一个点做个水平和竖直的线,右上和左下的矩形可用
可用dp求出一个矩阵中第一个点在角上,三点共线的个数
#include <cstdio> #include <algorithm> using namespace std; #define MAXN 1003 #define MAXP 170 int gcd(int x, int y) { while (y) { int t = x % y; x = y; y = t; } return x; } int n, m; int a[MAXN][MAXN], ct[MAXN][MAXN]; long long f[MAXN][MAXN]; int main() { scanf("%d%d", &n, &m); ++n, ++m; for (int i = 1; i <= n; ++i) for (int j = 1; j <= m; ++j) { int g = gcd(i, j); int x = i / g, y = j / g; a[i][j] = ct[x][y]++; } for (int i = 1; i <= n; ++i) f[i][1] = f[i - 1][1] + a[i][1]; for (int i = 1; i <= m; ++i) f[1][i] = f[1][i - 1] + a[1][i]; for (int i = 2; i <= n; ++i) for (int j = 2; j <= m; ++j) f[i][j] = f[i - 1][j] + f[i][j - 1] - f[i - 1][j - 1] + a[i][j]; int N = n * m; long long ans = (long long)N * (N - 1) * (N - 2) / 6; ans -= (long long)m * (m - 1) * (m - 2) / 6 * n + (long long)n * (n - 1) * (n - 2) / 6 * m; for (int i = 1; i <= n; ++i) for (int j = 1; j <= m; ++j) ans -= f[n - i][j - 1] + f[n - i][m - j]; printf("%lld\n", ans); return 0; }
相关文章推荐
- [Bzoj3505][Cqoi2014]数三角形
- BZOJ 3505: [Cqoi2014]数三角形|组合数学
- 【bzoj3505】[Cqoi2014]数三角形 容斥原理
- BZOJ-3505-数三角形-CQOI2014
- bzoj3505[CQOI2014]数三角形
- BZOJ3505 [Cqoi2014]数三角形
- 【bzoj3505】 Cqoi2014—数三角形
- BZOJ-3505-数三角形-CQOI2014
- BZOJ 3505 [Cqoi2014]数三角形
- BZOJ 3505 CQOI 2014 数三角形 数学
- 【BZOJ】3505: [Cqoi2014]数三角形
- [bzoj3505][CQOI2014]数三角形_组合数学
- 【BZOJ 3505】[Cqoi2014]数三角形 组合数学
- BZOJ 3505 【CQOI2014】 数三角形
- bzoj3505 [Cqoi2014]数三角形
- BZOJ3505 [Cqoi2014]数三角形
- 【组合数学】BZOJ3505 [Cqoi2014]数三角形
- BZOJ 3505 [Cqoi2014]数三角形(组合数学)
- bzoj 3505: [Cqoi2014]数三角形 组合数学