SRM 572 DIV 2 1000

DP+数学  恶心死我了 DP那部分没什么 数学的那部分各种细节各种繁琐呀

在比赛中就可以做出来这种题的人果然不一般 自己还需锻炼呀

代码：

#include<iostream>

#include<cstdio>

#include<cstring>

#include<cmath>

#include<algorithm>

#include<vector>

#include<set>

#include<map>

#include<string>

#include<queue>

#include<stack>

#include <iomanip>

using namespace std;

#define LL long long

const double eps=1e-6;

const int INF=0x3f3f3f3f;

const int N=52;

const int M=1501;

const long long MOD=1000000007;

long long sum[N][N][M];

long long b[N];

bool prime(int n)

{

    for(int i=2;i*i<=n;++i)

    if(n%i==0)

    return false;

    return true;

}

long long C(long long n,int m)

{

     int a[N],I=0;

     for(int i=2;i<=m;++i)

     if(prime(i))

     a[I++]=i;

     for(int i=1;i<=m;++i)

     b[i]=n-m+i;

     for(int i=2;i<=m;++i)

     {

         int k=i;

         for(int j=0;j<I&&k>1;++j)

         if(k%a[j]==0)

         {

             k=k/a[j];

             for(int l=1;l<=m;++l)

             if(b[l]%((long long)(a[j]))==0)

             {b[l]=b[l]/((long long)(a[j]));break;}

             --j;

         }

     }

     long long tmp=(long long)(1);

     for(int i=1;i<=m;++i)

     tmp=(tmp*(b[i]%MOD))%MOD;

     return tmp;

}

class DistinctRemainders

{

    public :

    int howMany(long long n, int m)

    {

        memset(sum,0,sizeof(sum));

        for(int i=0;i<=m;++i)

        sum[i][0][0]=(long long)(1);

        int a=0;

        for(int i=1;i<=m;++i)

        {

            for(int j=1;j<=i;++j)

            {

                for(int l=0;l<=a;++l)

                {

                    sum[i][j][l]=sum[i-1][j][l];

                    if(l-(i-1)>=0)

                    sum[i][j][l]=(sum[i][j][l]+sum[i-1][j-1][l-(i-1)])%MOD;

                }

            }

            a=a+i;

        }

        long long ans=0;

        long long lm=(long long)(m);

        for(int j=1;j<=m;++j)

        for(int l=(int)(n%lm);l<M&&((long long)(l)<=n);l+=m)

        {

            if(sum[m][j][l]>0)

            {

            	long long h=(n-(long long)(l))/lm;

            	if(h<0) continue;

                long long tmp=C(h+j-1,j-1);

                for(int w=1;w<=j;++w)

                tmp=(tmp*(long long)(w))%MOD;

                tmp=(tmp*sum[m][j][l])%MOD;

                ans=(ans+tmp)%MOD;

            }

        }

        return (int)(ans);

    }

};