HDOJ 5793 A Boring Question（快速幂+逆元+数学推导）

最新推荐文章于 2019-04-01 20:02:36 发布

原创最新推荐文章于 2019-04-01 20:02:36 发布 · 675 阅读

0 ·

CC 4.0 BY-SA版权

------题解------ 专栏收录该内容

179 篇文章

订阅专栏

本文介绍了一种结合快速排序和逆元解决特定数学问题的方法。通过实例解析，阐述了如何利用这些技术来提高算法效率并确保计算准确性。同时提供了完整的代码实现。

思路：

据说这是道水题…然后只需要推导出一个数学公式...

表示看了之后觉得推不出来...但是涌神告诉我们说，必须要用到快速排序和逆元。快速排序是为了更快的降低时间复杂度，而逆元则是为了在运算中不会出错。比如，两个int类型的数相乘，可能会爆掉，所以这时候只能用逆元来做。类似于之前做过的一道题，求（a/b）%9973的那道题其实是一个道理。至于说找规律，不妨举个例子，a2=3a1+1，等比为m，那么配凑一下，1+x=m*x，x=1/(m-1)，上面那个例子就凑成了(a2+1/2)/(a1+1/2)=3。所以首项是：m+1+[1/(m-1)]，第n项是[m+1+1/(m-1)]*m^(n-1)，那么答案就是：[m^(n+1)-1]/(m-1)。

代码：

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#define mod 1000000007
using namespace std;

long long powt(long long a,long long b)
{
    long long r = 1;
    while(b)
    {
        if(b & 1) r = r * a % mod;
        a = a * a % mod;
        b >>= 1;
    }
    return r;
}

int main()
{
    long long t,n,m;
    scanf("%lld",&t);
    while(t--)
    {
        scanf("%lld%lld",&n,&m);
        printf("%lld\n",((powt(m,n + 1) - 1) * powt(m - 1,mod - 2) % mod + mod) % mod);
    }
    return 0;
}

打表代码：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <queue>
#include <map>
#include <set>
#include <vector>
#define LL long long
#define mid(a,b) ((a+b)>>1)
#define eps 1e-8
#define maxn 2100
#define mod 1000000007
#define inf 0x3f3f3f3f
#define IN freopen("in.txt","r",stdin);
using namespace std;

LL e[510][510];
void make(){
    for(int i=0;i<510;i++)
        e[i][0]=1;
    for(int i=1;i<510;i++)
        for(int j=1;j<510;j++)
            e[i][j]=(e[i-1][j-1]+e[i-1][j])%mod;
}

int n,m,ans;
void fun(int len, vector<int> cur, int last) {
    if(len == m) {
        int tmp = 1;
        for(int i=1; i<cur.size(); i++) {
            tmp *= e[cur[i]][cur[i-1]];
        }
        ans += tmp;
        return;
    }

    for(int i=last; i<=n; i++) {
        cur.push_back(i);
        fun(len+1, cur, i);
        cur.pop_back();
    }
}

int main(int argc, char const *argv[])
{
    //IN;

    make();

    for(n=0; n<=5; n++) {
        for(m=2; m<=5; m++) {
            ans = 0;
            vector<int> cur; cur.clear();
            fun(0,cur,0);
            printf("%d-%d : %d\n", n,m,ans);
        }
    }

    return 0;
}