BZOJ3505: [Cqoi2014]数三角形（洛谷P3166）

最新推荐文章于 2019-10-23 19:38:30 发布

forezxl

最新推荐文章于 2019-10-23 19:38:30 发布

阅读量292

点赞数 1

分类专栏：洛谷 BZOJ 组合---组合数学相关蒟蒻zxl的Blog专栏文章标签： BZOJ 洛谷组合

本文链接：https://blog.csdn.net/a1799342217/article/details/82972878

版权

蒟蒻zxl的Blog专栏同时被 3 个专栏收录

460 篇文章 6 订阅

订阅专栏

洛谷

292 篇文章 1 订阅

订阅专栏

BZOJ

281 篇文章 1 订阅

订阅专栏

组合

BZOJ题目传送门
 洛谷题目传送门

直接计算比较烦，考虑全部方案减不合法的方案。不合法的方案就是三个点都在一条直线上。除去横的和竖的线，其它的枚一枚求一求gcd乘一乘就好了。

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 1005
using namespace std;
typedef long long LL;
LL n,m,ans;
#define C(x) (x*(x-1)*(x-2)/6)
LL gcd(LL a,LL b){ return !b?a:gcd(b,a%b); }
int main(){
	scanf("%lld%lld",&n,&m),n++,m++;
	ans=C(n*m)-m*C(n)-n*C(m);
	for (LL i=1;i<n;i++)
		for (LL j=1,r;j<m;j++)
			if ((r=gcd(i,j)+1)>2) ans-=(r-2)*(n-i)*(m-j)*2;
	return printf("%lld\n",ans),0;
}