题解 P3166 [CQOI2014]数三角形

Posted suwakow

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了题解 P3166 [CQOI2014]数三角形相关的知识,希望对你有一定的参考价值。

蒟蒻的空间

做完之后看了看题解,怎么一篇和我思路一样的也没有...我好慌啊qwq(所以一定是窝太弱了看不懂dalao的思路)

好吧窝的方法确实很奇怪:

核心代码只有3行 输入 循环 输出 一气呵成 是题解中的豪杰

最重要的是

没有组合数 没有容斥 没有斜率 没有向量 DA☆ZE

(只有我们的好朋友gcd


咳咳 那么开始正题(敲黑板)

首先,我们定义一个网格被一个三角形完全覆盖,当且仅当这个三角形的三个顶点都在网格边界上,并且沿着网格内部任意一条线段把网格切开,一定会把三角形切成两部分。比如下面的例子就是一个完全覆盖(渣鼠绘):

技术图片

但是这个就不是:(因为沿红色竖线切开,并不能切到三角形)

技术图片

那么我们发现:

  • 每个顶点都在格点上的三角形,有且只有一个可以被它完全覆盖的网格。所以只要求出原矩形当中所有子矩形的完全覆盖三角形的数量,就可以不重不漏地找出顶点都在格点上的三角形。

——但是子矩形好多啊,枚举左上角和右下角,至少是\(n^2m^2\)的数量级哇qwq

——[恋符]MasterSpark.gif 注意到我们并不关心每个子矩形的位置,而只关心它们的长宽,以及长宽均相同的矩形的数目,所以枚举子矩形的长 \(i\) 和宽 \(j\),则\(i*j\)的矩形数量为\((n-i+1)*(m-j+1)\),子矩阵数量级降为\(nm\)

至此,原问题转化为 

**给定网格的长宽,迅速求解完全覆盖网格的三角形的数目**

继续观察,我们发现:

  • 如果某个三角形(暂且称之为\(ABC\))完全覆盖了某个网格(\(MNPQ\)),\(ABC\)一定有至少一个顶点在\(MNPQ\)的角上。

首先分析只有一个顶点在\(MNPQ\)角上的情况。不妨设\(A\)点与\(M\)点重合,为了使\(ABC\)完全覆盖\(MNPQ\)\(B\)\(C\)必须分别在\(NP\)\(PQ\)边上(如下图):

技术图片

显然对于一个\(i*j\)的网格(这里\(i\) \(j\)指的是空格的数量而非格点,上图\(i=6,j=10\)),固定顶点的位置有四种,每种对应的另外两个顶点的位置有\((i-1)*(j-1)\)种(\(B\) \(C\)不能与\(N\) \(P\) \(Q\)重合),共\(4*(i-1)*(j-1)\)种。

分析两个顶点在\(MNPQ\)角上的情况。不妨设\(A\)\(M\)重合。此时另外一个角上的点(不妨设为点\(B\))有三种情况:

1、\(B\)\(N\)重合。此时\(C\)一定在\(QP\)上。共\((i-1)\)种情况。

2、\(B\)\(Q\)重合。此时\(C\)一定在\(NP\)上。共\((j-1)\)种情况。

3、\(B\)\(P\)重合。

这是比较麻烦的一种状态,因为此时\(C\)点可以在网格中能构成三角形的任意一处。但是我们注意到,如果线段\(AB\)除了经过\(M\) \(P\)之外,还经过了一些其他格点,\(C\)是不能与它们重合的。

那么有多少个格点被\(AB\)穿过呢qwq?

显然,不包括\(AB\)本身,有\(gcd(i,j)-1\)个(至于为什么,请读者自己思考(明明就是你自己也不会证吧kora

所以第三种情况的方案数是\((i+1)*(j+1)-4-gcd(i,j)+1\)(这里\(-4\)是因为C点不能放在网格的四个角上)。

注意到以上三种情况都可以反转,从而得到另一组与其一一对应的方案。

分析三个顶点在\(MNPQ\)角上的情况。显然只有四种。

综上,对于一个长宽为\(i,j\)的网格,可以把它完全覆盖的三角形的个数

\(S=4*(i-1)*(j-1)+2*[(i-1)+(j-1)+(i+1)*(j+1)-4-gcd(i,j)+1]+4\)

\(=6*i*j-2*gcd(i,j)\)

枚举子矩阵的复杂度为\(mn\),单次求解\(gcd\)的复杂度为\(log(m+n)\),总复杂度\(O(mnlog(m+n))\),实际运行跑的飞起。

下面是AC代码:

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define R register
#define ll long long
using namespace std;

int getgcd(int a, int b)

    if (!b) return a;
    return getgcd(b, a % b);


int main()

    ll n, m, ans = 0;
    cin >> n >> m;
    for (R int i = 1; i <= n; ++i)
        for (R int j = 1; j <= m; ++j)
            ans += (n - i + 1) * (m - j + 1) * (6 * i * j - 2 * getgcd(i, j));
    cout << ans;
    return 0;

那么这篇文章就到这里,希望对您能有帮助。ありがとナスます~

以上是关于题解 P3166 [CQOI2014]数三角形的主要内容,如果未能解决你的问题,请参考以下文章

bzoj3505 / P3166 [CQOI2014]数三角形

洛谷P3166 [CQOI2014]数三角形

P3166 [CQOI2014]数三角形

Luogu P3166 [CQOI2014]数三角形 组合数学

BZOJ3505[Cqoi2014]数三角形 组合数

[计数][CQOI2014]数三角形