hdu 3652

题目：http://acm.hdu.edu.cn/showproblem.php?pid=3652

目要求包含13且被13整除，设计状态 dp[i,j,k,l]: i位数中第一位是j，是否包含13（k=0或1）,模13余数是l的数有几个。

 

知道怎么处理f数组即可

#include <cstdio>
#include <algorithm>
#include <cstring>

using namespace std;
long long f[11][10][2][13];
int n;
int a[10];
//f[i,j,k,l] 首位为j的i位数，是否含13 (k=0、1)，余数为l 
//k=0为辅助数组 
long long jc(int x)//10^x 
{
	if (x==1)
		return 10;
	if (x==0)
		return 1;
	if (x%2==1)
		return jc(x/2)*jc(x/2)*10;
	else
		return jc(x/2)*jc(x/2);
}

void dp(int x)
{
	int i,j,len=0,temp=0,t=0;
	long long ans=0;
	
	while (x)
	{
		a[++len]=x%10;
		x/=10;
	}
	for (i=len;i>=1;i--)
	{
		for (j=0;j<a[i];j++)
		{
			ans+=f[i][j][1][(13-temp)%13];//首位为j的i位数前面含13的个数 
			if ((a[i+1]==1 && j==3)|| t==1)//当前位为13的个数 
				ans+=f[i][j][0][(13-temp)%13];
		}
		if (t==0 && a[i+1]==1 && a[i]==3)//前面已经有13 
			t=1;
			
		temp=(temp+a[i]*jc(i-1))%13;//求余 
	}
	printf("%lld\n",ans);
}

void baoli()
{
	int n,i,i1,len;
	long long tot=0;
	while (scanf("%d",&n)!=EOF)
	{
		tot=0;
		for (i=13;i<=n;i+=13)
		{
			i1=i;
			len=0;
			while (i1)
			{
				a[++len]=i1%10;
				i1/=10;
				if (a[len]==1 && a[len-1]==3)
				{
					tot++;
					break;
				}
			}
		}
		printf("%d\n",tot);
	}
	return;
}

int main()
{
	int i,j,k,l,l1;
	
	freopen("a.txt","r",stdin);
//	baoli();
	for (i=1;i<=10;i++)//预处理 
	{
		for (j=0;j<=9;j++)
		{
			if (i==1)//1位数 
			{
				f[i][j][0][j]=1;
				continue;
			}
			for (l=0;l<13;l++)
			{
				l1=(l-j*jc(i-1))%13;
				if (l1<0)
					l1+=13;
				for (k=0;k<=9;k++)
				{
					f[i][j][1][l]+=f[i-1][k][1][l1];
					if (j==1 && k==3)
						f[i][j][1][l]+=f[i-1][k][0][l1];
					else
						f[i][j][0][l]+=f[i-1][k][0][l1];	
				}
			}
		}
	}
	while (scanf("%d",&n)!=EOF)
	{
		memset(a,0,sizeof(a));
		dp(n+1);
	}
	
	
	return 0;
}