D. Magic Gems（矩阵快速幂）

 题意：每颗魔法宝石可以分解成 m 颗普通宝石，魔法宝石和普通宝石都占据 1 体积的空间，但普通宝石不能再被分解。想要使一些魔法宝石分解，使得所有宝石占据的空间恰好为 n单位体积。显然，一个魔法宝石分解后会占据 m 体积空间，不分解的魔法宝石仍占据 1体积空间。想要求出有多少种分解方案，可以让最后得到的宝石恰好占据 n 单位体积。两种分解方案不同当且仅当分解的魔法宝石数量不同，或者是所用的宝石的编号不同。

题解：令f(n,0/1）=长度为n的数量，最后一位为0或1.

f(n,0)=f(n-m,0)+f(n-m,1);

f(n,1)=f(n-1,0)+f(n-1,1);

所以得到的状态转移方程：f(n)=f(n,0)+f(n,1)=f(n-m)+f(n-1).

最后一个魔法宝石没有分解的方案数+最后一个魔法宝石分解的方案数

下面就是构建矩阵

#include<bits/stdc++.h>
using namespace std;
#define int long long
int n,m;
const int N=105;
int f[N][N],a[N][N];
const int mod=1e9+7;
int c[N][N];
void mul(int c[],int a[],int b[][N]){
	int s[N]={0};
	for(int i=0;i<=N-1;i++){
		for(int j=0;j<=N-1;j++){
			s[i]=(s[i]+a[j]*b[j][i])%m;
		}
	}
	memcpy(c,s,sizeof(s));
}
void mul(int c[N][N],int a[N][N]){
	int p[N][N]={0};    
	for(int i=1;i<=m;i++){
		for(int j=1;j<=m;j++){
			for(int k=1;k<=m;k++){
			p[i][j]=(p[i][j]+c[i][k]*a[k][j])%mod;				
			}
		}
	}
	memcpy(a,p,sizeof(p));
}
signed main(){
	cin>>n>>m;
	if(n<m) {
		cout<<"1\n";
		return 0;
	}
	for(int i=2;i<=m;i++) f[i][1]=a[i][i-1]=1;
	f[1][1]=2,a[1][1]=a[1][m]=1,n-=m;
	while(n){
		if(n&1) mul(a,f);
		mul(a,a); 
		n/=2;
	} 
	cout<<f[1][1]<<"\n";
}