[BZOJ3505][Cqoi2014]数三角形（组合数学+gcd）

最新推荐文章于 2019-05-29 15:53:53 发布

Clove_unique

最新推荐文章于 2019-05-29 15:53:53 发布

阅读量567

点赞数

分类专栏：题解数学相关组合数学省选

本文链接：https://blog.csdn.net/Clove_unique/article/details/54576420

版权

题解同时被 3 个专栏收录

1097 篇文章 4 订阅

订阅专栏

省选

211 篇文章 0 订阅

订阅专栏

数学相关

78 篇文章 0 订阅

订阅专栏

题目描述

传送门

题解

三角形数=随便选出3个点-3点共线的
总数就是 $C_{n*m}^3$
首先减去在横竖直线上的3点共线的
考虑如何计算斜线上三点共线的
然后n*m枚举一条斜线，这条斜线经过的整点就是(i,j)-1，所以斜线上有(i,j)-1中选择
将这条直线左右平移n-i，上下平移m-j都是合法的，所以总的选择应该是((i,j)-1) * (n-i)*(m-j)
向左倾斜和向右倾斜有两种情况，再*2

代码

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
using namespace std;
#define LL long long
#define N 1000000

LL n,m,c[N+1][5],ans;

LL gcd(LL a,LL b)
{
    if (!b) return a;
    else return gcd(b,a%b);
}
int main()
{
    for (int i=0;i<=N;++i) c[i][0]=1LL;
    for (int i=1;i<=N;++i)
        for (int j=1;j<=3;++j)
            c[i][j]=c[i-1][j]+c[i-1][j-1];
    scanf("%lld%lld",&n,&m);++n,++m;
    ans=c[n*m][3];
    ans-=c[n][3]*m+c[m][3]*n;
    for (LL i=1;i<n;++i)
        for (LL j=1;j<m;++j)
            ans-=(gcd(i,j)-1)*(n-i)*(m-j)*2;
    printf("%lld\n",ans);
}