整数拆分

//整数分拆：给定一个数，问他拆分个数是多少。拆分是一系列整数，他们的和为给定的数
//复杂度O(N*sqrt(N))
//HDU 4651
//将一个数用一个或多个正整数的无序和来表示
//分拆数中全部是奇数的个数和两两不相同的个数是相同的

/*
一些有关限制分拆的结论：
    n的分拆数中最大部分为m的个数=把n分拆成m部分的个数
    n的分拆数中每部分的数都相等的个数=n 的因子个数
        Eg. 6=2+2+2, 6=3+3,6=1+1+1+1+1+1
    n的分拆数中每部分都是1或2(或者把n分拆成1或2部分)的个数=floor(n/2+1);
        Eg. 6=1+1+1+1+1+1, 6=1+1+1+1+2, 6=1+1+2+2, 6=2+2+2
    n的分拆数中每部分都是1或2或3(或者把n分拆成1或2或3部分)的个数=(n+3)^2/12;
*/ 


#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std ;

const int MOD = 1e9 + 7 ;
const int maxn = 1e5 + 10  ;

int dp[ maxn ]  ;

void init(){
    memset(dp , 0 , sizeof( dp )) ;
    dp[0] = 1 ;
    for( int i = 1 ; i<= 100000 ; i++){
        for( int j = 1 , r = 1 ; i - ( 3* j * j - j ) / 2 >= 0  ; j ++ , r *= -1 ){
            dp[i] += dp[ i - ( 3 * j * j - j ) / 2 ] * r ;
            dp[i] %= MOD ;
            dp[i] = ( dp[i] + MOD )%MOD ;
            if( i - ( 3 * j * j + j ) / 2 >= 0 ){
                dp[i] += dp[i - ( 3 * j * j + j ) / 2] * r ;
                dp[i] %= MOD ;
                dp[i] = ( dp[i] + MOD ) %MOD ;
            }
        }
    }
}

int main(){
    int T , n ;
    init() ;
    scanf("%d" , & T) ;
    while( T -- ){
        scanf("%d" , & n ) ;
        printf("%d\n" , dp[n]) ;
    }
    return 0 ;

} 

/*
 * Author:  kuangbin
 * Created Time:  2013/8/8 11:53:35
 * File Name: 1004.cpp
 */
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <string>
#include <vector>
#include <stack>
#include <queue>
#include <set>
#include <time.h>
using namespace std;
const int MOD = 1e9+7;
int dp[100010];
void init()
{
    memset(dp,0,sizeof(dp));
    dp[0] = 1;
    for(int i = 1;i <= 100000;i++)
    {
        for(int j = 1, r = 1; i - (3 * j * j - j) / 2 >= 0; j++, r *= -1)
        {
            dp[i] += dp[i -(3 * j * j - j) / 2] * r;
            dp[i] %= MOD;
            dp[i] = (dp[i]+MOD)%MOD;
            if( i - (3 * j * j + j) / 2 >= 0 )
            {
                dp[i] += dp[i - (3 * j * j + j) / 2] * r;
                dp[i] %= MOD;
                dp[i] = (dp[i]+MOD)%MOD;
            }

        }
    }
}

int solve(int n,int k)
{
    int ans = dp[n];
    for(int j = 1, r = -1; n - k*(3 * j * j - j) / 2 >= 0; j++, r *= -1)
    {
        ans += dp[n -k*(3 * j * j - j) / 2] * r;
        ans %= MOD;
        ans = (ans+MOD)%MOD;
        if( n - k*(3 * j * j + j) / 2 >= 0 )
        {
            ans += dp[n - k*(3 * j * j + j) / 2] * r;
            ans %= MOD;
            ans = (ans+MOD)%MOD;
        }

    }
    return ans;
}

int main()
{
    init();
    int T;
    int n,k;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&n,&k);
        printf("%d\n",solve(n,k));
    }
    return 0;
}