50 years, 50 colors

这个题是最小点覆盖

根据定理也就是求二分图的最大匹配(证明没看懂)

就是我们把行拿出来作为左部分,列拿出来作为右部分,然后两者相连的边就代表了原图中的点

#include<iostream>
#include<cstdio>
#include<map>
#include<string.h>
using namespace std;
int un,uv;
map<int,int>m;
const int maxn = 100+10;
int g[maxn][maxn];
int linker[maxn];
int used[maxn];
int n,k;
int res[100];
bool dfs(int u,int t){
    for(int v = 0;v<uv;++v){
        if(g[u][v]==t&&!used[v]){
            used[v] = true;
            if(linker[v]==-1||dfs(linker[v],t)){
            linker[v] = u;
            return true;}
        }
    }
    return false;
}

int hungry(int t){
    int res = 0;
    memset(linker,-1,sizeof(linker));
    for(int u =0;u<un;++u){
        memset(used,false,sizeof(used));
        if(dfs(u,t))res++;
    }
    return res;
}
int main(){
    while(~scanf("%d%d",&n,&k)&&n+k!=0){
            m.clear();
            un = n,uv= n;
        for(int i = 0;i<n;++i){
            for(int j =0;j<n;++j){
                scanf("%d",&g[i][j]);
                m[g[i][j]]=1;
            }
        }

        int h = 0;
        map<int,int>::iterator it;
        for(it = m.begin();it!=m.end();++it){
            if(hungry(it->first)>k){
                res[h++]=it->first;
            }
        }
        if(h==0)
            cout<<-1<<endl;
        else{
                for(int i =0;i<h;++i){
                    if(i)
                        cout<<" "<<res[i];
                    else
                        cout<<res[i];
                }
                cout<<endl;
            }

    }

}

 

以下是可能的 Matlab 代码实现: ``` % 读取 NDBI 的 CSV 文件 ndbi_data = csvread('ndbi.csv'); % 获取年份数据和 NDBI 数据 years = ndbi_data(:, 1); ndbi_values = ndbi_data(:, 2:end); % 使用似然比检验方法分割时间序列 [num_years, num_locations] = size(ndbi_values); threshold = 0.05; % 设置显著性水平为 0.05 colors = {'r', 'g', 'b', 'm', 'k'}; % 不同的颜色 scatter_data = cell(1, length(colors)); num_clusters = 0; for loc_idx = 1:num_locations ndbi_loc = ndbi_values(:, loc_idx); [change_points, ~] = findchangepts(ndbi_loc, 'MaxNumChanges', num_years-2); if isempty(change_points) continue; end p_values = zeros(1, length(change_points)); for cp_idx = 1:length(change_points) cp = change_points(cp_idx); [h, p_values(cp_idx)] = vartest2(ndbi_loc(1:cp), ndbi_loc(cp+1:end), 'Alpha', threshold); end significant_changes = find(p_values < threshold); if ~isempty(significant_changes) num_clusters = num_clusters + 1; scatter_data{num_clusters} = [years, ndbi_loc]; for sc_idx = 1:length(significant_changes) change_idx = significant_changes(sc_idx); change_year = years(change_points(change_idx)); plot(change_year, ndbi_loc(change_points(change_idx)), 'ko', 'MarkerSize', 10, 'LineWidth', 2); end cluster_color = colors{mod(num_clusters-1, length(colors))+1}; plot(years, ndbi_loc, '.', 'Color', cluster_color, 'MarkerSize', 20); end end % 输出分割后散点图 figure; hold on; for cluster_idx = 1:num_clusters cluster_data = scatter_data{cluster_idx}; cluster_color = colors{mod(cluster_idx-1, length(colors))+1}; scatter(cluster_data(:, 1), cluster_data(:, 2), 50, cluster_color, 'filled'); end xlabel('Year'); ylabel('NDBI value'); ``` 请注意,此代码仅提供了一种可能的实现方法。具体实现可能因数据格式、分析目标等因素而异。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值