dp的练习总结(5)

P6187 最小环

1.求 n 和 k 最大公约数，记为 t，那么 t 就是互不相关的环的个数。

2.那么答案最大的情况就是大的放在一个环中，小的放在一个环中，从而实现积的不平均。

或者说，记 n/t=p，那么第 1 到第 p 大的放在一个环中，第p+1 到第 2p 大的放在一个环中，以此类推。注意一下 k=0 的情况和边界。

3.考虑一下优化：到一个数的因子很少，这道题 n 又是固定的，也就是说 t,p 的个数很少，所以只用记忆化就行了。只要 p 相同结果就相同，如果以前算过 p，就直接输出答案。

PS：（记得开long long）

代码呈上：

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int MAXN=2e5+10;
ll a[MAXN],f[MAXN];
int read(){
	int x=0,f=1;
	char c=getchar();
	while(c<'0'||c>'9'){
		if(c=='-') f=-1;
		c=getchar();
	}
	while(c>='0'&&c<='9'){
		x=x*10+c-'0';
		c=getchar();
	}
	return x*f;
}
int gcd(int a,int b){
	if(b==0) return a;
	return gcd(b,a%b);
}
int main(){
	int n=read(),m=read();
	for(int i=1;i<=n;i++) a[i]=read();
	sort(a+1,a+n+1);
	while(m--){
		ll ans=0;
		int k=read();
		if(k==0||n==1){
			for(int i=1;i<=n;i++) ans+=a[i]*a[i];
			printf("%lld\n",ans);
			continue;
		}
		int t=gcd(n,k),p=n/t;
		if(f[p]){
			printf("%lld\n",f[p]);
			continue;
		}
		for(int i=1;i<=n;i+=p){
			for(int j=0;j<p-2;j++) ans+=a[i+j]*a[i+j+2];
			ans+=(a[i]*a[i+1]+a[i+p-1]*a[i+p-2]);
		}
		printf("%lld\n",ans);
		f[p]=ans;
	}
	return 0;
}

P2470 压缩

1.定义 dp[i][j] 表示前 i 个字符最后在第 j-1 到 j 位插入 M 得到的最短长度。

2.由于是求最短长度，于是我们将 dp 数组初始化为极大值，方便之后的 min 操作。边界值 dp[1][1]=1 ，前 1 个字符最短长度一定是 1 。

之后 i 从 1 到字符串长度， dp[i][i] 可初始化为：

for(int j=1;j<i;j++)dp[i][i]=min(dp[i][i],dp[i-1][j]+2);

3.之后 dp[i][j] 状态转移分为两种情况，一是直接插入当前字符，即是 dp[i][j]=dp[i−1][j]+1，而另一种方式首先判断满足"abab"格式，如下：

        l∼(l+r)/2与 (l+r)/2+1∼r 部分完全相等，可另开函数判断；

        (r−l+1) mod 2=0

代码如下：

#include <bits/stdc++.h>
char a[505];
int dp[505][505];
using namespace std;
bool check(int l, int r) {
    int mid = (l + r) / 2;
    for (int i = l, j = mid + 1; i <= mid && j <= r; i++, j++) {
        if (a[i] != a[j]) {
            return false;
        }
    }
    return true;
}
int main() {
    while (scanf("%s", a + 1) != EOF) {
        int n = strlen(a + 1);
        memset(dp, 127, sizeof(dp));
        dp[1][1] = 1;
        for (int i = 1; i <= n; i++) {
            for (int j = 1; j < i; j++) dp[i][i] = min(dp[i][i], dp[i - 1][j] + 2);
            for (int j = 1; j < i; j++) {
                dp[i][j] = dp[i - 1][j] + 1;
                if (check(j, i) && (j - i + 1) % 2 == 0)
                    dp[i][j] = min(dp[i][j], dp[(i + j - 1) / 2][j] + 1);
            }
        }
        int ans = 1e9;
        for (int j = 0; j < n; j++) ans = min(ans, dp[n][j]);
        cout << ans << endl;
    }
    return 0;
}