HDOJ HDU 1141 Factstone Benchmark

Gy_Blog

于 2017-09-12 22:43:42 发布

阅读量337

点赞数

分类专栏： -----------HDOJ 题目----------- 数论题目文章标签： HDOJ

本文链接：https://blog.csdn.net/Gy_1096440940/article/details/77951617

版权

-----------HDOJ 题目----------- 同时被 2 个专栏收录

56 篇文章 0 订阅

订阅专栏

数论题目

7 篇文章 0 订阅

订阅专栏

HDOJ 1141 Factstone Benchmark

题目

点此查看 HDOJ 1141 Factstone Benchmark

分类

数论对数化简

题意

amtel每10年芯片位数翻一倍，把位数和阶乘做对比作为Facfstone基准
举例 1960 年位数为4, 2^4 < 4!, 2^4 > 3! 所以1960-1969年Facfstone 为 3

题解

设y为(年数/10 - 196)，设Facfstone 基准为f
即：对比 2^(4 * 2^y) 和 f!
由于指数函数和阶乘增长太快，可以用对数函数化简
$\ln2^{4 \times 2^y} = 4 \times 2^y \times \ln(2)$
$\ln f! = \ln1 + \ln2 + ... + \ln f$
对比化简结果即可

表

提供表用于大家查看对拍

1960-1969 : 3
1970-1979 : 5
1980-1989 : 8
1990-1999 : 12
2000-2009 : 20
2010-2019 : 34
2020-2029 : 57
2030-2039 : 98
2040-2049 : 170
2050-2059 : 300
2060-2069 : 536
2070-2079 : 966
2080-2089 : 1754
2090-2099 : 3210
2100-2109 : 5910
2110-2119 : 10944
2120-2129 : 20366
2130-2139 : 38064
2140-2149 : 71421
2150-2159 : 134480
2160-2169 : 254016

代码

#include <iostream>
#include <cmath>

#define fmax 1000000
#define smax 216 - 196 + 1

using namespace std;

double f[fmax];
double s[smax];

int main()
{
    int y;
    double t,sum;
    f[0] = 0;
    for(int i = 1;i < fmax;i++)
    {
        sum += i;
        f[i] = f[i-1] + log(i);
    }
    sum = 4;
    for(int i = 0;i < smax;i++)
    {
        t = sum * log(2);
        for(int j = 1;j < fmax;j++)
            if(t < f[j])
            {
                s[i] = j - 1;
                cout << (196+i) * 10 << '-' << (196+i) * 10 + 9 << " : " << s[i] << endl;
                break;
            }
        sum *= 2;
    }
    while(cin >> y && y)
    {
        y /= 10;
        y -= 196;
        cout << s[y] << endl;
    }
    return 0;
}