ABC300E
题面
思路
对于投骰子,最后一步可能得到1、2、3、4、5、6,那么对应的最后一步之前的数是n/1、n/2、n/3、n/4、n/5,n/6。并且每个数字(1、2、3、4、5、6)得到的概率是一样的,即为1/6。
也就是F(n)=(1/6)(F(n/1)+F(n/2)+F(n/3)+F(n/4)+F(n/5)+F(n/6))*{只有n%4==0,F(n/4)才能产生贡献,其余同理}
移项得:F(n)=1/5(F(n/1)+F(n/2)+F(n/3)+F(n/4)+F(n/5)),如果一位数组能开下,就可以直接线性dp进行状态转移,但是这个n<=1e18,所以用记忆化搜索的方式实现dp,记录mp[x]为得到x的概率,由于在取模过程中进行了/5操作,所以求一下5在 mod 99824435意义下的乘法逆元即可
细节:mp[x]在调用之前如果不存在会创建mp[x]=0,对记忆化搜索的过程有影响,所以记忆化x是否保存结果的时候应该使用mp.count(x)
代码;
ABC300F
题面
#include using namespace std;#define FAST ios::sync_with_stdio(false), cin.tie(0), cout.tie(0)#define PII pair#define de(a) cout << #a << " = " << a << "\n";#define deg(a) cout << #a << " = " << a <>= 1;}res %= mod;return res;}map mp;int dp(int n, int inv){if (mp.count(n))return mp[n];int ans = 0;for (int i = 2; i > n;mp[1] = 1;cout <> t;while (t--)solve();return 0;}
思路
先抛出结论:起点一定在1-N某个点(假设字符串第一个字符的下标为1)
证明:首先T是由M个S串拼接而成,假设最优方案具有从第N+1或者更晚的某个字符开始的最长连续的o,然后对于第i个字符,我们决定使用第i-N个字符去替代它(他们是相同的,在不考虑修改的情况下,而且这种情况必定存在,因为i>=N+1)。
得到这个结论之后,就可以O(N)枚举起点,那怎么得到终点呢?
??暴力枚举,从i(1<=i<=N)开始,枚举j(N<=j<=NM)直到i-j这个区间内的x的个数大于k停止,此时答案即为j-i+1 (i-j范围内的x最多出现k次)。思路没错,但是这样跑直接T飞。
写到这一步就可以二分了,sum[i]表示从1-i这个区间内x的个数是多少。只需要求出1-N的即可,因为后面的NM-N个字符,每N个都与前1-N个字符相同,那对于一个大于N的下标i可以这样得到1-i中的x的个数
int f(int x, int n, vector &rw) // 返回从1-x有多少个x{//x为终点下表,n即为上文的N,也就是S串的长度,rw是前缀和数组,rw[i]代表1-i有多少个xint res = (x / n) * rw[n];int rem = x % n;res += rw[rem];return res;}
总结:
枚举起点i(1<=i<=N),二分终点j(i<=j<=NM),[i,j]这个闭区间中x的个数<=k
代码:
#include using namespace std;#define FAST ios::sync_with_stdio(false), cin.tie(0), cout.tie(0)#define PII pair#define de(a) cout << #a << " = " << a << "\n";#define deg(a) cout << #a << " = " << a << " ";#define endl "\n"#define int long long#define LL long longconst int mod = 1e9 + 7;const int N = 1e6 + 5;int dx[4] = {1, 0, -1, 0}, dy[4] = {0, 1, 0, -1};int f(int x, int n, vector &rw) // 返回从1-x有多少个x{int res = (x / n) * rw[n];int rem = x % n;res += rw[rem];return res;}void solve(){int n, m, k;cin >> n >> m >> k;string s;cin >> s;vector rw(n + 1, 0);for (int i = 0; i < n; i++){rw[i + 1] = rw[i];if (s[i] == 'x')rw[i + 1]++;}int res = 0;for (int i = 1; i <= n; i++){int fbeg = f(i - 1, n, rw); // 1-i-1位置有多少xint l = i, r = n * m;while (l > 1;if (f(mid, n, rw) - fbeg <= k){// f(mid, n, rw) - fbeg从i-mid有多少xl = mid + 1;}else{r = mid - 1;}}res = max(r - i + 1, res);}cout << res <> t;while (t--)solve();return 0;}
ABC300G
参考:AtCoder Beginner Contest 300——A-G题讲解atcoder比赛阿史大杯茶的博客-CSDN博客
官方题解:Editorial – UNIQUE VISION Programming Contest 2023 Spring(AtCoder Beginner Contest 300)
#include using namespace std;#define FAST ios::sync_with_stdio(false), cin.tie(0), cout.tie(0)#define PII pair#define de(a) cout << #a << " = " << a << "\n";#define deg(a) cout << #a << " = " << a << " ";#define endl "\n"#define int long long#define LL long longconst int mod = 1e9 + 7;const int N = 1e6 + 5;int dx[4] = {1, 0, -1, 0}, dy[4] = {0, 1, 0, -1};int n, p;void push(vector &a, int num) // 将与num结合合法的数全部加进来{int sz = a.size();for (int i = 0; i n)break;a.push_back(t);}}}void solve(){vector prime = {2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97};cin >> n >> p;while (p < prime.back()) // 将大于p的质因数删去prime.pop_back();vector frt = {1}, bck = {1};for (auto &c : prime){if (frt.size() < bck.size())push(frt, c); // 采用这种方式降低push函数里面的时间复杂度elsepush(bck, c);}sort(frt.begin(), frt.end());sort(bck.begin(), bck.end());int res = 0;for (int i = 0, j = bck.size() - 1; i = 0 && left < bck[j])j--;if (j < 0)break;res += j + 1;}cout << res <> t;while (t--)solve();return 0;}