BZOJ3505&&洛谷P3166 [CQOI2014]数三角形

最新推荐文章于 2021-08-25 18:56:43 发布

AcerMo

最新推荐文章于 2021-08-25 18:56:43 发布

阅读量264

点赞数

分类专栏：数论&&组合数学

本文链接：https://blog.csdn.net/ACerAndAKer/article/details/81333931

版权

数论&&组合数学专栏收录该内容

49 篇文章 0 订阅

订阅专栏

发现除了三点共线，其他都行，那就减去三点共线的情况就好了

代码

//By AcerMo
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define lli long long int
using namespace std;
const lli M=1000500;
lli n,m;
lli f[M][4],ans;
inline lli gcd(lli a,lli b)
{
    if (!b) return a;
    return gcd(b,a%b);
}
signed main()
{
    cin>>n>>m;f[0][0]=1;n++;m++;
    for (lli i=1;i<=n*m;i++)
    {
        f[i][0]=1;
        for (lli k=1;k<=3;k++) f[i][k]=(f[i-1][k-1]+f[i-1][k]);
    }
    ans=f[n*m][3];
   	ans-=n*f[m][3];
   	ans-=m*f[n][3];
    for (lli i=1;i<n;i++)
    for (lli k=1;k<m;k++)
    {
        lli d=gcd(i,k);
        if (d>1) ans-=(d-1)*2*(n-i)*(m-k);	
    } 
    cout<<ans;
    return 0;
}