Skip to main content
  1. Posts/

(BC 一周年) hdu 5312 Sequence

·2 mins
Note: This article is available in Chinese only. 本文暂无英文版本。 View original

比赛的时候没做出来.这道题需要用到的一个重要的性质是,任意一个自然数可以表示成至多三个三角形数(1,3,6,10,15…..)的和(orz高斯)然后也有推广到任意自然数可以表示成k个k角形数的和的结论(费马提出了猜想,柯西给了证明)然后官方题解说的比较好:

这个题看上去是一个贪心, 但是这个贪心显然是错的. 事实上这道题目很简单, 先判断1个是否可以, 然后判断2个是否可以. 之后找到最小的k (k > 2)k(k>2), 使得(m - k) mod 6 = 0(m−k)mod6=0即可.

证明如下: 3n(n-1)+1 = 6(n(n-1)/2)+13n(n−1)+1=6(n∗(n−1)/2)+1, 注意到n(n-1)/2n∗(n−1)/2是三角形数, 任意一个自然数最多只需要3个三角形数即可表示. 枚举需要kk个, 那么显然m=6(km=6(k个三角形数的和)+k)+k, 由于k ge 3k≥3, 只要m-km−k是6的倍数就一定是有解的.**

事实上, 打个表应该也能发现规律.

另外还有一点,特判一个和两个的情况时,一个的好判断,扫一遍就好了

两个的话,由于这个数列是递增的,我们可以从两边往中间,算是一个不错的优化,具体见代码.

 1    /*************************************************************************
 2    	> File Name: code/nv/#ann/1003.cpp
 3    	> Author: 111qqz
 4    	> Email: rkz2013@126.com
 5    	> Created Time: 2015年07月28日 星期二 23时03分09秒
 6     ************************************************************************/
 7
 8    #include<iostream>
 9    #include<iomanip>
10    #include<cstdio>
11    #include<algorithm>
12    #include<cmath>
13    #include<cstring>
14    #include<string>
15    #include<map>
16    #include<set>
17    #include<queue>
18    #include<vector>
19    #include<stack>
20    #define y0 abc111qqz
21    #define y1 hust111qqz
22    #define yn hez111qqz
23    #define j1 cute111qqz
24    #define tm crazy111qqz
25    #define lr dying111qqz
26    using namespace std;
27    #define REP(i, n) for (int i=0;i<int(n);++i)
28    typedef long long LL;
29    typedef unsigned long long ULL;
30    const int inf = 0x7fffffff;
31    const int N=1E5+7;
32    int k,m,f[N];
33    void init()
34    {
35        for ( int i = 1 ; i <N; i++)
36        {
37    	f[i]=3*i*(i-1)+1;
38    	if (f[i]>1000000000)
39    	{
40    	    k = i-1;
41    	    break;
42    	}
43        }
44    }
45    int solve (int x)
46    {
47        for ( int i = 1 ; f[i]<=x ; i++ )
48        {
49    	if (x==f[i])
50    	    return 1;
51        }
52        int j = k;
53        for ( int i = 1 ; i <= k-1&&f[i]<x ; i++)
54        {
55    	while(f[i]+f[j]>x) j--;
56    	if (f[i]+f[j]==x) return 2;
57        }
58        for ( int i = 3 ; i <= m ; i++ )
59        {
60    	if ((m-i)%6==0)
61    	    return i;
62        }
63    }
64    int main()
65    {
66        int T;
67        init();
68        cin>>T;
69        int ans;
70        while (T--)
71        {
72
73    	scanf("%d",&m);
74    	cout<<solve(m)<<endl;
75
76        }
77
78    	return 0;
79    }

Related

codeforces 535 C.Tavas and karafs (解方程)

·2 mins
http://codeforces.com/problemset/problem/535/C 题读了好几遍才读懂。 题意是给出一个等差数列,操作严格要求从最左边不为零的连续m个数减去1,最多执行t次后问离最左边最远的位置在哪里。 有两个限制条件…一个是本身的si不能大于t,否则无法吃完。 还有一个是从sl到sr的和不能超过m*t (比赛的时候考虑的不周到。。实际上只有当r-l+1比m大的时候才是m,也就是说要取min(m,l-r+1)) 这题正解应该是二分….直接Lower_bound。。。看到也有人用前缀和搞的。 我是解方程了(貌似是个傻逼做法)…. 可以列出一个关于r的一元二次方程。。。然后求根公式2333 方程是:

codeforces 534 C Polycarpus' Dice

·1 min
http://codeforces.com/problemset/problem/534/C 题意是说一共有N个骰子,第I个筛子一共有di面…现在知道这些骰子的点数之和,问对于每一个骰子不能取得值有多少个。

cf 534B. Covered Path

·1 min
http://codeforces.com/problemset/problem/534/B 题意是说一辆车,每秒内的速度恒定…第I秒到第I+1秒的速度变化不超过D。初始速度为V1,末速度为V2,经过时间t,问最远能走多远。

hdu 2138 How many prime numbers

ACM STEPS里的…这题前面一道是求LCM….结果接下来就是这么一道。。。 朴素会超….筛法会爆….题目顺序真是按照难度来的? 于是想到 Miller-Rabin素数测试……. 这个方法是基于费马小定理 我的理解就是… 如果我要判断n是否为素数 只要取k个数 如果满足 a^(n-1)mod n =1 那么n就很可能为素数。 证明什么的…暂时还是算了吧…论文里貌似扯了一大堆 第一次用,竟然真的A了。。。。 感觉更好的办法也许是先打一个比较小的素数表,然后每次random选取若干个进行判断…那样应该更可靠些? 本来想WA掉之后再改的。。。没想到这么写就A掉了。。。。杭电数据略水?

三角形数_百度百科

·2 mins
它有一定的规律性,排列如下(构成图),像上面的1、3、6、10、15等等这些能够表示成三角形的形状的总数量的数,叫做三角形数。