Code forces 93E Lostborn （dp+递归）*

题目链接：http://codeforces.com/problemset/problem/93/E

#include<bits/stdc++.h>
using namespace std;

#define debug puts("YES");
#define rep(x,y,z) for(int (x)=(y);(x)<(z);(x)++)
#define ll long long

#define lrt int l,int r,int rt
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define root l,r,rt
#define mst(a,b) memset((a),(b),sizeof(a))
const int  maxn =100+5;
const int mod=1e9+7;
const int ub=2e5;

ll powmod(ll x,ll y){ll t; for(t=1;y;y>>=1,x=x*x%mod) if(y&1) t=t*x%mod; return t;}
ll gcd(ll x,ll y){return y?gcd(y,x%y):x;}
/*
题意很简单，
给定一个n，和k个数字，
求在1到n这个区间内，不被k个数字任意一个整除的数字个数。

这道题刚开始一看以为是，，，弱智容斥题，
但交了几次全部超时的我不得不反思。。。
看题解才知道还有dp的这一套路，对范围较小的n做记忆化，
对范围较大的n直接递归求解。
复杂度上面应该是记忆化dp更胜一筹，
虽然我最初的容斥我自己算了下扩展曾应该不多，
但是稍微复杂起来复杂度也不能保证多稳，
dp的思想就是，dp（i，j）代表前j个数中不被ai，ai+1.。。等数所整除的
个数，那么dp（1，n）就是所求答案。
dp递推式子就是：dp（i，j）=dp（i+1，j）-dp（i+1，j/a[i]).

时间复杂度：O（能过）。
*/

ll n,k,ans;
ll a[maxn];
ll dp[maxn][ub];

ll Solve(int m,ll n)
{
    if(m>k||n==0) return n;
    if(n<=ub&&dp[m][n]!=-1) return dp[m][n];
    ll ans=Solve(m+1,n)-Solve(m+1,n/a[m]);
    if(n<ub) dp[m][n]=ans;
    return ans;
}

bool cmp(ll x,ll y){return x>y;}

int main(){
    scanf("%lld%lld",&n,&k);
    for(int i=1;i<=k;i++) scanf("%lld",&a[i]);
    sort(a+1,a+k+1,cmp);
    mst(dp,-1);
    printf("%I64d\n",Solve(1,n));
    return 0;
}