[DP] ARC074 E - RGB Sequence

本文介绍了一种使用动态规划解决特定计数问题的方法,通过定义状态dp数组来记录方案数,并排除非法组合实现最终答案的计算。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

Solution

dpr,g,b表示考虑前max{r,g,b}个格子,三种颜色最后一次放在r,g,b位置的方案数。
把不合法的(r,g,b)三元组ban掉就好了。

#include <bits/stdc++.h>
using namespace std;

const int MOD = 1000000007;
const int N = 305;
typedef pair<int, int> Pairs;

inline char get(void) {
    static char buf[100000], *S = buf, *T = buf;
    if (S == T) {
        T = (S = buf) + fread(buf, 1, 100000, stdin);
        if (S == T) return EOF;
    }
    return *S++;
}
template<typename T>
inline void read(T &x) {
    static char c; x = 0; int sgn = 0;
    for (c = get(); c < '0' || c > '9'; c = get()) if (c == '-') sgn = 1;
    for (; c >= '0' && c <= '9'; c = get()) x = x * 10 + c - '0';
    if (sgn) x = -x;
}

int n, m, x, y, z, w, ans;
vector<Pairs> v[N];
int dp[N][N][N];
int ban[N][N][N];

inline void Add(int &x, int a) {
    x = (x + a) % MOD;
}
inline void Update(int r, int g, int b, int w) {
    if (ban[r][g][b]) return;
    Add(dp[r][g][b], w);
}
inline void Check(int r, int g, int b) {
    int k = max(r, max(g, b));
    ban[r][g][b] = 0;
    for (auto x: v[k])
        if ((r >= x.first) + (g >= x.first) + (b >= x.first) != x.second)
            return (void)(ban[r][g][b] = 1);
}

int main(void) {
    freopen("1.in", "r", stdin);
    freopen("1.out", "w", stdout);
    read(n); read(m);
    for (int i = 0; i < m; i++) {
        read(x); read(y); read(z);
        v[y].push_back(Pairs(x, z));
    }
    for (int r = 0; r <= n; r++)
        for (int g = 0; g <= n; g++)
            for (int b = 0; b <= n; b++)
                Check(r, g, b);
    dp[0][0][0] = 1;
    for (int r = 0; r <= n; r++)
        for (int g = 0; g <= n; g++)
            for (int b = 0; b <= n; b++)
                if (w = dp[r][g][b]) {
                    int k = max(r, max(g, b));
                    Update(k + 1, g, b, w);
                    Update(r, k + 1, b, w);
                    Update(r, g, k + 1, w);
                    if (k == n) Add(ans, w);
                }
    cout << ans << endl;
    return 0;
}
import cv2 import numpy as np def is_approx_rect(contour, epsilon_factor=0.02): peri = cv2.arcLength(contour, True) approx = cv2.approxPolyDP(contour, epsilon_factor * peri, True) return (4 <= len(approx) <= 5 and cv2.isContourConvex(approx)), approx def calc_center(approx): M = cv2.moments(approx) if M["m00"] == 0: return None return int(M["m10"] / M["m00"]), int(M["m01"] / M["m00"]) def distance(p1, p2): return np.sqrt((p1[0]-p2[0])**2 + (p1[1]-p2[1])**2) def main(): cap = cv2.VideoCapture("222.mp4") if not cap.isOpened(): print("打开视频失败") return prev_center = None while True: ret, frame = cap.read() if not ret: break gray = cv2.cvtColor(frame, cv2.COLOR_BGR2GRAY) _, binary = cv2.threshold(gray, 120, 255, cv2.THRESH_BINARY_INV) closed = cv2.morphologyEx(binary, cv2.MORPH_CLOSE, cv2.getStructuringElement(cv2.MORPH_RECT, (50, 50))) contours_data = cv2.findContours(closed, cv2.RETR_TREE, cv2.CHAIN_APPROX_SIMPLE) contours = contours_data[1] if len(contours_data) == 3 else contours_data[0] candidates = [] for cnt in contours: is_rect, approx = is_approx_rect(cnt) if is_rect: center = calc_center(approx) if center: candidates.append((approx, center, cv2.contourArea(approx))) if not candidates: selected = None elif prev_center is None: selected = max(candidates, key=lambda x: x[2]) else: candidates.sort(key=lambda x: distance(x[1], prev_center)) top_n = [candidates[0]] for c in candidates[1:]: if distance(c[1], prev_center) - distance(candidates[0][1], prev_center) < 50: top_n.append(c) else: break selected = max(top_n, key=lambda x: x[2]) display_frame 将上述代码改成适用于 openmv4 h7 plus 的代码要求给出完整代码
最新发布
08-03
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值