HDU 5051 Fraction 本福特定律

ACMmaxx

于 2014-09-27 22:10:36 发布

阅读量1.7k

点赞数

CC 4.0 BY-SA版权

分类专栏：数学和杂题

本文链接：https://blog.youkuaiyun.com/ACMmaxx/article/details/39619435

该博客主要探讨了HDU 5051 Fraction题目，涉及本福特定律的应用。博主指出，当数列a[]满足ai = b*q^i，且q不为特定值时，数列符合本福特定律，可以用来计算n作为a[]前缀的概率。文章表达了对出题者使用非常见知识点的调侃。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

【题目大意】

给定n，b，q。对于无限长数列的a[]，ai = b*q^i。问n是ai的前缀的概率。

【思路】

本题需要的知识点只有一个：本福特定律。（p.s. 我真想打死出题人，平时出出物理题就够了，现在出了个从来没有听过的本福特定律。。。）

只要q不为1,10,100,1000，就可以认为a[]数列是“一堆从实际生活得出的数据”。。。

//#pragma comment(linker, "/STACK:102400000,102400000")
#include<cstdio>
#include<cstring>
#include<vector>
#include<queue>
#include<cmath>
#include<cctype>
#include<string>
#include<algorithm>
#include<iostream>
#include<ctime>
#include<map>
#include<set>
using namespace std;
#define MP(x,y) make_pair((x),(y))
#define PB(x) push_back(x)
typedef __int64 LL;
//typedef unsigned __int64 ULL;
/* ****************** */
const int INF=100011122;
const double INFF=1e100;
const double eps=1e-8;
//const LL mod=1000000007;
const int NN=100010;
const int MM=2000010;
/* ****************** */

double fun(LL n,LL x)
{
    char s1[50], s2[50];
    int i;
    sprintf(s1, "%I64d", n);