Luogu1169BZOJ1057[ZJOI2007]棋盘制作

一、悬线法

原作者博客 https://rpdreamer.blog.luogu.org/p1169

#include <bits/stdc++.h>
using namespace std;
namespace fdata
{
inline char nextchar()
{
    static const int BS = 1 << 21;
    static char buf[BS], *st, *ed;
    if (st == ed)
        ed = buf + fread(st = buf, 1, BS, stdin);
    return st == ed ? -1 : *st++;
}
inline int poread()
{
    int ret = 0;
    char ch;
    while (!isdigit(ch = nextchar()))
        ;

    do
        ret = ret * 10 + ch - '0';
    while (isdigit(ch = nextchar()));
    return ret;
}
} // namespace fdata
using fdata::poread;
const int MAXN = 2005;
int res[MAXN][MAXN], left_[MAXN][MAXN], right_[MAXN][MAXN], up_[MAXN][MAXN];
int n, m, ans1, ans2;
int main()
{
#ifdef lky233
    system("color FD");
    freopen("testdata.in", "r", stdin);
    freopen("testdata.out", "w", stdout);
#endif
    n = poread();
    m = poread();
    for (register int i = 1; i <= n; ++i)
    {
        for (register int j = 1; j <= m; ++j)
        {
            res[i][j] = poread();
            left_[i][j] = right_[i][j] = j;
            up_[i][j] = 1;
        }
    }
    for (register int i = 1; i <= n; ++i)
        for (register int j = 2; j <= m; ++j)
        {
            if (res[i][j] != res[i][j - 1])
                left_[i][j] = left_[i][j - 1];
        }
    for (register int i = 1; i <= n; ++i)
        for (register int j = m - 1; j > 0; --j)
            if (res[i][j] != res[i][j + 1])
                right_[i][j] = right_[i][j + 1];
    for (register int i = 1; i <= n; ++i)
        for (register int j = 1; j <= m; ++j)
        {
            if (i > 1 && res[i][j] != res[i - 1][j])
            {
                left_[i][j] = max(left_[i][j], left_[i - 1][j]);
                right_[i][j] = min(right_[i][j], right_[i - 1][j]);
                up_[i][j] = up_[i - 1][j] + 1;
            }
            int a = right_[i][j] - left_[i][j] + 1;
            int b = min(a, up_[i][j]);
            ans1 = max(ans1, b * b);
            ans2 = max(ans2, a * up_[i][j]);
        }
    printf("%d\n%d", ans1, ans2);
}

思路:预处理左右边界然后计算长度。

二、单调栈

#include <bits/stdc++.h>
using namespace std;
namespace fdata
{
inline char nextchar()
{
    static const int BS = 1 << 21;
    static char buf[BS], *st, *ed;
    if (st == ed)
        ed = buf + fread(st = buf, 1, BS, stdin);
    return st == ed ? -1 : *st++;
}
inline int poread()
{
    int ret = 0;
    char ch;
    while (!isdigit(ch = nextchar()))
        ;

    do
        ret = ret * 10 + ch - '0';
    while (isdigit(ch = nextchar()));
    return ret;
}
} // namespace fdata
using fdata::poread;
const int MAXN = 2005;
int n, m, top, cur, ans1, ans2;
int stack[MAXN], map[MAXN][MAXN], h[MAXN];
int main()
{
#ifdef lky233
    system("color FD");
    freopen("testdata.in", "r", stdin);
    freopen("testdata.out", "w", stdout);
#endif
    n = poread();
    m = poread();
    for (register int i = 1; i <= n; ++i)
        for (register int j = 1; j <= m; ++j)
            ::map[i][j] = poread();
    for (register int j = 1; j <= n; ++j)
    {
        for (register int i = 1; i <= m; ++i)
            if (j > 1 && ::map[j][i] != ::map[j - 1][i])
                ++h[i];
            else
                h[i] = 1;
        cur = 1;
        while (cur <= m)
        {
            ::stack[0] = cur - 1;
            ::stack[top = 1] = cur++;
            while (cur <= m && ::map[j][cur] != ::map[j][cur - 1])
            {

                while (top && h[::stack[top]] > h[cur])
                {
                    int x = h[::stack[top]];
                    int y = cur - ::stack[top - 1] - 1;
                    int z = min(x, y);
                    ans1 = max(ans1, z * z), ans2 = max(ans2, x * y), --top;
                }
                ::stack[++top] = cur++;
            }
            while (top)
            {
                int x = h[::stack[top]];
                int y = cur - ::stack[top - 1] - 1;
                int z = min(x, y);
                ans1 = max(ans1, z * z);
                ans2 = max(ans2, x * y);
                --top;
            }
        }
    }
    printf("%d\n%d\n", ans1, ans2);
    return 0;
}

 

转载于:https://www.cnblogs.com/Shiina-Rikka/p/11263334.html

基于粒子群优化算法的p-Hub选址优化(Matlab代码实现)内容概要:本文介绍了基于粒子群优化算法(PSO)的p-Hub选址优化问题的研究与实现,重点利用Matlab进行算法编程和仿真。p-Hub选址是物流与交通网络中的关键问题,旨在通过确定最优的枢纽节点位置和非枢纽节点的分配方式,最小化网络总成本。文章详细阐述了粒子群算法的基本原理及其在解决组合优化问题中的适应性改进,结合p-Hub中转网络的特点构建数学模型,并通过Matlab代码实现算法流程,包括初始化、适应度计算、粒子更新与收敛判断等环节。同时可能涉及对算法参数设置、收敛性能及不同规模案例的仿真结果分析,以验证方法的有效性和鲁棒性。; 适合人群:具备一定Matlab编程基础和优化算法理论知识的高校研究生、科研人员及从事物流网络规划、交通系统设计等相关领域的工程技术人员。; 使用场景及目标:①解决物流、航空、通信等网络中的枢纽选址与路径优化问题;②学习并掌握粒子群算法在复杂组合优化问题中的建模与实现方法;③为相关科研项目或实际工程应用提供算法支持与代码参考。; 阅读建议:建议读者结合Matlab代码逐段理解算法实现逻辑,重点关注目标函数建模、粒子编码方式及约束处理策略,并尝试调整参数或拓展模型以加深对算法性能的理解。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值