首页 > 其他分享> 文章详细

CQOI2014 数三角形

2019-02-13 08:53:08 阅读：257 来源： 互联网

标签：ch int long while RG CQOI2014 三角形 define

数三角形

直接正着做不好做，那么考虑整体减去部分。

首先算出如果可以3点共线，总共多少个，显然是C[n*m][3]。

接下来考虑共线的情况：

横排的显然是n*C[m][3]，竖排的就是m*C[n][3]。

比较麻烦一点的是斜着的，乍一看好像统计斜着的n³？

实际上，有这样一个结论：

对于点(a,b) (x,y)连成的线段而言(其中a>x,b>y)，

在它们中间有gcd(a-x,b-y)-1个整点。

所以我们不妨把一个点恒定为（1,1）点。（整体+1之后）

然后枚举每个点，计算贡献就好了。

#include<bits/stdc++.h>
#define RG register
#define IL inline
#define LL long long 
using namespace std;

IL int gi() {
    RG int x=0,w=0; char ch=0;
    while (ch<'0'||ch>'9') {if (ch=='-') w=1;ch=getchar();}
    while (ch>='0'&&ch<='9') x=(x<<3)+(x<<1)+(ch^48),ch=getchar();
    return w?-x:x;
}

const int N=1e3+10;
const int MAX=1002005;

int n,m;
LL rep,C[MAX+10][4];

int main ()
{
    //freopen ("count.in","r",stdin);
    //freopen ("count.out","w",stdout);
    RG int i,j;
    n=gi()+1,m=gi()+1;
    for (i=0;i<=MAX;++i) C[i][0]=1; 
    for (i=1;i<=MAX;++i)
        for (j=1;j<=3&&j<=i;++j)
            C[i][j]=C[i-1][j]+C[i-1][j-1];
    for (i=2;i<=n;++i)
        for (j=2;j<=m;++j) 
            rep+=(n-i+1)*(m-j+1)*(__gcd(i-1,j-1)-1);
    printf("%lld\n",C[n*m][3]-(n*C[m][3]+m*C[n][3]+rep*2));
    return 0;
}

BY BHLLX

标签：ch,int,long,while,RG,CQOI2014,三角形,define
来源： https://www.cnblogs.com/Bhllx/p/10367991.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

ICode9

CQOI2014 数三角形