Ural 1586 Threeprime Numbers

题意：构造一个长度为n 的数，其中任意取三个连续的数，保证这三个连续的数是一个三位数的素数，求长度为n 的数有多少种可能。

链接：http://acm.timus.ru/problem.aspx?space=1&num=1586

思路：首先预处理凷100-1000以内所有的素数，保存在num[j][k][l]中，再进行动态规划，设计状态dp[i][j][k]，表示长度为i时，第一位为j，第二位为k的数种类数，状态转移方程：dp[i][x][y] = ( dp[i-1][y][z] + dp[i][x][y] ) % mod

注意点：当数只有三位时候的状态数需要初始化

以下为AC代码：

ID	Date	Author	Problem	Language	Judgement result	Test #	Execution time	Memory used
6257899	19:57:32 21 May 2015	luminous11	1586. Threeprime Numbers	G++ 4.9 C++11	Accepted		0.062	9 202 KB

/* 
***********************************************
*# @Author  : Luminous11 (573728051@qq.com)
*# @Date    : 2015-05-21 22:02:02
*# @Link    : http://blog.csdn.net/luminous11
*********************************************** 
*/

#include <bits/stdc++.h>
#define clr(a, v) memset( a , v , sizeof(a) )
using namespace std;
const int mod = 1e9 + 9;
const double eps = 1e-10;
int dp[10005][15][15];
bool prime[3005];
bool vis[10][10][10];
int main()
{
//	ios::sync_with_stdio ( false );
	int n;
	memset ( dp, 0, sizeof ( dp ) );
	memset ( vis, 0, sizeof ( vis ) );
	memset ( prime, 1, sizeof ( prime ) );
	for ( int i = 2; i < 1005; i ++ ){
		if ( prime[i] ){
			for ( int j = i + i; j < 1005; j += i ){
				prime[j] = 0;
			}
		}
	}
	int x, y, z;
	for ( int i = 100; i < 1000; i ++ ){
		if ( prime[i] ){
			x = i / 100;
			y = i / 10 % 10;
			z = i % 10;
			vis[x][y][z] = 1;
		}
	}
	for ( int x = 1; x < 10; x ++ ){
		for ( int y = 0; y < 10; y ++ ){
			for ( int z = 0; z < 10; z ++ ){
				if ( vis[x][y][z] ){
					dp[3][x][y] ++;
				}
			}
		}
	}
	for ( int i = 4; i < 10005; i ++ ){
		for ( int x = 1; x < 10; x ++ ){
			for ( int y = 0; y < 10; y ++ ){
				for ( int z = 0; z < 10; z ++ ){
					if ( vis[x][y][z] ){
						dp[i][x][y] = ( dp[i-1][y][z] + dp[i][x][y] ) % mod;
					}
				}
			}
		}
	}
	while ( scanf ( "%d", &n ) != EOF ){
		int ans = 0;
		for ( int i = 1; i < 10; i ++ ){
			for ( int j = 0; j < 10; j ++ ){
				ans += dp[n][i][j];
				ans %= mod;
			}
		}
		printf ( "%d\n", ans );
	}
    return 0;
}