数位dp:bzoj1799
Posted qq936584671
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了数位dp:bzoj1799相关的知识,希望对你有一定的参考价值。
一开始没有想到枚举所有模数求解,卡了很久。
写法见代码注释。
#include<iostream> #include<cstdio> #include<cstring> #include<algorithm> #include<cmath> #include<stack> #include<map> #include<vector> #include<queue> #include<set> #include<iomanip> #include<cctype> #include<ctime> using namespace std; #define ll long long #define edl putchar(‘ ‘) #define sscc ios::sync_with_stdio(false);cin.tie(0);cout.tie(0); #define FOR(i,a,b) for(int i=a;i<=b;i++) #define ROF(i,a,b) for(int i=a;i>=b;i--) #define FORLL(i,a,b) for(ll i=a;i<=b;i++) #define ROFLL(i,a,b) for(ll i=a;i>=b;i--) #define mst(a) memset(a,0,sizeof(a)) #define mstn(a,n) memset(a,n,sizeof(a)) #define zero(x)(((x)>0?(x):-(x))<eps) inline ll gcd(ll a,ll b){return b?gcd(b,a%b):a;} inline ll lcm(ll a,ll b){return a/gcd(a,b)*b;} int month[15]={0,31,28,31,30,31,30,31,31,30,31,30,31}; const int MAXN=1e5+5; const int INF=1<<30; const long long mod=1e9+7; const double eps=1e-8; ll dp[25][170][170]; /*dp数组记录答案,第一维为第几位,第二维为当前数字模p的结果,第三维为数位和*/ int a[25],p,vis[25][170][170]; /*a数组记录每一位的值,p为模数,vis数组对应dp,代表在模p时是否访问过*/ ll dfs(int pos,int m,int sta,int limit)/*第几位 模后的数 数位和 上限*/ { if(pos==0) return sta==0&&m%p==0; if(!limit&&vis[pos][m][sta]==p) { return dp[pos][m][sta]; } vis[pos][m][sta]=p;/*更新vis访问标记*/ int up=limit?a[pos]:9; ll ans=0; FOR(i,0,up) { if(sta>=i) ans+=dfs(pos-1,(m*10+i)%p,sta-i,limit&&i==up); /*由于枚举p时数位和是固定的,所以直接用p减去数位和*/ } if(!limit) dp[pos][m][sta]=ans; return ans; } ll solve(ll x) { int pos=0,t; while(x) { a[++pos]=x%10; x/=10; } a[pos+1]=0; t=pos*9; ll ans=0; FOR(i,1,t) p=i,ans+=dfs(pos,0,i,1); return ans; } int main() { ll n,m; { scanf("%lld%lld",&n,&m); mst(vis); n=solve(n-1); mst(vis);/*每次处理都要对vis初始化*/ m=solve(m); printf("%lld ",m-n); } return 0; }
以上是关于数位dp:bzoj1799的主要内容,如果未能解决你的问题,请参考以下文章
bzoj 1799 [Ahoi2009]self 同类分布(数位DP)
BZOJ 1799: [Ahoi2009]self 同类分布 ( 数位dp )