数三角形 (组合数学+分类讨论)

数三角形

[Link](信息学奥赛一本通（C++版）在线评测系统 (ssoier.cn))

题意

给定一个$n\times m$的网格，请计算三点都在格点上的三角形共有多少个。

题解

首先我们先将长度都加一变成格点数，发现直接去统计合法的方案不是很方便，我们可以容斥一下总方案-不合法方案，一共有$n\times m$个点，从里面选三个点构成三角形，总方案为$C_{n\times m}^3$。

不合法的方案分为四类:

1. 于$x$轴平行：$n\times C_m^3$ （一共$n$行选一行，从$m$列中选三列）

2. 于$y$轴平行：$m\times C_n^3$

3. 三点在某一条斜率为$k$(正数)*的线段上上，假设左下角和右上角分别为$A,B$，且$A,B$横坐标之差为$j$,纵坐标之差为$i$,这样的线段一共有$(n-i)\ast (m-j)$条，然后线段$A,B$之间的点$C$一共有$gcd(i,j)-1$种情况。

   为什么有$gcd(i,j)-1$种情况？

   首先考虑特列从点$(0,0)$到点$(n,m)$，所有经过的点都应该在$y=\frac{m}{n}x$上，因为要经过格点，所以$x,y$应均为整数，若要$y$为整数则$\frac{m}{n}x$必是整数，即$n\mid mx$, 相当于$x$要包含$n除去gcd(n,m)$以外的$n$的约数才可以，所以$x$只能是

   $\frac{n}{(n,m)},2\frac{n}{(n,m)},3\frac{n}{(n,m)},...,(n,m)\frac{n}{(n,m)}$,所以一共有$gcd(n,m)$个，我们再除去右上角的格点，即一共有$gcd(n,m)-1$个情况，等价到线段情况即可。

​ 因此这种不合法方案有$(n-i)\ast (m-j)\ast (gcd(i,j)-1)$

4. 三点在某一条斜率为$k(负数)$与情况三同理，方案数一样

Code

#include <iostream>
#include <algorithm>
#include <cstring>
#include <cstdio>
#include <set>
#include <queue>
#include <vector>
#include <map>
#include <bitset>
#include <unordered_map>
#include <cmath> 
#include <stack>
#include <iomanip>
#include <deque> 
#include <sstream>
#define x first
#define y second
#define debug(x) cout<<#x<<":"<<x<<endl;
using namespace std;
typedef long double ld;
typedef long long LL;
typedef pair<int, int> PII;
typedef pair<double, double> PDD;
typedef unsigned long long ULL;
const int N = 1010, M = 2 * N, INF = 0x3f3f3f3f, mod = 1e9 + 7;
const double eps = 1e-8, pi = acos(-1), inf = 1e20;
#define tpyeinput int
inline char nc() {static char buf[1000000],*p1=buf,*p2=buf;return p1==p2&&(p2=(p1=buf)+fread(buf,1,1000000,stdin),p1==p2)?EOF:*p1++;}
inline void read(tpyeinput &sum) {char ch=nc();sum=0;while(!(ch>='0'&&ch<='9')) ch=nc();while(ch>='0'&&ch<='9') sum=(sum<<3)+(sum<<1)+(ch-48),ch=nc();}
int dx[] = {-1, 0, 1, 0}, dy[] = {0, 1, 0, -1};
int h[N], e[M], ne[M], w[M], idx;
void add(int a, int b, int v = 0) {
    e[idx] = b, w[idx] = v, ne[idx] = h[a], h[a] = idx ++;
}
int n, m, k;
int a[N];
int gcd(int a, int b) {
    return b ? gcd(b, a % b) : a;
}
LL C(int a, int b) {
    return (LL)a * (a - 1) * (a - 2) / 6;
}
int fact[N], infact[N];
int main() {
    ios::sync_with_stdio(false), cin.tie(0);
    cin >> n >> m;
    n ++, m ++;
    LL res = 0;
    res = (LL)C(n * m, 3)  - (LL)n * C(m, 3) - (LL)m * C(n, 3);
    for (int i = 1; i <= n; i ++)
        for (int j = 1; j <= m; j ++)
            res -= 2ll * (gcd(i, j) - 1) * (n - i) * (m - j);
    cout << res << endl;
    return 0;
}