1006: [HNOI2008]神奇的国度

9 篇文章 0 订阅
0 篇文章 0 订阅

1006: [HNOI2008]神奇的国度

Time Limit: 20 Sec   Memory Limit: 162 MB
Submit: 3494   Solved: 1592
[ Submit][ Status][ Discuss]

Description

  K国是一个热衷三角形的国度,连人的交往也只喜欢三角原则.他们认为三角关系:即AB相互认识,BC相互认识,CA
相互认识,是简洁高效的.为了巩固三角关系,K国禁止四边关系,五边关系等等的存在.所谓N边关系,是指N个人 A1A2
...An之间仅存在N对认识关系:(A1A2)(A2A3)...(AnA1),而没有其它认识关系.比如四边关系指ABCD四个人 AB,BC,C
D,DA相互认识,而AC,BD不认识.全民比赛时,为了防止做弊,规定任意一对相互认识的人不得在一队,国王相知道,
最少可以分多少支队。

Input

  第一行两个整数N,M。1<=N<=10000,1<=M<=1000000.表示有N个人,M对认识关系. 接下来M行每行输入一对朋

Output

  输出一个整数,最少可以分多少队

Sample Input

4 5
1 2
1 4
2 4
2 3
3 4

Sample Output

3

HINT

  一种方案(1,3)(2)(4)

Source

[ Submit][ Status][ Discuss]

基于图中的一个环,我们定义以下概念:
弦:连接一个简单环中不相邻的两个点的边。
弦图:任意长度大于3的简单环都至少包含一条弦的图。
单纯点:若点x和它的邻域组成的点集的生成子图为一个团,则称x是单纯点。
显然,基于题设的限制,给出的图一定是一张弦图
原题等价于,给每个点染色,要求相邻的点不同色
对这张弦图,求出一个完美消除序列,即一个点集序列,满足:
对于序列中任意一个后缀点集的生成子图该后缀的第一个点是单纯点
有个结论,在完美消除序列里,从后往前染色,每个点取当前能染的编号最小的颜色,这样就是最优了
然后完美消除序列可用最大势算法求解,从n~1给每个点标号
记label[i]:与i相邻的已经标号的点的数量
每次取label值最大的点给出当前标号
最后,完美消除序列的第k个点就是标号为k的点
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<vector>
#include<queue>
#include<set>
#include<map>
#include<stack>
#include<bitset>
#include<ext/pb_ds/priority_queue.hpp>
using namespace std;
 
const int maxn = 1E4 + 10;
 
struct data{
    int Num,w; data(){}
    data(int Num,int w): Num(Num),w(w){}
    bool operator < (const data &B) const {return w < B.w;}
};
typedef __gnu_pbds::priority_queue<data,less<data>,__gnu_pbds::pairing_heap_tag> Heap;
 
int n,m,cnt,label[maxn],Num[maxn],c[maxn],bo[maxn];
bool vis[maxn];
 
Heap Q;
vector <int> v[maxn];
Heap::point_iterator id[maxn];
 
int getint()
{
    char ch = getchar(); int ret = 0;
    while (ch < '0' || '9' < ch) ch = getchar();
    while ('0' <= ch && ch <= '9')
        ret = ret * 10 + ch - '0',ch = getchar();
    return ret;
}
 
int main()
{
    #ifdef DMC
        freopen("DMC.txt","r",stdin);
    #endif
     
    n = getint(); m = getint();
    while (m--)
    {
        int x = getint(),y = getint();
        v[x].push_back(y);
        v[y].push_back(x);
    }
    for (int i = 1; i <= n; i++) id[i] = Q.push(data(i,0));
    for (int i = n; i; i--)
    {
        data K = Q.top(); Q.pop();
        int k = K.Num; Num[i] = k; vis[k] = 1;
        for (int j = 0; j < v[k].size(); j++)
        {
            int to = v[k][j];
            if (vis[to]) continue; ++label[to];
            Q.modify(id[to],data(to,label[to]));
        }
    }
    for (int i = n; i; i--)
    {
        int k = Num[i]; ++cnt;
        for (int j = 0; j < v[k].size(); j++)
            bo[c[v[k][j]]] = cnt;
        for (int j = 1; ; j++)
            if (bo[j] != cnt) {c[k] = j; break;}
    }
    int Ans = 0;
    for (int i = 1; i <= n; i++) Ans = max(Ans,c[i]);
    cout << Ans << endl;
    return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
在信号处理领域,DOA(Direction of Arrival)估计是一项关键技术,主要用于确定多个信号源到达接收阵列的方向。本文将详细探讨三种ESPRIT(Estimation of Signal Parameters via Rotational Invariance Techniques)算法在DOA估计中的实现,以及它们在MATLAB环境中的具体应用。 ESPRIT算法是由Paul Kailath等人于1986年提出的,其核心思想是利用阵列数据的旋转不变性来估计信号源的角度。这种算法相比传统的 MUSIC(Multiple Signal Classification)算法具有较低的计算复杂度,且无需进行特征值分解,因此在实际应用中颇具优势。 1. 普通ESPRIT算法 普通ESPRIT算法分为两个主要步骤:构造等效旋转不变系统和估计角度。通过空间平移(如延时)构建两个子阵列,使得它们之间的关系具有旋转不变性。然后,通过对子阵列数据进行最小二乘拟合,可以得到信号源的角频率估计,进一步转换为DOA估计。 2. 常规ESPRIT算法实现 在描述中提到的`common_esprit_method1.m`和`common_esprit_method2.m`是两种不同的普通ESPRIT算法实现。它们可能在实现细节上略有差异,比如选择子阵列的方式、参数估计的策略等。MATLAB代码通常会包含预处理步骤(如数据归一化)、子阵列构造、旋转不变性矩阵的建立、最小二乘估计等部分。通过运行这两个文件,可以比较它们在估计精度和计算效率上的异同。 3. TLS_ESPRIT算法 TLS(Total Least Squares)ESPRIT是对普通ESPRIT的优化,它考虑了数据噪声的影响,提高了估计的稳健性。在TLS_ESPRIT算法中,不假设数据噪声是高斯白噪声,而是采用总最小二乘准则来拟合数据。这使得算法在噪声环境下表现更优。`TLS_esprit.m`文件应该包含了TLS_ESPRIT算法的完整实现,包括TLS估计的步骤和旋转不变性矩阵的改进处理。 在实际应用中,选择合适的ESPRIT变体取决于系统条件,例如噪声水平、信号质量以及计算资源。通过MATLAB实现,研究者和工程师可以方便地比较不同算法的效果,并根据需要进行调整和优化。同时,这些代码也为教学和学习DOA估计提供了一个直观的平台,有助于深入理解ESPRIT算法的工作原理。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值