BZOJ 3505 数三角形

原创 2017年08月13日 06:16:21

Description

给定一个nxm的网格,请计算三点都在格点上的三角形共有多少个。下图为4x4的网格上的一个三角形。

注意三角形的三点不能共线。

Input

输入一行,包含两个空格分隔的正整数m和n。

Output


输出一个正整数,为所求三角形数量。

Sample Input


2 2

Sample Output

76


数据范围
1<=m,n<=1000

如果不考虑三点共线的限制,那么方案数就是C(n*m,3) 然后来考虑不符合条件的数目


如果是在同一行或同一列上好说,就是C(n,3)*m+C(m,3)*n


那么在对角线上的呢? .. 用gcd搞一下,大致可理解为枚举每种矩形,然后再计算这种矩形


内的对角线在整个棋盘中有多少种,再用Ans减去就好了

#include <iostream>
#include <cstring>
#include <algorithm>
#include <stdio.h>
typedef long long ll;
#define MAXN 1000005
ll n,m,tot,C[MAXN][4],Ans;
 
 
template<typename _t>
inline _t read(){
    _t x=0,f=1;
    char ch=getchar();
    for(;ch>'9'||ch<'0';ch=getchar())if(ch=='-')f=-f;
    for(;ch>='0'&&ch<='9';ch=getchar())x=x*10+(ch^48);
    return x*f;
}
inline ll gcd(ll x,ll y){return y==0?x:gcd(y,x%y);}
 
void Get_C(){
    tot = n * m;
    C[0][0]=1;
    for(int i=1;i<=tot;i++){
        C[i][0]=1;
        for(int j=1;j<=3;j++)
            C[i][j]=C[i-1][j-1]+C[i-1][j];
    }
}
 
int main(){
    n=read<int>();m=read<int>();
    ++n;++m;Get_C();
    Ans += C[tot][3];Ans -= C[n][3] * m + C[m][3] * n;
    for(int i=1;i<n;i++)
        for(int j=1;j<m;j++){
            register ll kk = gcd(i,j);
            if(kk > 1)Ans -= 2*(kk-1)*(n-i)*(m-j);
        }
    printf("%lld\n",Ans);
}







版权声明:传送门:blog.csdn.net/zyf_2002 举报

相关文章推荐

bzoj 3505: [Cqoi2014]数三角形

3505: [Cqoi2014]数三角形 Time Limit: 10 Sec  Memory Limit: 128 MB Submit: 1547  Solved: 949 [Submit][Sta...

【BZOJ 3505】 [Cqoi2014]数三角形

思路题~

我是如何成为一名python大咖的?

人生苦短,都说必须python,那么我分享下我是如何从小白成为Python资深开发者的吧。2014年我大学刚毕业..

BZOJ 3505 数三角形 (数论 组合数 gcd)

BZOJ 3505 数三角形 (数论 组合数)Description 给定一个nxm的网格,请计算三点都在格点上的三角形共有多少个。下图为4x4的网格上的一个三角形。注意三角形的三点不能共线。Inp...
  • w4149
  • w4149
  • 2017-07-05 16:51
  • 40

bzoj 3505 数三角形 组合数 解题报告

数三角形

BZOJ 3251: 树上三角形

Description 给定一大小为n的有点权树,每次询问一对点(u,v),问是否能在u到v的简单路径上取三个点权,以这三个权值为边长构成一个三角形。同时还支持单点修改。   Inp...

BZOJ 3251 树上三角形 数学

BZOJ 3251 树上三角形 数学

BZOJ 1914: [Usaco2010 OPen]Triangle Counting 数三角形

Description 在一只大灰狼偷偷潜入Farmer Don的牛群被群牛发现后,贝西现在不得不履行着她站岗的职责。从她的守卫塔向下瞭望简直就是一件烦透了的事情。她决定做一些开发智力的小练习,...

三角形数

三角形数 时间限制(普通/Java):1000MS/3000MS          运行内存限制:65536KByte 总提交:460            测试通过:114 描述...

三角形套三角形

/*   * Copyright (c) 2011, 烟台大学计算机学院   * All rights reserved.   * 文件名称:test.cpp   * 作 ...

HDU 2393 Higher Math

Higher Math Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) To...
返回顶部
收藏助手
不良信息举报
您举报文章:深度学习:神经网络中的前向传播和反向传播算法推导
举报原因:
原因补充:

(最多只允许输入30个字)