部分转载 携程A聪明的猴子

以下是poj1091的题解 转载

这题考虑到了数论和组合数学的相关知识，题目大意是给定给定两个整数n和m, 求出长度为n+1的满足条件的数列data的个数，数列要求如下：

1)1 <= data[i] <= m, for 1<=i<=n

2)data[n+1] = m;

3)这个n+1个数满足: 存在整数 x1, x2, ..., xn, xn + 1, 满足 x1 * data[1] + x2 * data[2] + ... + xn+1 * data[n + 1] = 1;

根据数论的知识，n个数的最大公约数规定为这n个数线性和的最小自然数，所以此题就是要求最大公约数为1的数列的个数

由于data[n + 1] = m, 是已知的，所以这题可以用组合数学鸽巢原理的知识来解，步骤如下：

1）求出m的所有质约数，假设有k个如果这n + 1个数的最大公约数不是1，m至少含有一个质约数，这个质约数也是前n个数的约数

2）序列的总个数为 m ^ n;

3）设t(n)表示 能够整除 任意n个m的质约数且<=m的数的个数

利用容斥原理可知公约数不为1的数列的个数为 f = t(1) - t(2) + t(3) - ... + (-1) ^ (k - 1) t(k)

4) 符合要求的序列个数为 m ^ n - f

由于题目数据量小，还可以用dp

 d[i][j]表示i个数最大公约数为j的数的数目

枚举第i+1个数k，那么d[i][gcd(j,k)]=sum{d[i-1][j]}

代码是别人的，比较简短

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
long long dp[30][30];
int gcd(int a,int b){
    return b?gcd(b,a%b):a;
}
int main()
{
    int T,a,b,i,j,k;
    scanf("%d",&T);
    while(T--){
        memset(dp,0,sizeof(dp));
        scanf("%d%d",&a,&b);
        for (i=1;i<=b;i++) dp[1][i]=1;
        for (i=1;i<a;i++){
            for (j=1;j<=b;j++){
                for (k=1;k<=b;k++){
                    dp[i+1][gcd(j,k)]+=dp[i][j];
                }
            }
        }
        long long ans=0;
        for (i=1;i<=b;i++){
            if (gcd(i,b)==1) ans+=dp[a][i];
        }
        printf("%I64d\n",ans);
    }
    return 0;
}