这套题做的蛋疼菊紧
250
简单题。 问一个数能否被表示 成 某个素数的若干次方
我用了一个很损精度得法
立即学习“前端免费学习笔记（深入）”；
其实只要判平方完了直接枚举素数就OK
vector<int>ans;bool check(int x) { int m = (int)sqrt(x * 1.0) + 1; if(x == 2) return true; for(int i = 2; i <= m; i++) { if(x % i == 0) return false; } return true;}void gao(long long x) { int x1 = -1, x2 = -1; for(int i = 2; i < 60; i++) { int f = (int)(pow((double)x, 1.0 / i) + eps); long long tmp = 1; for(int j = 0; j < i; j++) tmp = tmp * (long long)f; if(tmp == x) { if(check(f)) { x1 = f, x2 = i; } } } if(x1 != -1) { ans.push_back(x1); ans.push_back(x2); }}500
区间DP
题目意思是说，给一个A串，一个B串
都是只包含0和1，然后用一些列reverse操作，将A变成B
reverse(i,j)表示把i,j这个区间反转
然后这系列操作有个限制
就是进行完一个操作之后，下一个操作必须在这个操作的区间中的子区间中进行，每个操作都是如此
然后这肯定是方便进行区间DP的
看有人写了一个很暴力的DFS， 没敢尝试，感觉复杂度没法算
dp[k][i][j][0]代表a串i位置开始长度为k的子串 不翻转 变成b串j位置开始长度为k的子串 需要的步数
dp[k][i][j][1]代表a串i位置开始长度为k的子串 翻转 变成b串j位置开始长度为k的子串 需要的步数
int n = a.size(); memset(dp, 0x3f, sizeof(dp)); for(int j = 0; j <= n; j++) for(int k = 0; k <= n; k++) dp[0][j][k][0] = dp[0][j][k][1] = 0; for(int i = 1; i <= n; i++) for(int j = 0; j + i <= n; j++) for(int k = 0; k + i <= n; k++) { if(a[j] == b[k]) { dp[i][j][k][0] = min(dp[i][j][k][0], dp[i - 1][j + 1][k + 1][0]); } if(a[j + i - 1] == b[k + i - 1]) { dp[i][j][k][0] = min(dp[i][j][k][0], dp[i - 1][j][k][0]); } if(a[j] == b[k + i - 1]) { dp[i][j][k][1] = min(dp[i][j][k][1], dp[i - 1][j + 1][k][1]); } if(a[j + i - 1] == b[k]) { dp[i][j][k][1] = min(dp[i][j][k][1], dp[i - 1][j][k + 1][1]); } dp[i][j][k][0] = min(dp[i][j][k][0], dp[i][j][k][1] + 1); dp[i][j][k][1] = min(dp[i][j][k][1], dp[i][j][k][0] + 1); } return dp[n][0][0][0] >= 1000 ? -1: dp[n][0][0][0];1000
这题公式很简单
n*(1/n+1/(n - 1) + 1/ (n - 2) +...+ 1/(n - k + 1) )
关键问题来了
n ,k都巨大
然后发现这个是个调和级数求和
数字大的时候只有近似公式
那么试试呗
(1/n+1/(n - 1) + 1/ (n - 2) +...+ 1/(n - k + 1) ) 约等于 log(n + 1) + R
R是欧拉常数
完了k大的时候就用这个公式去搞。不然直接for了
但是wa出翔了
最后发现别人这么干的 本来求出来的公式是log((n + 1) / (n - k + 1))
然后有个函数叫log1p ，是干什么的呢 log1p(x)返回的就是log(x + 1)
但是问题来了，当x巨小的时候,log1p的精度比较高，用log的时候x+1就丢精度了
然后就凑呗，凑着用log1p还不行，分母减个0.5，就是用来调控精度的。
这给我蛋疼的。
完了发现房里好多不用log1p的， 我全给cha掉了
double expectedBuy(string n, string k) { long long x = gao(n); long long y = gao(k); double ans = 0; long long s = x - y + 1; long long mx = 10000000; while(s <= mx) { ans += 1.0 / s; if(s == x) return x * ans; s++; } ans += log1p((double)(x - s + 1) / (s - 0.5)); return ans * x; }
