hdu4427 Math Magic dp

点击打开链接~~

给出N,M,K

就是求能满足以下条件的的方法数，（1，2）与（2，1）是两种方式

lcm为最小公倍数

gcd为最大公约数

lcm（i，j）=i＊j／gcd（i，j）

1.  SUM (A ₁ , A ₂ , ..., A _i , A _i+1 ,..., A _K ) = N
2.  LCM (A ₁ , A ₂ , ..., A _i , A _i+1 ,..., A _K ) = M

dp[i][j][k]=方法数   

i  ---- i个数

j  ---- 和为j

k ---- 最小公倍数为k

有一个常识，如果i为x的因子，则i与x最小公倍数一定为x

所以首先预先处理lcm（i，j），再把x的因子预先求出来，这样就可以在这些因子中枚举1～1000之内最多的因子数很小，所以可以得到优化

#include 
   
   
    
    
#include 
    
    
     
     
#include 
     
     
      
      
#include 
      
      
       
       
using namespace std;
const int N=1005;
const int mod=1e9+7;
int dp[2][1005][1005];
int lcm[N][N];
int fat[N];
int gcd(int x,int y)
{
    if(y==0)
        return x;
    return gcd(y,x%y);
}
int main()
{
    int n,m,k;
    for(int i=1;i
       
       
         n) continue; dp[e^1][tem][lcm[ll][fat[h]]]+=dp[e][j][ll]; int tt=dp[e^1][tem][lcm[ll][fat[h]]]%=mod; //printf("%d %d %d\n",tem,lcm[ll][fat[h]],tt); } } e^=1; } printf("%d\n",dp[e][n][m]); } }