Problem:
ジグザグ数 (Zig-Zag Numbers)

問題

正の整数を (先頭に 0 をつけずに) 10 進法で書いて桁の数字を順番に見ていくと増加と減少を交互に繰り返しているとき，その数は「ジグザグ数」であると呼ぶことにしよう．たとえば，2947 は桁の数字が 2 → 9 → 4 → 7 と，増加 → 減少 → 増加 の順になっているのでジグザグ数である．また，71946 は 減少 → 増加 → 減少 → 増加 の順なのでジグザグ数である．一方，123 や 71446 や 71442 や 88 はジグザグ数ではない．なお，1 桁の正の整数はジグザグ数であると考えることとする．

A 以上 B 以下の M の倍数のうち，ジグザグ数の個数を 10000 で割った余りを求めるプログラムを作成せよ．

入力

入力は 3 行からなり，1 行に 1 つずつ正の整数が書かれている．

1 行目の整数は A を，2 行目の整数は B を，3 行目の整数は M を表す．これらは 1 ≦ A ≦ B ≦ 10500，1 ≦ M ≦ 500 を満たす．

※ A や B の値は，通常の整数を表すデータ型には収まらない可能性があることに注意せよ．

出力

A 以上 B 以下の M の倍数のうち，ジグザグ数の個数を 10000 で割った余りを 1 行で出力せよ．

入出力例

入力例 1

100
200
5

出力例 1

13

入出力例 1 において，100 以上 200 以下の 5 の倍数であるジグザグ数は，105, 120, 130, 140, 150, 160, 165, 170, 175, 180, 185, 190, 195 の 13 個である．

入力例 2

6
1234567
3

出力例 2

246

入出力例 2 において，6 以上 1234567 以下の 3 の倍数であるジグザグ数は 50246 個あるので，それを 10000 で割った余りである 246 を出力する．

問題文と自動審判に使われるデータは、情報オリンピック日本委員会が作成し公開している問題文と採点用テストデータです。

Solution in C (judge verdict: Memory Limit Exceeded)

#include<stdio.h>
#include<string.h>
#define R 10000
int dp[501][11][501][3][3];//[今の桁][今の数字][mod][上下][flag]
int func(int v[],int vn,int m){
	int i,j,k,l,w,x,y,z,ans;
	memset(dp,0,sizeof(dp));
	ans=0;
	if(vn==1){
		for(i=1;i<=v[0];i++)if(i%m==0)ans++;
		return ans;
	}
	for(i=0;i<=v[0];i++){
		if(i==v[0])z=1;
		else z=0;
		if(i==0){
			dp[0][0][0][2][2]=1;
		}
		else{
			dp[0][i][i%m][0][z]++;
			dp[0][i][i%m][1][z]++;
		}
	}
	for(i=1;i<vn;i++){
		for(l=0;l<m;l++){
			for(j=0;j<10;j++){
				if(j==0){
					dp[i][0][0][2][2]=dp[i-1][0][0][2][2];
				}
				if(j!=v[i-1]){
					if(v[i-1]<j)x=1,y=0;
					else if(v[i-1]>j)x=0,y=1;
					if(j==v[i]){
						dp[i][j][(l*10+j)%m][x][1]+=dp[i-1][v[i-1]][l][y][1];
						dp[i][j][(l*10+j)%m][x][1]%=R;
					}
					else if(j<v[i]){
						dp[i][j][(l*10+j)%m][x][0]+=dp[i-1][v[i-1]][l][y][1];
						dp[i][j][(l*10+j)%m][x][0]%=R;
					}
				}
				for(k=0;k<10;k++){
					if(j==k)continue;
					else if(k<j)x=1,y=0;
					else x=0,y=1;
					if(k==0){
						//dp[i][j][(l*10+j)%m][1][0]+=dp[i-1][0][l][2][2];
						//dp[i][j][(l*10+j)%m][1][0]%=R;
						//dp[i][j][(l*10+j)%m][0][0]+=dp[i-1][0][l][2][2];
						//dp[i][j][(l*10+j)%m][0][0]%=R;
						dp[i][j][(l*10+j)%m][2][0]+=dp[i-1][0][l][2][2];
						//dp[i][j][(l*10+j)%m][2][0]%=R;
					}
					dp[i][j][(l*10+j)%m][x][0]+=dp[i-1][k][l][y][0]+dp[i-1][k][l][2][0];
					//dp[i][j][(l*10+j)%m][x][0]%=R;
				}
				dp[i][j][(l*10+j)%m][2][0]%=R;
				dp[i][j][(l*10+j)%m][x][0]%=R;
				//if(dp[i][j][(l*10+j)%m][1][0]>0)printf("dp[%d][%d][%d][1][0]=%d\n",i,j,(l*10+j)%m,dp[i][j][(l*10+j)%m][1][0]);
				//if(dp[i][j][(l*10+j)%m][0][0]>0)printf("dp[%d][%d][%d][0][0]=%d\n",i,j,(l*10+j)%m,dp[i][j][(l*10+j)%m][0][0]);
				//if(dp[i][j][(l*10+j)%m][1][1]>0)printf("dp[%d][%d][%d][1][1]=%d\n",i,j,(l*10+j)%m,dp[i][j][(l*10+j)%m][1][1]);
				//if(dp[i][j][(l*10+j)%m][0][1]>0)printf("dp[%d][%d][%d][0][1]=%d\n",i,j,(l*10+j)%m,dp[i][j][(l*10+j)%m][0][1]);
				//printf("\n");
			}
		}
	}
	for(i=0;i<10;i++){
		for(j=0;j<2;j++){
			//if(dp[vn-1][i][0][2][j]>0)printf("dp[%d][%d][0][2][%d]=%d\n",vn-1,i,j,dp[vn-1][i][0][2][j]);
			//if(dp[vn-1][i][0][1][j]>0)printf("dp[%d][%d][0][1][%d]=%d\n",vn-1,i,j,dp[vn-1][i][0][1][j]);
			//if(dp[vn-1][i][0][0][j]>0)printf("dp[%d][%d][0][0][%d]=%d\n",vn-1,i,j,dp[vn-1][i][0][0][j]);
			ans+=dp[vn-1][i][0][2][j]+dp[vn-1][i][0][1][j]+dp[vn-1][i][0][0][j];
			ans%=R;
		}
	}
	return ans;
}
int main(){
	int i,j,k,l,m,x[510],y[510],z,xn,yn,ans,c;
	char a[510],b[510];
	
	scanf("%s%s",a,b);
	scanf("%d",&m);
	xn=strlen(a);
	yn=strlen(b);
	for(i=0;i<xn;i++)x[i]=a[i]-'0';
	for(i=0;i<yn;i++)y[i]=b[i]-'0';
	c=1;
	for(i=xn-1;i>=0;i--){
		x[i]=x[i]-c;
		c=0;
		if(x[i]<0){
			c=1;
			x[i]+=10;
		}
	}
	
	ans=(10000+func(y,yn,m)-func(x,xn,m))%R;
	printf("%d\n",ans);
	return 0;
}