洛谷-3166 [CQOI2014]数三角形

最新推荐文章于 2025-07-28 12:46:22 发布

原创最新推荐文章于 2025-07-28 12:46:22 发布 · 596 阅读

1 ·

本内容遵循CC 4.0 BY-SA版权协议

收录于

算法题专栏收录该内容

310 篇文章

订阅专栏

本文介绍了一种算法，用于计算在给定的nxm网格上，三点都在格点上且不共线的三角形数量。通过计算所有点的组合并减去共线点的方法，实现了高效准确的计算。

题目描述
给定一个nxm的网格，请计算三点都在格点上的三角形共有多少个。下图为4x4的网格上的一个三角形。注意三角形的三点不能共线。
输入格式
输入一行，包含两个空格分隔的正整数m和n。
输出格式
输出一个正整数，为所求三角形数量。

输入输出样例
输入 #1
2 2

输出 #1
76

说明/提示
数据范围
bzoj上是1<=m,n<=1000

解释：首先算出所有点的组合，再减去直线就好了，即 $C (n * m, 3)$ $-$ 横直线 $-$ 竖直线 $-$ 斜直线,竖和横都好计算，斜的话，我们固定一点原点，在枚举另一个端点，在利用结论 $g c d (x, y) - 1$ 为线段中点个数，再平移线段统计即 $* (n - x) * (m - y)$ ，考虑的矩阵的对称性，对于斜率k>0和k<0只需要计算一个再乘2就好了。


#include<iostream>
using namespace std;
long long n=0,m=0;
long long ret=0;
long long C(long long x){
    if(x<3) return 0;
    return (x-2)*(x-1)*x/6;
}
long long gcd(long long a,long long b){
    if(b==0) return a;
    return gcd(b,a%b);
}
int main(){
    cin>>m>>n;
    n++;m++;
    ret+=C(n*m);
    ret-=n*C(m);
    ret-=m*C(n);
    for(int i=1;i<n;i++){
        for(int j=1;j<m;j++){
            ret-=(n-i)*(m-j)*(gcd(i,j)-1)*2;
        }
    }
    cout<<ret<<endl;
    return 0;
}