1423 J Bubble Cup hypothesis

1423 J Bubble Cup hypothesis

题面传送门

Solution 1 : DP

这个做法不想多讲，毕竟是过不去的。但是可能可以用于更广的$c_i$取值。

大概就是$d{p}_{i,j}$表示$c_i$需要是$j$的方案数，其中$c_i$可以分配一些到后面处理。

时间复杂度$O({\log }_{2}m\ast K)$的。$K$表示$c_i$要$\le K$。

不过这题有$5e5$个数据，且卡$1s$，所以不太好卡过去。（但是仍然有不少神仙用dp的这个做法过了，具体的可以看看code里$900ms+$的代码）

这里放上我的巨慢代码(4s+)

using namespace std;
const int INF=0x3f3f3f3f;
typedef pair<int,int> mp;
/*}
*/
int vis[15][62];
int dp[15][62];
const int MOD=1e9+7;
LL m;
int t;
inline int calc(int j,int i){
	if(j==0){
		return dp[i][j]=(i<=7);
	}
	if(vis[i][j]==t) return dp[i][j];
	vis[i][j]=t;
	dp[i][j]=0;
	int addition=((m>>(j-1))&1ll);
	if(i==0){
		dp[i][j]=calc(j-1,addition);
	}
	else{
		dp[i][j]=calc(j,i-1);
		if(i>7){
			int k=i-7-1;
			if(k+k+addition<=14)
			dp[i][j]+=MOD-calc(j-1,k+k+addition);
			if(dp[i][j]>=MOD) dp[i][j]-=MOD;
		}
		int k=i;
		if(k+k+addition<=14)
		dp[i][j]+=calc(j-1,k+k+addition);
		if(dp[i][j]>=MOD) dp[i][j]-=MOD;
	}
	return dp[i][j];
}
int main(){
	scanf("%d",&t);
	while(1){
		scanf("%I64d",&m);
		printf("%d\n",calc(61,0));
		t--;
		if(t==0) break;
	} 
	return 0;
}

Solution 2 : 数学

$P(x)=c_0+c_{1}x...+c_n{x}^n+...=(c_0+8\ast c_3+64\ast c_6...)+2\ast (c_1+8\ast c_4...)+4\ast (c_2+8\ast c_5)$

$P(x)=x+2\ast y+4\ast z$

一个数$x$与$c_0+8\ast c_1+...8^k\ast c_k+...$是一一对应的（$c_0{c}_1{c}_2..$）是$x$的8进制。

方程$x+2y=m$的解的个数=$\lfloor \frac{m}{2}\rfloor +1$

则
$$Ans=\sum _{i=0}^{\lfloor \frac{m}{4}\rfloor }(\lfloor \frac{m-4i}{2}\rfloor +1)=\sum _{i=0}^{\lfloor \frac{m}{4}\rfloor }(\lfloor \frac{m}{2}\rfloor -2z+1)=(\lfloor \frac{m}{4}\rfloor +1)(\lfloor \frac{m}{2}\rfloor -\lfloor \frac{m}{4}\rfloor +1)$$
这样就可以$O(1)$算了。