Sicily 3394 Pair

本文探讨了一道有趣的图论问题,并提供了三种不同的解决方案,包括类似拓扑排序的方法、深度优先搜索(DFS)以及匈牙利算法。通过这些方法展示了如何有效地解决特定类型的图论问题。

http://soj.me/3394

一道很有趣的题目,居然写了三种做法。可以用类似拓扑排序的思想,先把入度为一的点入队,每次pop出2个,并且把图中所有指向第二个点的边删掉.如果删掉边的点入度为一了则入队.


#include <vector>
#include <queue>
#include <algorithm>
#include <iostream>
#include <cstdio>
#include <cstring>

using namespace std;

#define FILL(x) memset(x, 0, sizeof(x));
#define PB push_back
const int M = 10010;

int main()
{
    int T;
    cin>>T;
    while(T--) {
        vector<int> g[M];
        int f, t;
        int nn;
        int deg[M];
        FILL(deg);
        scanf("%d", &nn);
        for (int i=0; i<nn-1; i++) {
            scanf("%d%d", &f, &t);
            g[f].PB(t);
            g[t].PB(f);
            deg[f] ++;
            deg[t] ++;

        }
        queue<int> q;
        bool vis[M];
        FILL(vis);
        for (int i=1; i<=nn; i++) {
            if (deg[i] == 1) 
                q.push(i);
        }
        while(!q.empty()) {
            int f = q.front(); q.pop();
            if (deg[f] == 1 && !vis[f]) {
                int t = g[f][0];
                vis[t] = vis[f] = 1;
                g[f].clear();
                for (int i=0; i<g[t].size(); i++) {
                    int v = g[t][i];
                    for (int j=0; j<g[v].size(); j++) {
                        int x = g[v][j];
                        if (x == t) {
                            g[v].erase(g[v].begin() + j);
                            deg[v] --;
                        }
                    }
                    if (deg[v] == 1) 
                        q.push(v);
                }
            }
        }
        bool flag = 1;
        for (int i=1; i<=nn; i++) {
            if (vis[i] != 1)
                flag = 0;
        }
        if (flag) cout<<"Yes\n";
        else cout<<"No\n";
    }
}                                 
好吧,看错了,其实是水题. 一次dfs
#include <vector>
#include <queue>
#include <algorithm>
#include <iostream>
#include <cstdio>
#include <cstring>

using namespace std;
const int M = 10010;
int vis[M];
vector<int> g[M];
int dfs(int f);

int main()
{
    int T;
    cin>>T;
    while(T--) {
        for (int i=0; i<M; i++) {
            g[i].clear();
        }
        int f, t, nn;
        scanf("%d", &nn);
        memset(vis, 0, sizeof(vis));
        for (int i=0; i<nn-1; i++) {
            scanf("%d%d", &f, &t);
            g[f].push_back(t);
            g[t].push_back(f);
        }
        if (dfs(1) == 0) 
            cout<<"Yes"<<endl;
        else 
            cout<<"No"<<endl;
    }
}
int dfs(int f) {
    vis[f] = 1;
    int cnt = 0;
    int flag = 0;
    for (int i=0; i<g[f].size(); i++) {
        int t=g[f][i];
        if (!vis[t])  {
            flag ++;
            cnt += dfs(t);
        }
    }
    if (flag == 0) return 1; // leaf
    if (cnt >= 2) return 2;
    else if (cnt == 1) return 0;
    else if (cnt == 0) return 1;
}

匈牙利算法

#include <vector>
#include <list>
#include <map>
#include <set>
#include <queue>
#include <deque>
#include <stack>
#include <bitset>
#include <algorithm>
#include <iostream>
#include <iomanip>
#include <cstdio>
#include <cmath>
#include <cstdlib>
#include <ctime>
#include <cstring>

using namespace std;
#ifdef DEBUG
#define cvar(x) cerr << "<" << #x << ": " << x << ">"
#define evar(x) cvar (x) << endl
#define debug(...) printf( __VA_ARGS__) 
template<class T> void DISP(const char *s, T x, int n) {cerr << "[" << s << ": "; for (int i = 0; i < n; ++i) cerr << x[i] << " "; cerr << "]" << endl;}
#define disp(x,n) DISP(#x " to " #n, x, n)
#else
#define debug(...) 
#define cvar(...) ({})
#define evar(...) ({})
#define disp(...) ({})
#endif

#define FILL(x) memset(x, 0, sizeof(x));
#define PB push_back

template<class T> inline T cub(T a){return a*a*a;}
template <typename T> T gcd(T x, T y) {for (T t; x; t = x, x = y % x, y = t); return y; }
typedef long long int64;


const int INF = 0x3f3f3f3f;
const int M = 10010;
const int MOD = int(1e9) + 7;
const double EPS = 1E-9;
const double  PI = acos(-1.0); //M_PI;
const int dx[] = {-1, 0, 1, 0};
const int dy[] = {0, 1, 0, -1};
int vis[M], mat[M];
vector<int> g[M];
bool find(int i);

int main()
{
#ifdef LOCAL 
    freopen("data.in","r",stdin);  
    //freopen("data.out", "w", stdout):
#endif
    int T;
    cin>>T;
    while(T--) {
        for (int i=0; i<M; i++) {
            g[i].clear();
        }
        int f, t;
        int nn;
        scanf("%d", &nn);
        FILL(vis);
        FILL(mat);
        for (int i=0; i<nn-1; i++) {
            scanf("%d%d", &f, &t);
            g[f].PB(t);
            g[t].PB(f);
        }
        int cnt=0;
        for (int i=1; i<=nn; i++) {
            FILL(vis);
            if (!mat[i])
            if (find(i)) cnt++;
        }
        if (nn % 2 == 0 && cnt == nn/2) {
            cout<<"Yes\n";
        }
        else cout<<"No\n";
    }
    return 0;
}
bool find(int f) {
    for (int i=0; i<g[f].size(); i++) {
        int t = g[f][i];
        if (!vis[t]) {
            vis[t] = 1;
            if (mat[t] == 0 || find(mat[t]))  {
                mat[t] = f;
                return 1;
            }

        }
    }
    return 0;
}                                 






                
内容概要:本文系统研究了电力系统短期负荷预测问题,提出并实现了基于极限学习机(ELM)及其智能优化改进模型的预测方法。研究涵盖标准ELM、白鲸优化算法(BWO)优化ELM和鹭鹰优化算法(IBOA)优化ELM三种模型,重点通过智能优化算法对ELM的输入权重与偏置参数进行全局寻优,有效克服了传统ELM因参数随机初始化导致的不稳定性和泛化能力不足的问题。文章完整呈现了从数据预处理、特征选择、模型构建、参数优化到预测结果对比分析的全流程,利用Matlab编程实现各模型的仿真验证,显著提升了预测精度与模型鲁棒性,为电力系统调度决策提供了可靠的技术支撑。; 适合人群:具备电力系统基础知识、时间序列预测理论及Matlab编程能力的高校研究生、科研机构研究人员以及电力公司从事负荷预测、电网调度与规划工作的技术人员。; 使用场景及目标:①应用于实际电力系统短期负荷预测业务中,提升电网运行调度的精细化与智能化水平;②作为智能优化算法与神经网络融合的经典案例,服务于学术论文撰写、科研项目申报及算法性能对比研究;③应对新能源大规模接入背景下负荷波动加剧的挑战,为构建高精度、强鲁棒性的现代负荷预测体系提供解决方案。; 阅读建议:建议读者结合所提供的Matlab代码进行动手实践,深入理解ELM网络结构与优化算法的集成机制,重点对比分析不同优化策略在收敛速度、预测误差(如MAE、RMSE、MAPE)等方面的性能差异,进而掌握智能优化技术在提升预测模型性能方面的关键作用。
内容概要:本文研究了基于Benders分解与输电网运营商(TSO)和配电网运营商(DSO)协调机制的不确定环境下输配电网双层优化模型,旨在提升高比例可再生能源接入背景下电网系统的协调性与鲁棒性。模型上层以系统整体经济性为目标进行优化调度,下层采用Benders分解实现TSO与DSO之间的信息交互与协同决策,通过引入割平面迭代机制保障求解的收敛性与全局最优性。研究充分考虑新能源出力与负荷需求的不确定性,构建了具有强适应性的双层优化框架,并基于Matlab完成了模型的编程实现与仿真验证,有效解决了多主体、多层级、多不确定性因素耦合下的电力系统优化调度难题。; 适合人群:具备电力系统分析、运筹学与优化理论基础,熟悉Matlab编程环境,从事智能电网、能源互联网、分布式能源集成、电力市场等方向的研究生、科研人员及工程技术人员。; 使用场景及目标:①研究高渗透率可再生能源条件下输配电网协同优化调度策略;②掌握Benders分解在电力系统双层优化建模中的应用方法与实现技巧;③构建TSO-DSO多主体协调机制,实现跨层级电网资源的高效互动与决策解耦;④提升对不确定性建模、分解算法设计及大规模优化问题求解能力。; 阅读建议:建议读者结合Matlab代码逐模块剖析模型构建流程,重点理解Benders割的生成逻辑、主从问题的信息传递机制及收敛判据设定,推荐在标准IEEE测试系统上复现实验以深入掌握模型特性与算法性能。
内容概要:本文提出了一种基于断线解环思想的配电网辐射状拓扑约束建模方法,旨在通过Matlab代码实现确保配电网在重构或运行过程中始终保持辐射状结构,防止环路形成,从而提升系统的安全性与稳定性。该方法通过系统性地识别网络中的潜在环路,并依据拓扑规则自动切断特定支路,有效处理配电网在优化调度、故障恢复及网络重构中的拓扑约束问题。文中详细阐述了算法的核心逻辑、数学模型构建过程、实现步骤及关键判据,并结合标准测试系统进行了仿真验证,充分证明了该方法在复杂配电网络中的有效性与实用性,尤其适用于含分布式电源接入的智能配电网场景。; 适合人群:具备一定电力系统分析基础和Matlab编程能力的高校研究生、科研人员,以及从事配电网自动化、智能电网优化、电力系统运行与控制等相关领域的工程技术人员。; 使用场景及目标:①解决配电网重构过程中的辐射状拓扑可行性验证与约束建模问题;②支撑含高比例分布式电源的配电网在故障恢复、动态重构中的安全运行分析;③为相关高水平EI期刊论文的模型复现、算法验证及科研项目申报提供可靠的代码实现与技术参考。; 阅读建议:建议读者结合Matlab代码与电力网络拓扑理论进行同步学习,重点理解断线解环的图论基础、环路搜索算法及支路断开逻辑的实现机制,并尝试在不同规模的测试系统(如IEEE 33节点系统)上进行仿真调试,以深入掌握该方法的应用技巧与优化潜力。
内容概要:本文围绕基于元模型优化算法的主从博弈多虚拟电厂动态定价与能量管理展开研究,提出了一种结合主从博弈理论与元模型优化方法的协同决策框架,通过Matlab代码实现,旨在解决高比例可再生能源接入背景下多虚拟电厂在复杂电力市场环境中的协调优化难题。研究构建了上层领导者(如主网或运营商)与下层跟随者(各虚拟电厂)之间的非对称互动模型,实现了动态电价制定与多主体能量调度的联合优化,有效提升了系统整体运行效率、经济收益与市场公平性。文中详细阐述了模型构建过程、算法设计思路及仿真验证方案,重点突出了元模型在降低计算复杂度、处理不确定性因素以及加速求解收敛方面的优势,具有较强的工程复现价值与理论参考意义。; 适合人群:具备一定电力系统运行、博弈论基础、优化建模能力及Matlab编程技能的研究生、科研人员,以及从事虚拟电厂运营、能源互联网规划、智能电网调度等相关领域的技术人员。; 使用场景及目标:①用于多主体能源系统中市场机制设计与竞价策略分析;②支撑含分布式能源的主动配电网协同优化调度研究;③为虚拟电厂参与电力市场的动态定价、需求响应与能量管理提供仿真验证平台与解决方案参考。; 阅读建议:建议读者结合Matlab代码逐模块理解算法实现流程,重点关注主从博弈架构的数学建模方式与元模型近似优化技巧的应用细节,同时可通过调整市场参数、负荷场景或可再生能源出力数据进行拓展性实验,以深化对模型鲁棒性与泛化能力的理解。
内容概要:本文围绕列车-轨道-桥梁耦合系统开展动力学交互仿真研究,基于Matlab平台构建多体动力学数值模型,综合考虑列车移动荷载、轨道结构特性与桥梁动态响应之间的耦合作用,实现对列车通过桥梁过程中振动传递规律、结构受力特性和动力响应行为的精确模拟。研究涵盖系统建模、运动方程求解、关键参数设定及仿真结果分析全过程,提供完整的Matlab代码实现方案,有助于深入理解轨道交通基础设施在运营条件下的动力性能,为桥梁结构安全性评估、轨道平顺性优化及减振设计提供理论支持和技术手段。; 适合人群:具备一定结构动力学、振动力学基础知识及Matlab编程能力的研究生、高校教师、科研机构研究人员以及从事铁路与桥梁工程设计、运维的工程技术人才。; 使用场景及目标:①用于高速铁路桥梁在列车荷载作用下的动力响应仿真与安全评估;②支撑轨道-桥梁系统减振降噪设计与结构优化;③作为高等教学与科研中的典型案例,辅助讲授多体系统动力学建模与数值仿真方法; 阅读建议:建议读者结合结构动力学相关理论教材,逐步运行并调试所提供的Matlab代码,重点关注质量-刚度-阻尼矩阵的构建、轮轨接触关系处理、时间积分算法实现等核心模块,深入理解仿真结果的物理含义及其工程应用价值。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值