BZOJ 2999 inint

从起点1开始，每次选择当前数的任意一位上加上去，问得到n的最小步数以及方案数。多组数据。

例如，从1开始得到100，有很多方法，其中有下面两种方式：
A. 1-2-4-8-16-17-18-19-20-22-24-28-36-39-48-56-62-68-76-83-91-100
B. 1-2-4-8-16-17-24-28-36-39-48-56-62-68-76-83-91-100
显然，B只需要17步。
而事实上，有两种17步的方法。
C. 1-2-4-8-16-22-24-28-36-39-48-56-62-68-76-83-91-100

发现状态只和有哪些1~9的数字有关。
状态：$f[d][S][i][j]$表示对于任意一个满足只存在S状态集合中的数x，我们从$x\ast 10^d+i$转移到$(x+1)\ast 10^d+j$且$(x+1)\ast 10^d+j$是转移过程中第一个大于等于$(x+1)\ast 10^d$的点。
然后可以逐层转移。
最后可以用这个数组求出对于n的答案。
时间复杂度$O(10^5{2}^9)$BZOJ最慢。

AC Code:

#include<bits/stdc++.h>
#define inf 0x3f3f3f3f3f3f3fll
#define mod 1000000007
#define LL long long
using namespace std;

LL n;
struct data{
	LL f;
	int g;
	data(LL f=inf,int g=0):f(f),g(g){}
	inline void check(data B,data C){
		if(f > B.f + C.f) f = B.f + C.f , g = 0;
		if(f == B.f + C.f) g = (g + 1ll * B.g * C.g) % mod;
	} 
}num[10],f[14][512][10][10],t[12][10];

int main(){
	/*
	freopen("inint.in","r",stdin);
	freopen("inint.out","w",stdout);
	*/
	for(int i=1;i<=9;i++) num[i] = data(1,1);
	
	data b[20];
	for(int sta=0;sta<512;sta++)
		for(int i=0;i<=9;i++){
			for(int j=0;j<20;j++) b[j]=data(inf,0);
			b[i] = data(0,1);
			for(int j=0;j<10;j++)
				if(b[j].g)
					for(int k=1;k<=9;k++)
						if(j == k || (sta>>(k-1)&1))
							b[j+k].check(b[j] , num[k]);
			for(int j=0;j<10;j++)
				f[1][sta][i][j] = b[j+10];
		}
		
	for(int bit=2;bit<13;bit++){
		for(int sta=0;sta<512;sta++)
			for(int i=0;i<=9;i++){
				for(int j=0;j<=10;j++)
					for(int k=0;k<10;k++)
						t[j][k] = data(inf,0);
				t[0][i] = data(0,1);
				for(int j=0;j<10;j++)
					for(int k=0;k<10;k++) if(t[j][k].g)
						for(int p=0;p<10;p++)
							t[j+1][p].check(t[j][k] , f[bit-1][sta|(j?(1<<j-1):0)][k][p]);
				for(int j=0;j<10;j++)
					f[bit][sta][i][j]  = t[10][j];
			}
		}
	int T;scanf("%d",&T);
	while(~scanf("%lld",&n)){
		data b[10],c[10];
		for(int i=0;i<10;i++) b[i] = data(inf,0);
		b[1] = data(0,1);
		int S = 0 , S2 = 0;
		for(LL i=1000000000000ll,dep=12;dep;i/=10,dep--){
			int p=0;
			for(;n>=i;n-=i){
				S2 = (S | (p?(1<<p-1):0));
				memcpy(c,b,sizeof b);
				for(int j=0;j<10;j++) b[j] = data(inf,0);
				for(int j=0;j<10;j++) if(c[j].g)
					for(int k=0;k<10;k++)
						b[k].check(c[j],f[dep][S2][j][k]);
				p++;
			}
			S |= ((p?(1<<p-1):0));
		}
		for(int j=0;j<10;j++) if(b[j].g)
			for(int k=1;k<10;k++)
				if(j+k<10 && (j == k || (S>>(k-1)&1)))
					b[j+k].check(b[j],num[k]);
		
		if(b[n%10].g) printf("%lld %d\n",b[n%10].f,b[n%10].g);
		else puts("IMPOSSIBLE");
	}
}