【省内训练2018-10-26】网友数

【思路要点】

  • 首先,当 k ≥ 9 k≥9 k9 ,任何 ≥ 28 ≥28 28 的数 x x x 均为网友数。 ( 1 ) (1) (1)
  • ( 1 ) (1) (1) 的证明:
    x ≤ 1000 x≤1000 x1000 ,我们可以暴力验证命题的正确性。
    否则,考虑 x x x 的末位,我们可以用 9 9 9 4 4 4 或是 7 7 7 造出 0 0 0 9 9 9 以内的任何末位,这样,我们可以在 x x x 中将造出的数减去, x x x 将减少一位。
  • 假设 k ≤ 8 k≤8 k8 ,考虑如何判断一个数 x x x 是否为网友数。
  • 考虑枚举每一位的网友数分布,记 d p i , j , k dp_{i,j,k} dpi,j,k 表示考虑了 x x x 最高的 i i i 位,余数为 j j j ,已用了 k k k 个网友数的状态是否能够达到。转移时首先枚举下一位使用的网友数的数量 k ′ k' k ,应当保证 k ′ ≥ k k'≥k kk ,再枚举 4 4 4 的数量 f o u r four four ,从而算出 7 7 7 的数量 s e v e n = k ′ − f o u r seven=k'-four seven=kfour ,新的余数 j ′ j' j 即为 j ∗ 10 + a i + 1 − 4 ∗ f o u r − 7 ∗ s e v e n j*10+a_{i+1}-4*four-7*seven j10+ai+14four7seven 。最后判断是否能够达到 i i i 为位数, j = 0 j=0 j=0 的状态即可。注意到 j j j 每过一位就会 ∗ 10 *10 10 ,因此 j > 6 j>6 j>6 的状态均为无效状态,可以剪枝。
  • 考虑原题,我们显然需要数位 d p dp dp ,但是直接数位 d p dp dp 会重复计算一个数。因此,我们不妨把前面的 d p dp dp 作为状态压入最终的数位 d p dp dp 中。即将上述的 d p i , ∗ , ∗ dp_{i,*,*} dpi,, 作为一个整体计入状态,这大约是 7 ∗ 9 = 63 7*9=63 79=63 0 / 1 0/1 0/1 位,看似状态量十分庞大。
  • 但实际上,从初始状态出发,经过有限步合法转移,能够到达的状态是十分有限的,取 k = 8 k=8 k=8 时,能够到达的状态仅有 1654 1654 1654 种,因此,直接按照前述 d p dp dp 即可。
  • 时间复杂度 O ( 10 ∗ C n t ∗ ∣ V ∣ + 10 ∗ C n t L o g C n t ) O(10*Cnt*|V|+10*CntLogCnt) O(10CntV+10CntLogCnt) ,其中 C n t = 1654 Cnt=1654 Cnt=1654 ∣ V ∣ |V| V 表示 L , R L,R L,R 的位数。
  • 我们甚至可以将本题改为多组询问,可以得到 O ( Q ∗ ∣ V ∣ + 10 ∗ C n t ∗ ∣ V ∣ + 10 ∗ C n t L o g C n t ) O(Q*|V|+10*Cnt*|V|+10*CntLogCnt) O(QV+10CntV+10CntLogCnt) 的时间复杂度。

【代码】

#include<bits/stdc++.h>
using namespace std;
const int MAXN = 2e5 + 5;
typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
template <typename T> void chkmax(T &x, T y) {x = max(x, y); }
template <typename T> void chkmin(T &x, T y) {x = min(x, y); } 
template <typename T> void read(T &x) {
	x = 0; int f = 1;
	char c = getchar();
	for (; !isdigit(c); c = getchar()) if (c == '-') f = -f;
	for (; isdigit(c); c = getchar()) x = x * 10 + c - '0';
	x *= f;
}
template <typename T> void write(T x) {
	if (x < 0) x = -x, putchar('-');
	if (x > 9) write(x / 10);
	putchar(x % 10 + '0');
}
template <typename T> void writeln(T x) {
	write(x);
	puts("");
}
namespace force {
	const int MAXN = 38;
	bool valid[MAXN];
	vector <int> good;
	ll calc(ll r) {
		ll ans = 0;
		for (int i = 1; i <= 37 && i <= r; i++)
			ans += valid[i];
		if (r >= 37) ans += r - 37;
		return ans;
	}
	void work(ll l, ll r) {
		valid[0] = true;
		for (int i = 1; i <= 9; i++) {
			for (int j = 37; j >= 0; j--)
				if (valid[j]) {
					if (j + 4 <= 37) valid[j + 4] = true;
					if (j + 7 <= 37) valid[j + 7] = true;
				}
		}
		writeln(calc(r) - calc(l - 1));
	}
}
ll ql, qr, k;
namespace DynamicProgramming {
	const int states = 1654;
	const int MAXS = 1800;
	const int MAXN = 25;
	unordered_map <bitset <128>, int> mp;
	bitset <128> trans(bitset <128> now, int nxt) {
		bitset <128> res; res.reset();
		for (int i = 0; i < 128; i++) {
			if (!now[i]) continue;
			int Residue = i / 10, Used = i % 10;
			for (int used = Used; used <= k; used++)
			for (int four = used; four >= 0; four--) {
				int residue = Residue * 10 + nxt - four * 4 - (used - four) * 7;
				if (residue < 0) break;
				if (residue <= 6) res.set(residue * 10 + used);
			}
		}
		return res;
	}
	int l, r, tot, bits, a[MAXN], edge[MAXS][10];
	ll dp[MAXN][2][MAXS]; bitset <128> q[MAXS];
	ll getans(int pos, bool type, int state) {
		if (pos == 0) {
			bool ans = false;
			for (int i = 0; i <= k; i++)
				ans |= q[state][i];
			return ans;
		}
		if (dp[pos][type][state] != -1)
			return dp[pos][type][state];
		ll ans = 0;
		for (int i = 0; i <= (type ? a[pos] : 9); i++)
			ans += getans(pos - 1, type && (i == a[pos]), edge[state][i]);
		return dp[pos][type][state] = ans;
	}
	ll work(ll r) {
		bits = 0;
		while (r != 0) {
			a[++bits] = r % 10;
			r /= 10;
		}
		memset(dp, -1, sizeof(dp));
		return getans(bits, 1, 1);
	}
	void main() {
		l = r = tot = 1, q[1] = 1;
		mp[1] = tot++;
		while (l <= r) {
			bitset <128> now = q[l++];
			for (int i = 0; i <= 9; i++) {
				bitset <128> dest = trans(now, i);
				if (!mp.count(dest)) {
					mp[dest] = tot++;
					q[++r] = dest;
				}
			}
		}
		//cerr << r << endl;
		for (int i = 1; i <= r; i++)
		for (int j = 0; j <= 9; j++)
			edge[i][j] = mp[trans(q[i], j)];
		writeln(work(qr) - work(ql - 1));
	}
}
int main() {
	read(ql), read(qr), read(k);
	if (k >= 9) force :: work(ql, qr);
	else DynamicProgramming :: main();
	return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值