CodeForces 1409F : Subsequences of Length Two DP

这篇博客探讨了一种字符串操作问题,给定一个长度为n的字符串S和长度为2的字符串T,以及整数k,通过最多k次修改将S中的字符转换,求解最多可以得到多少T的子序列。博主使用动态规划的方法,通过状态转移来解决这个问题,最终得出答案。

传送门

题目描述

给定长度为n的串S和长度为2的串T,和一个整数k,
一次操作你可以将S串中任意一个字符修改为其他字符,最多进行k次操作,
问最后S串中最多有多少个子序列T。

分析

我们分类讨论一下第i位字符修改,不修改的情况,令f[i][j][l]表示前i个字母中修改j次,包含有l个op[0]的情况下,答案总数,然后进行状态转移即可

代码

#pragma GCC optimize(3)
#include <bits/stdc++.h>
#define debug(x) cout<<#x<<":"<<x<<endl;
#define dl(x) printf("%lld\n",x);
#define di(x) printf("%d\n",x);
#define _CRT_SECURE_NO_WARNINGS
#define pb push_back
#define mp make_pair
#define all(x) (x).begin(),(x).end()
#define fi first
#define se second
#define SZ(x) ((int)(x).size())
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> PII;
typedef vector<int> VI;
const int INF = 0x3f3f3f3f;
const int N = 220;
const ll mod= 1000000007;
const double eps = 1e-9;
const double PI = acos(-1);
template<typename T>inline void read(T &a){char c=getchar();T x=0,f=1;while(!isdigit(c)){if(c=='-')f=-1;c=getchar();}
while(isdigit(c)){x=(x<<1)+(x<<3)+c-'0';c=getchar();}a=f*x;}
int gcd(int a,int b){return (b>0)?gcd(b,a%b):a;}
char str[N],op[3];
int f[N][N][N];
int n,k;

int main(){
    read(n),read(k);
    scanf("%s %s",str,op);
    memset(f,-0x3f,sizeof f);
    f[0][0][0] = 0;
    int s = (op[0] == op[1]);
    for(int i = 0;i < n;i++){
        int x = (str[i] == op[0]),y = (str[i] == op[1]);
        for(int j = 0;j <= k;j++)
            for(int p = 0;p <= i;p++){
                f[i + 1][j][p + x] = max(f[i + 1][j][p + x],f[i][j][p] + y * p);
                f[i + 1][j + 1][p + 1] = max(f[i + 1][j + 1][p + 1],f[i][j][p] + s * p);
                f[i + 1][j + 1][p + s] = max(f[i + 1][j + 1][p + s],f[i][j][p] + p);
            }
    }
    int ans = 0;
    for(int i = 0;i <= k;i++)
        for(int j = 0;j <= n;j++)
            ans = max(ans,f[n][i][j]);
    di(ans);
    return 0;
}

/**
*  ┏┓   ┏┓+ +
* ┏┛┻━━━┛┻┓ + +
* ┃       ┃
* ┃   ━   ┃ ++ + + +
*  ████━████+
*  ◥██◤ ◥██◤ +
* ┃   ┻   ┃
* ┃       ┃ + +
* ┗━┓   ┏━┛
*   ┃   ┃ + + + +Code is far away from  
*   ┃   ┃ + bug with the animal protecting
*   ┃    ┗━━━┓ 神兽保佑,代码无bug 
*   ┃        ┣┓
*    ┃        ┏┛
*     ┗┓┓┏━┳┓┏┛ + + + +
*    ┃┫┫ ┃┫┫
*    ┗┻┛ ┗┻┛+ + + +
*/



一种基于有效视角点方法的相机位姿估计MATLAB实现方案 该算法通过建立三维空间点与二维图像点之间的几何对应关系,实现相机外部参数的精确求解。其核心原理在于将三维控制点表示为四个虚拟基点的加权组合,从而将非线性优化问题转化为线性方程组的求解过程。 具体实现步骤包含以下关键环节:首先对输入的三维世界坐标点进行归一化预处理,以提升数值计算的稳定性。随后构建包含四个虚拟基点的参考坐标系,并通过奇异值分解确定各三维点在该基坐标系下的齐次坐标表示。接下来建立二维图像点与三维基坐标之间的投影方程,形成线性约束系统。通过求解该线性系统获得虚拟基点在相机坐标系下的初步坐标估计。 在获得基础解后,需执行高斯-牛顿迭代优化以进一步提高估计精度。该过程通过最小化重投影误差来优化相机旋转矩阵和平移向量。最终输出包含完整的相机外参矩阵,其中旋转部分采用正交化处理确保满足旋转矩阵的约束条件。 该实现方案特别注重数值稳定性处理,包括适当的坐标缩放、矩阵条件数检测以及迭代收敛判断机制。算法能够有效处理噪声干扰下的位姿估计问题,为计算机视觉中的三维重建、目标跟踪等应用提供可靠的技术基础。 资源来源于网络分享,仅用于学习交流使用,请勿用于商业,如有侵权请联系我删除!
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值