HDU3652 含13並且能被13整除

 

待完成

#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 


//計算10的x次方 
long long jc(int 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];
			if ((a[i+1]==1 && j==3)|| t==1)
				ans+=f[i][j][0][(13-temp)%13];
		}
		if (t==0 && a[i+1]==1 && a[i]==3)
		{
			t=1;
		}
//		printf("%d\n",jc(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位數，全部不含13 
				//餘數爲j,餘數即爲j，滿足條件有1個 
				f[i][j][0][j]=1;  //f[i,j,k,l] 首位爲j的i位數，是否含13 (k=0、1)，餘數爲l 
				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];  //1表示含有13 
					}
					else
					{
						f[i][j][0][l]+=f[i-1][k][0][l1];	//0表示不含13 
					}
				}
			}
		}
	}
	while (scanf("%d",&n)!=EOF)
	{
		memset(a,0,sizeof(a));
		dp(n+1);
	}
	
	
	return 0;
}