2021/01/17训练总结

前言

放假&开始训练的第一天。寒假应该会每天都会训练了，记录一下一天收获和写一些题解，现在洛谷1160题，看寒假结束能有多少吧。今天主要是写写了两题数位dp然后把牛客小白赛31补完了。(我会记住上海站面对数位dp时的无力)

P4317 花神的数论题

题目链接：P4317 花神的数论题
题目大意：

数据范围:$1\le N\le 10^{15}$
题解：我们考虑去计算$1-N$中$sum(i)=x$的个数，然后累乘起来就行了。现在问题就是如何求出$sum(i)=x$的$i$的个数。我们采取记忆化搜索(数位dp)的方法来进行。重要讲解一下$dfs$函数,参数里面$cur$表示剩余位数,$up$表示是否顶着上界，$now$表示现在已经遍历出来的1的个数，$que$表示要求的1的个数。
AC代码:

#include<bits/stdc++.h>

#define ld long double
#define ll long long
using namespace std;
template<class T>
void read(T& x)
{
	T res = 0, f = 1; char c = getchar();
	while (!isdigit(c)) {
		if (c == '-')f = -1; c = getchar();
	}
	while (isdigit(c)) {
		res = (res << 3) + (res << 1) + c - '0'; c = getchar();
	}
	x = res * f;
}
const ll N = 200000 + 10;
const int mod = 10000007;
#define int long long
ll fpow(int x, int y)
{
	int ans = 1;
	while (y)
	{
		if (y & 1)ans = 1ll * ans * x % mod;
		x = 1ll * x * x % mod; y >>= 1;
	}
	return ans;
}
int w[55];
int f[55][2][55][55];
int dfs(int cur, int up, int now, int que)
{
	if (!cur)
	{
		return now == que;
	}

	if (~f[cur][up][now][que])return f[cur][up][now][que];
	int tp = up ? w[cur] : 1;
	int ans = 0;
	for (int i = 0; i <= tp; i++)
	{
		ans = (ans + dfs(cur - 1, up && i == tp, now + (i == 1), que));
	}
	return f[cur][up][now][que] = ans;
}
signed main()
{
	//ios::sync_with_stdio(false);
#ifndef ONLINE_JUDGE
	freopen("test.in", "r", stdin);
#endif // ONLINE_JUDGE
	ll n;
	read(n);
	while (n)
		w[++w[0]] = n & 1, n >>= 1;
	ll ans = 1;
	memset(f, -1, sizeof(f));
	for (int i = 1; i <= w[0]; i++)
	{
		ans = ans * fpow(i, dfs(w[0], 1, 0, i)) % mod;
	}
	printf("%lld\n", ans);
	return 0;
}

P4127 [AHOI2009]同类分布

题目链接:P4127 [AHOI2009]同类分布
题目大意：

数据范围:$1\le a\le b\le 10^{18}$
题解：很明显的数位dp，基本上数位dp都是要求$1-n$中有多少个数满足xx要求，我们用$find(n)$函数来代替，这题答案就是$find(b)-find(a-1)$,我们考虑$find(n)$这个函数要怎么求，一般来说数位dp的搜索函数都要有这几个参数,$cur$表示还剩多少位，$limit$表示是否顶着上界，然后有几个参数用来判断是否满足条件和剪枝。这题里面$sum$表示搜索的数的数位和是多少,$rem$表示搜索的数模$mod$后的值，其中$mod$表示我们要求的数位和。
AC代码：

#include<bits/stdc++.h>

#define ld long double
#define ll long long
using namespace std;
template<class T>
void read(T& x)
{
	T res = 0, f = 1; char c = getchar();
	while (!isdigit(c)) {
		if (c == '-')f = -1; c = getchar();
	}
	while (isdigit(c)) {
		res = (res << 3) + (res << 1) + c - '0'; c = getchar();
	}
	x = res * f;
}
const ll N = 200 + 10;
#define int long long
ll f[20][163][163], w[N];
int mod;
ll dfs(int cur, int sum, int rem, int limit)
{
	if (!cur)
	{
		return (mod == sum) && (rem == 0);
	}
	if (!limit && ~f[cur][sum][rem])return f[cur][sum][rem];
	if (sum > mod || sum + 9 * cur < mod)return 0;
	int tp = limit ? w[cur] : 9;
	ll ans = 0;
	for (int i = 0; i <= tp && sum + i <= mod; i++)
	{
		ans += dfs(cur - 1, sum + i, (rem * 10 + i) % mod, limit && (i == tp));
	}
	return f[cur][sum][rem] = ans;
}
ll find(ll x)
{
	w[0] = 0;
	while (x)
		w[++w[0]] = x % 10, x /= 10;
	ll ans = 0;
	for (mod = 1; mod <= w[0] * 9; mod++)
	{
		memset(f, -1, sizeof(f));
		ans += dfs(w[0], 0, 0, 1);
	}
	return ans;
}
signed main()
{
	//ios::sync_with_stdio(false);
#ifndef ONLINE_JUDGE
	freopen("test.in", "r", stdin);
#endif // ONLINE_JUDGE
	ll a, b;
	read(a), read(b);
	printf("%lld\n", find(b) - find(a - 1));

	return 0;
}

小白赛31A

题目链接:A|B
题目大意：

数据范围:$t<=1e5$,$1\le a,x\le 10^9$
题解:题目一看还以为数位dp,在看$t$是$1e5$数量级，明显排除数位dp。我们来考虑，没有第一个限制的情况,假设a是$10101010$,b只要对于a的二进制为1的位填上0即可，而a的二进制为0的位，b可以选择0|1,所以答案就是$2^{a二进制数中0的个数}$。我们现在再加上限制，假设x是$10101010$,对于我们假设当前遍历选择的b前缀都和x相同，例如我们现在遍历到了第4位，那么我们选择的b前缀都是$1010$。我们考虑x每一个为1的位，我们新增答案，我们可以假设b在当前位置为0，那么我们后面的位就相当于没有x的限制了，因为无论我们怎么选择都不会超过x。所以新增$2^{a剩余二进制数中0的个数}$。这样我们遍历了$0\sim x-1$,再加上x并且去除0,就可以得到答案
AC代码:

#include<bits/stdc++.h>

#define ld long double
#define ll long long
using namespace std;
template<class T>
void read(T& x)
{
	T res = 0, f = 1; char c = getchar();
	while (!isdigit(c)) {
		if (c == '-')f = -1; c = getchar();
	}
	while (isdigit(c)) {
		res = (res << 3) + (res << 1) + c - '0'; c = getchar();
	}
	x = res * f;
}
//#define int long long
const ll N = 200000 + 10;
const int mod = 1e9 + 7;
ll t, a, x;

signed main()
{
	//ios::sync_with_stdio(false);
#ifndef ONLINE_JUDGE
	freopen("test.in", "r", stdin);
#endif // ONLINE_JUDGE
	read(t);
	while (t--)
	{

		read(a), read(x);
		ll ans = 0;
		for (int i = 31; i >= 0; i--)
		{
			if (((1ll << i) & x) == 0)continue;
			int c = 0;
			for (int j = i - 1; j >= 0; j--)
				if (((1ll << j) & a)==0)c++;
			ans += (1ll << c);
			if (((1ll << i) & a))break;
		}
		if ((a | x) == (a + x))ans++;
		printf("%lld\n", ans - 1);

	}

	return 0;
}

小白赛31E 解方程

题目链接:解方程
题目大意：

数据范围：$1\le t\le 10^3$,$1\le a,b\le 10^6$
题解:其实之前就遇见过类似的，不过还是没记住这种题目的套路，一般看见x*y这种式子就要想着如何去进行因式分解，这题就是两边加上$a\ast b$得到$$a\ast x+b\ast y+a\ast b=x\ast y+a\ast b$$,移项得到$$a\ast b=x\ast y-a\ast x-b\ast y+a\ast b$$即$$a\ast b=(x-a)\ast (y-b)$$答案就是求$a\ast b$的因子数。暴力$(O\sqrt{n})$求解肯定是不行的，拿唯一分解定理搞一下就行
AC代码:

#include<bits/stdc++.h>

#define ld long double
#define ll long long
using namespace std;
template<class T>
void read(T& x)
{
	T res = 0, f = 1; char c = getchar();
	while (!isdigit(c)) {
		if (c == '-')f = -1; c = getchar();
	}
	while (isdigit(c)) {
		res = (res << 3) + (res << 1) + c - '0'; c = getchar();
	}
	x = res * f;
}
const ll N = 2000000 + 10;
const int mod = 1e9 + 7;
int t, a, b;
int prime[N];
bool fprime[N];
void init(int n)
{
	for (int i = 2; i <= n; i++)
	{
		if (!fprime[i])prime[++prime[0]] = i;
		for (int j = 1; j <= prime[0] && i * prime[j] <= n; j++)
		{
			fprime[prime[j] * i] = 1;
			if (i % prime[j] == 0)break;
		}
	}
}
int solve(ll x)
{
	int ans = 1;
	for (int j = 1; j <= prime[0] && prime[j] * prime[j] <= x; j++)
	{
		int cnt = 0;
		while (x % prime[j] == 0)x /= prime[j], cnt++;
		ans *= cnt + 1;
	}
	if (x != 1)ans *= 2;
	return ans;
}
int main()
{
	//ios::sync_with_stdio(false);
#ifndef ONLINE_JUDGE
	freopen("test.in", "r", stdin);
#endif // ONLINE_JUDGE
	init(1000000);
	read(t);
	while (t--)
	{
		read(a), read(b);
		
		printf("%d\n", solve(a*b));
	}

	return 0;
}