据说这是道水题…然后只需要推导出一个数学公式...
表示看了之后觉得推不出来...但是涌神告诉我们说,必须要用到快速排序和逆元。快速排序是为了更快的降低时间复杂度,而逆元则是为了在运算中不会出错。比如,两个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;
}