[打基础]luogu2181对角线——计数原理
啦啦啦我ysw又回来啦!之后大概会准备打acm,暑假尽量复习复习,因为已经快两年没碰oi了,最多也就高三noip前学弟学妹出题讲题,所以从这一篇blog开始大概会有一系列“打基础”的blog,既是复习也是重新学一遍叭~高中的时候学得就不太扎实
题意:对于一个任何三条对角线不会交于同一点的凸$n$边形,问对角线交点个数,$n\leq 10^5$。

例如$n=6$的情形下答案为15
似乎我们可以随意给一个端点编号为1,按照逆时针/顺时针顺序依次编号$1,2,3……n$,这样一来我们可以先只考虑以点1为其中一个端点出发的所有对角线(图中红线)对应的交点,这个问题似乎好处理了许多,因为如果另一条对角线要和它有交点(蓝线)当且仅当蓝线的两个端点在红线的两侧!(嗯!这不就是计算几何里面判断直线是否有交点的方法嘛!?)
于是我们想,对于编号为$k$的点,1~k内从点2到点$k-1$是可以作为其中一个起点的,一共是$k-1-2+1=k-2$个点,同样的另一侧从$k+1$到$n$可以作为对角线另一点,一共是$n-(k+1)+1=n-k$个点
乘法原理得到$(k-2)(n-k)$,接着只要对所有情形进行求和就行了,对于$k$来说最小取到3(2的话就不是对角线了),最大取到$n-1$
所以对于以任意一点为端点的所有对角线的交点便是:$m=\sum_{k=3}^{n-1}(n-k)(k-2)$,对$n$个点进行同样的操作就得到$nm$,而注意到如果这样写,每个交点对应两条对角线的四个端点都作为了一次起点,也就是说每个交点都被算了四次,所以我们要的答案应该是$ans=\frac{1}{4}nm$,式子推出来啦!
接着就可以编程求解了…首先整个问题的时间复杂度是$O(n)$,轻轻松松通过,不过我们注意到对应答案的上界是$O(n^3)$级别的也就是$10^{15}$级别的差不多,一开始我只开了个long long交上去然后挂了两个点…后面改成unsigned long long 就过了…emm好吧
话说回来我们再回头看看那个和式$m=\sum_{k=3}^{n-1}(n-k)(k-2)$,好像我们并不陌生呀!
里面展开就是$nk-2n-k^2-2k$我们对k求和就稍微整理一下变成$\sum_{k=3}^{n-1}[k(n-2)-2n-k^2]$,嗯第一项是等差数列,第二项是对于$k$的常数列,第三项是平方和,我们有结论$\sum_{k=1}^{n}k^2=\frac{1}{6}n(n+1)(2n+1)$,对于这题只要改个$n$再扣掉前两项就行
整个和式拆开再化简就得到这题的答案可以写成:
$\frac{n}{4}[\frac{1}{2}(n+2)^2-2(n^2-3n)-\frac{1}{6}n(n-1)(2n-1)+5]$
整个问题变成$O(1)$的复杂度啦!~(虽然好像没什么用,哈哈哈)
[打基础]luogu2181对角线——计数原理的更多相关文章
- STM32F4_TIM基本延时(计数原理)
Ⅰ.概述 STM32的TIM定时器分为三类:基本定时器.通用定时器和高级定时器.从分类来看就知道STM32的定时器功能是非常强大的,但是,功能强大了,软件配置定时器就相对复杂多了.很多初学者甚至工作了 ...
- Luogu 1351 NOIP 2014 联合权值(贪心,计数原理)
Luogu 1351 NOIP 2014 联合权值(贪心,计数原理) Description 无向连通图 G 有 n 个点,n-1 条边.点从 1 到 n 依次编号,编号为 i 的点的权值为 Wi, ...
- 【计数原理】【UVA11538】 Chess Queen
传送门 Description 给你一个n*m的棋盘,在棋盘上放置一黑一白两个皇后,求两个皇后能够互相攻击的方案个数 Input 多组数据,每组数据包括: 一行,为n和m 输入结束标志为n=m=0. ...
- Objective-C 引用计数原理
http://www.cocoachina.com/ios/20160112/14933.html 引用计数如何存储 有些对象如果支持使用 TaggedPointer,苹果会直接将其指针值作为引用计数 ...
- oc引用计数原理-引用计数相关变化
http://blog.csdn.net/null29/article/details/71191044 在 32 位环境下,对象的引用计数都保存在一个外部的表中,每一个对象的 Retain 操作,实 ...
- polay计数原理
公式: Burnside引理: 1/|G|*(C(π1)+C(π2)+C(π3)+.....+C(πn)): C(π):指不同置换下的等价类数.例如π=(123)(3)(45)(6)(7),X={1, ...
- 数学计数原理(Pólya):POJ 1286 Necklace of Beads
Necklace of Beads Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 7763 Accepted: 3247 ...
- 组合数学(Pólya计数原理):UvaOJ 10601 Cubes
Cubes You are given 12 rods of equal length. Each of them is colored in certain color. Your task is ...
- 数学计数原理(Pólya,高精度):SGU 294 He's Circles
He's Circles He wrote n letters "X" and "E" in a circle. He thought that there ...
随机推荐
- windbg 分析cpu异常
1. !threadpool 查看当前CPU状况 线程数等等 2. !runaway 查看那几个线程使用的高 建议多抓几个dump 然后确定到底是哪个线程 3. ~线程IDs 跳转到那个 ...
- 金九银十想面BAT?那这些JDK 动态代理的面试点你一定要知道
一.什么是代理 代理是一种常用的设计模式,其目的就是为其他对象提供一个代理以控制对某个对象的访问.代理类负责为委托类预处理消息,过滤消息并转发消息,以及进行消息被委托类执行后的后续处理. 代理模式UM ...
- leetcode117. 填充每个节点的下一个右侧节点指针 II
给定一个二叉树struct Node { int val; Node *left; Node *right; Node *next;}填充它的每个 next 指针,让这个指针指向其下一个右侧节 ...
- PHP 统计文件数和文件大小
/** * 统计文件数和文件大小 */private function getFileCacheCount($pathName){ $data = [ 'num' => 0, 'size' =& ...
- 一口气带你读懂80年IT发展史
计算机的发展历史有多长?真正意义上的计算机诞生,距今也只有80多年的时间.80年,对于每一个人来说,是很长的时间,但对于整个历史来说,只是短短的一瞬间.这八十多年只是整段历史中的一粒尘埃罢了,但却对这 ...
- MFC二进制文件读取
1.mfc Document-vew doc类中读取 doc类中读取,在Vew类中可直接使用. 在菜单栏NewFile/OpenFile 后,系统自动调用Serialize()函数 if :写入文件 ...
- node-sass版本问题
node-sass sass-loader的问题 出现了版本的问题 版本太高 版本不兼容解决方法: cnpm i node-sass@4.14.1 cnpm i sass-loader@7.3.1 - ...
- PyQt学习随笔:Model/View开发时从Model相关类派生自定义类需要注意的问题
在<PyQt学习随笔:重写setData方法截获Model/View中视图数据项编辑的注意事项>介绍的方法,从Model相关类派生自定义类,通过重写setData方法以获取View中数据的 ...
- linux进程管理(linux命令安装、进程生命周期、进程状态)
1 linux下如何杀掉进程 1)找到包名所占用的端口: ps aux | grep cbs_portal-1.0.1.jar(包名) 2)杀掉进程: kill 10942(端口号) PS: //-- ...
- 深入剖析HashMap
前言 很高兴遇见你~ HashMap是一个非常重要的集合,日常使用也非常的频繁,同时也是面试重点.本文并不打算讲解基础的使用api,而是深入HashMap的底层,讲解关于HashMap的重点知识.需要 ...