数位dp小练

1.HDU3555 BOMB

我是超链接

题意：找0~n中含有“49”的个数

题解：

状态：f[i][j]表示i位数以j开头的数中不含49的个数 

转移：if (j!=4 || k!=9) f[i][j]+=f[i-1][k]; （不含49，最后用总数减去不含的就是含有的）

因为数据n+1怕爆longlong，手动+1

代码：

#include <cstdio>
#include <cstring>
#define LL long long
using namespace std;
//i位数以j开头的数中不含49的个数 
LL f[50][50];
int T,digi[50];LL n,ans;
void dp()
{
	int i,j,k;
	f[0][0]=1;
	for (i=1;i<=19;i++)
	  for (j=0;j<=9;j++)
	    for (k=0;k<=9;k++)
	      if (j!=4 || k!=9)
		    f[i][j]+=f[i-1][k]; 
}
LL find(LL n)
{
	LL ans=0;int cnt=0,i,j;
	memset(digi,0,sizeof(digi));
	while (n)
	{
		digi[++cnt]=n%10; n/=10;
	}
	digi[1]++;
	for (i=1;i<=cnt;i++)
	  if (digi[i]>=10) digi[i+1]+=digi[i]/10,digi[i]%=10;
	if (digi[cnt+1]) cnt++;
	for (i=cnt;i>=1;i--)  
	  {
	  	for (j=0;j<digi[i];j++)
	  	  if (j!=9 || digi[i+1]!=4) ans+=f[i][j];
	  	if (digi[i]==9 && digi[i+1]==4) break;
	  }
	return ans;
}
int main()
{
	dp();
	scanf("%d",&T);
	while (T--)
	{
		scanf("%I64d",&n);
		ans=find(n);
	    ans=n-ans+1;
	    printf("%I64d\n",ans);
	}	  
}

2、HDU2089 不要62

我是超链接

题意：不要“4” 和“62”

题解：

和上个题差不多，最后用前缀和差分做

代码：

#include <cstdio>
#include <cstring>
#define LL long long
using namespace std;
LL digi[20],f[20][20];
void dp()
{
	int i,j,k;
	f[0][0]=1;
	for (i=1;i<=7;i++)
	  for (j=0;j<=9;j++)
	    for (k=0;k<=9;k++)
	      if (j!=4 && (j!=6 || k!=2))
	        f[i][j]+=f[i-1][k];
}
LL find(int n)
{
	LL ans=0;int cnt=0,i,j;
	memset(digi,0,sizeof(digi));
	while (n)
	{
		digi[++cnt]=n%10; n/=10;
	}
	for (i=cnt;i>=1;i--)  
	  {
	  	for (j=0;j<digi[i];j++)
	  	  if ((j!=2 || digi[i+1]!=6)&&j!=4) ans+=f[i][j];
	  	if ((digi[i]==2 && digi[i+1]==6) || digi[i]==4) break;
	  }
	return ans;
}
int main()
{
	dp();
	int n,m;
	scanf("%d%d",&n,&m);
	while (n!=0 || m!=0)
	{
		LL ans1=find(n);
		LL ans2=find(m+1);
		printf("%I64d\n",ans2-ans1);
		scanf("%d%d",&n,&m);
	}
}

3.HDU4722 Good Numbers

我是超链接

题意：求区间[A,B]内各个数位相加能被10整除的数的个数。

题解：

状态：f(i,j)表示i位数，所有各个数位相加%10==j的数的个数。 

有毒的找规律题目

代码：

#include <cstdio>
#include <cstring>
#define LL long long
using namespace std;
//f(i,j)表示i位数，所有各个数位相加%10==j的数的个数。 
LL f[30][20],m,n;
int T,t,digi[20];
void dp()
{
	int i,j,k;
	f[0][0]=1;
	for (i=1;i<=19;i++)
	  for (j=0;j<=9;j++)
	    for (k=0;k<=9;k++)
	      f[i][j]+=f[i-1][k];
}
LL find(LL n)
{
	LL ans=0;int cnt=0,i,j,last=0;
	memset(digi,0,sizeof(digi));
	while (n)
	{
		digi[++cnt]=n%10; n/=10;
	}
	for (i=cnt;i>=1;i--)  
	{
		for (j=0;j<digi[i];j++)
	  	  ans+=f[i-1][(j+last)%10];
	  	last+=digi[i];
	}  
	return ans;
}
int main()
{
	dp();
	scanf("%d",&T);
	for (int t=1;t<=T;t++)
	{
		scanf("%I64d%I64d",&n,&m);
		printf("Case #%d: %I64d\n",t,find(m+1)-find(n));
	}
}

4、HDU3652 B-number

我是超链接

题意：求区间[0,n]中所有能被13整除并且含“13”的数的个数。

题解  
状态：f(i,j,k)表示i位数，mod 13==j, k=0 含“13” k=1 不含“13”并且开头不为“3” k=2 不含“13”并且开头为“3” 的数的个数 
其实这种状态的表示是比较经典的，并且前几道题都可以用这种方式来dp，但是由于前几道题比较简单，用其他方法更好解决，但是这道题用这种方法就非常优越了。 
转移： amod13+bmod13≡(a+b)mod13(mod13) 进行转移。 
注意：由于前面的和会对后面产生影响，所以要维护一个记录和的量。

代码：

#include <cstdio>
#include <cstring>
using namespace std;
//f(i,j,k)表示i位数，mod 13==j, k=0 含“13” k=1 不含“13”并且开头不为“3” 
//k=2 不含“13”并且开头为“3” 的数的个数 
int n,f[20][20][3],mi[20],digit[20];
void dp()
{
	int i;
    f[0][0][1]=1;
	for (i=0;i<=9;i++)
	{
		f[1][i%13][0]=0;
        if (i!=3) f[1][i%13][1]=1;
        if (i==3) f[1][i%13][2]=1;
	}	
    int a,aMod,bMod;
    for (int i=2;i<=10;++i)
      for (int j=0;j<13;++j){
        for (int k=0;k<10;++k){//k枚举的是这个i位数的第一位是哪个数
          a=k*mi[i-1]; aMod=a%13;
          //(a+b)%13==j a%13=aMod
          bMod=((j-aMod)%13+13)%13;

          f[i][j][0]+=f[i-1][bMod][0];
          if (k==1) f[i][j][0]+=f[i-1][bMod][2];

          if (k!=3){
            f[i][j][1]+=f[i-1][bMod][1];
            if (k!=1) f[i][j][1]+=f[i-1][bMod][2];
          }

          if (k==3) f[i][j][2]+=f[i-1][bMod][1]+f[i-1][bMod][2];
        } 
      }
}
int find(int n){
    int cnt=0,ans=0,a=0,aMod=0,bMod=0,last=0;bool flag=false;
	memset(digit,0,sizeof(digit));
    while (n) digit[++cnt]=n%10,n/=10;
    for (int i=cnt;i>=1;--i){
        //首先加上除去这一位含13的 
        for (int j=0;j<digit[i];++j){
          a=last+mi[i-1]*j;
          aMod=a%13;
          bMod=((0-aMod)%13+13)%13;
          ans+=f[i-1][bMod][0];
        }
        //如果前面已经符合条件，除去这一位所有的都可以加上 
        if (flag){
            for (int j=0;j<digit[i];++j){
                a=last+mi[i-1]*j;
                aMod=a%13;
                bMod=((0-aMod)%13+13)%13;
                ans+=f[i-1][bMod][1]+f[i-1][bMod][2];
            }
            last+=mi[i-1]*digit[i];
            continue;
        }
        //如果上一位是1，这一位>3的话，要加上不含13但是开头为3的
        if (digit[i+1]==1&&digit[i]>3){
            a=last;
            aMod=a%13;
            bMod=((0-aMod)%13+13)%13;
            ans+=f[i][bMod][2];
        }

        //如果这一位>1的话，要加上不含13但是开头为3的
        if (digit[i]>1){
            a=last+mi[i-1];
            aMod=a%13;
            bMod=((0-aMod)%13+13)%13;
            ans+=f[i-1][bMod][2];
        }

        if (digit[i]==3&&digit[i+1]==1) flag=true;
        //最后加上这一位的数 
        last+=mi[i-1]*digit[i];
    }
    return ans;
}
int main()
{
	int n,i;
	mi[0]=1;for (i=1;i<=10;i++) mi[i]=mi[i-1]*10;
	dp();
	while (scanf("%d",&n)!=EOF) printf("%d\n",find(n+1));
}