【百度之星2014~初赛(第二轮)解题报告】JZP Set
声明
笔者近期意外的发现 笔者的个人站点http://tiankonguse.com/ 的非常多文章被其他站点转载,可是转载时未声明文章来源或參考自 http://tiankonguse.com/ 站点,因此,笔者加入此条声明。
郑重声明:这篇记录《【百度之星2014~初赛(第二轮)解题报告】JZP Set》转载自 http://tiankonguse.com/ 的这条记录: id=668" target="_blank" style="word-wrap:break-word; word-break:normal; margin:0px; padding:0px; color:rgb(153,153,153); text-decoration:none">http://tiankonguse.com/record/record.php?id=668
前言
近期要毕业了,有半年没做比赛了.
这次參加百度之星第二轮娱乐一下.
如今写一下 JZP Set 这道题的的解题报告.
正文
题意
题意是:给你n个数(1到n),给你一个规则。问用这个规则能够得到多少个合法的集合.
详细规则是:一个合法集合里随意挑两个数,假设这两个数之和是偶数,这个偶数除以2得到的数也要在这个合法集合里.
比方: 3 和9 在集合里,3+9是偶数。所以 (3+9)/2 = 6 也要在这个集合里.然后 {3,6,9}就是一个合法的集合.
起初我理解错题意了,以为是随意的两个数字之和除以2.比方 (3 + 6)/2 = 4, 我以为4也要在集合里.
于是非常快得到一个 (n+1)*n/2 的公式.
后来学弟告诉我正确的题意.
分析
这个题的例子有10^5个,每一个例子最大是10^7.
错误题意
首先来看看我理解错误题意时,是怎么做的.
如果 F(n) 是 n 的答案的话, 则 F(n + 1) = F(n) + f(n+1).
f(n+1)里面肯定有 n+1, 我如果f(n+1) 这些集合的任一个集合 S 里的最小值是 a, 则 (a+n+1)/2 肯定在 S 里面。然后这样递归下去发现 a到n+1的全部数字都必须在 S 里面.
于是 f(n+1) 就是 n+ 1 了.
于是终于方程就是 F(n+1) = F(n) + n+ 1.
仅仅可惜这是错误的题意的做法.
正确题意
学弟告诉我正确的题意,还是能够非常快写出方程来
F(n+1) = F(n) + f(n + 1).
当中f(n+1) 是含有 n+ 1的合法的集合.
奇数偶数合法<
首先对于我上面找到的那些肯定是合法集合的一部分,仅仅是我遗漏了一些合法集合.
遗漏的肯定是非连续的了.
然后非常快能够想到 一个奇数和一个偶数构成的集合都是合法集合.
于是 f(n + 1) 就又加上含有一奇一偶的合法集合的数量了,仅仅是提交后WA了.
等差为奇数的组合
然后学弟随手写了一个集合 {3, 6, 9 } 发现也是合法集合.
然后我无意见把12加入进去后发现还是合法集合.
于是我们得出结论:等差为奇数的数列都是合法集合,对于连续的那个仅仅只是是等差为1罢了.
于是我们如果 i/x 的个数为 num(i/x),
于是有
则终于答案是
C( num (n / 1) , 2) + C( num (n / 3) , 2) + C( num (n / 5) , 2) + ...
也就是等差位 x 的数列里。我们随便挑一段都是合法集合.
仅仅是到这一步我们发现这样做还是超时.
递推的等差数列
因为F(n+1) 与 F(n) 有非常大的关系,所以我们尝试找找递推行不行.
还是这个公式
F(n) = F(n-1) + f(n).
f(n) 代表含有 n 的合法集合.
然后这些集合须要全是等差数列.
于是写了这么一个公式
f(n-1) =
n/1 + n/3 + n/5 + ....
然后就没什么想法了.
今天看了这个解题报告才知道能够这样做.
发现我们假设写出 f(n-2) 那一项的公式
(n-1)/1 + (n-1)/3 + (n-1)/5 + ...
这两个公式大部分项是相等的,仅仅有个别的几个不相等.
自己举了几个样例发现 n 整除 以下的项时,这个除式会多一个.
比方 n = 12
则 12/3 = 4, 11/3 = 3.
这样这个f(n)函数就能够转化为
f( n ) = f( n - 1) + Count( n -1 ).
当中 Count( n ) 代表 n 的奇数约数个数.
然后小舟学长的模板上刚好有求关于小于 n 的全部数的约数的个数.当中有 O( n*log(n) )
的模板,也有 O( n ) 的模板。
因为 n*log(n) 的模板比較简单,也能够过这道题,于是我就是用 n*log(n) 的模板了.
代码
/*************************************************************************
> File Name: 4.2.cpp
> Author: tiankonguse
> Mail: i@tiankonguse.com
> Created Time: Mon 26 May 2014 01:06:18 PM CST
***********************************************************************/
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<string>
#include<queue>
#include<map>
#include<cmath>
#include<stack>
#include<algorithm>
#include<functional>
#include<stdarg.h>
using namespace std;
#ifdef __int64
typedef __int64 LL;
#else
typedef long long LL;
#endif
const int N = 10000010;
LL nod[N];
LL count[N];
LL ans[N]; void __sieve_nod() {
memset(nod, 0,sizeof(nod));
for (int i = 1; i < N; i+=2) {
for (int j = i; j < N; j += i) {
++nod[j];
}
}
} void init(){
ans[0] = 1;
ans[1] = 2;
LL count = 1;
for(int i=2;i<N;i++){
count += nod[i-1];
ans[i] = ans[i-1] + count;
}
} int main() {
__sieve_nod();
init();
int t,n;
scanf("%d",&t);
for(int i=1; i<=t; i++) {
scanf("%d",&n);
printf("Case #%d:\n%I64d\n",i,ans[n]);
}
return 0;
}
參考
http://acm.hdu.edu.cn/showproblem.php?
pid=4834
【百度之星2014~初赛(第二轮)解题报告】JZP Set的更多相关文章
- 【百度之星2014~初赛(第二轮)解题报告】Chess
声明 笔者近期意外的发现 笔者的个人站点http://tiankonguse.com/ 的非常多文章被其他站点转载.可是转载时未声明文章来源或參考自 http://tiankonguse.com/ 站 ...
- 动态规划-hdoj-4832-百度之星2014初赛第二场
Chess Problem Description 小度和小良近期又迷上了下棋.棋盘一共同拥有N行M列,我们能够把左上角的格子定为(1,1),右下角的格子定为(N,M).在他们的规则中,"王 ...
- 百度之星2014初赛 - 1002 - Grids
先上题目: Grids Time Limit: 10000/5000 MS (Java/Others) Memory Limit: 65535/65535 K (Java/Others)Tota ...
- 百度之星2017初赛A轮 1001 小C的倍数问题
小C的倍数问题 Accepts: 1990 Submissions: 4931 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/3 ...
- 【百度之星2014~复赛)解题报告】The Query on the Tree
声明 笔者最近意外的发现 笔者的个人网站http://tiankonguse.com/ 的很多文章被其它网站转载,但是转载时未声明文章来源或参考自 http://tiankonguse.com/ 网站 ...
- 【百度之星2014~复赛 解题报告~正解】The Query on the Tree
声明 笔者最近意外的发现 笔者的个人网站http://tiankonguse.com/ 的很多文章被其它网站转载,但是转载时未声明文章来源或参考自 http://tiankonguse.com/ 网站 ...
- 百度之星2017初赛A-1005-今夕何夕
今夕何夕 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submis ...
- 百度之星2014复赛 - 1001 - Find Numbers
先上题目: Find Numbers Time Limit: 10000/5000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Othe ...
- 百度之星2014复赛 - 1002 - The Query on the Tree
先上题目: The Query on the Tree Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (J ...
随机推荐
- sql select case when 语句
有道笔试题: 服务器监控表server_status中,当服务器状态发生server_status变化时数据表中将被插入一条记录,状态0表示停机 1表示正常,用SQL查询Server A 的停机开始 ...
- Easy UI下拉列表默认选中(多行)与为文本框赋值
1.为单行文本框赋值 var data2 = $('#LoadArea').combobox("getData"); if (data2) { $('#id).combobox(' ...
- CopyOnWriteArrayList并发容器
CopyOnWriteArrayList并发容器 Copy-On-Write简称COW,是一种用于程序设计中的优化策略.其基本思路是,从一开始大家都在共享同一个内容,当某个人想要修改这个内容的时候,才 ...
- 前端框架:react还是vue?
之前写了一篇前端框架的大汇总,主要介绍了当下主流的框架和其特性.最近除了bootstrap,就属react和vue最为热门,这篇就主要拿这两个框架来做一下详细对比. 究竟如何正确使用?作为小白的我们从 ...
- pt-query-digest
pt-query-digest默认查询时间分布 # Query_time distribution # 1us # 10us ##################################### ...
- input[type='file']样式美化及实现图片预览
前言 上传图片是常见的需求,多使用input标签.本文主要介绍 input标签的样式美化 和 实现图片预览. 用到的知识点有: 1.input标签的使用 2.filelist对象 和 file对象 3 ...
- DLL的导出函数重定向机制
曾经,调试时跟进HeapAlloc,结果发现直接进入到ntdll的RtlAllocateHeap中,感到很有趣,就使用Dependency Walker查看kernel32.dll的导出函数,结果发现 ...
- AIO5系统中-打印样式常见问题处理
1.为什么我的报表没有数据显示? 答:请看一下报表是否有绑定数据源,有添加查询,如果这些都做了,请看下主项数据的数据源有没有绑定,这些都绑定了,就会有数据的显示. 2.为什么我做的报表分组很乱? 答: ...
- redux中间件的原理——从懵逼到恍然大悟
前言react已经出来很久了,其生态圈之庞大,一锅炖不下!各种react-xx,已让我们不堪重负,github上随便一个demo,引入的模块至少都是五指之数+.看着头疼,嚼之无味…….在此建议新学者, ...
- 74、django之ajax补充
之前的ajax使用都是依据jquery来使用的,本篇会先分析ajax的原生的js代码实现,还有jsonp的介绍,与OMR的一些遗漏补充. 本篇导航: js实现的ajax 同源策略与Jsonp 一.js ...