LiberOJ 6003. 「网络流 24 题」魔术球贪心或者最小路径覆盖

Posted 2020-10-02 GeekZRF

tags:

篇首语：本文由小常识网(cha138.com)小编为大家整理，主要介绍了LiberOJ 6003. 「网络流 24 题」魔术球贪心或者最小路径覆盖相关的知识，希望对你有一定的参考价值。

内存限制：256 MiB时间限制：1000 ms标准输入输出

题目类型：传统评测方式：Special Judge

上传者：匿名

提交提交记录统计讨论测试数据

假设有 $\cdots1,2,3,4,? 的球。$

每次只能在某根柱子的最上面放球。
在同一根柱子中，任何

试设计一个算法，计算出在

文件第

第一行是球数。接下来的

样例输入

样例输出

$\leq n \leq 551≤n≤55$

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<set>
#include<queue>
#include<stack>
#include<map>
#include<vector>
using namespace std;
typedef long long ll;
typedef pair<int,int> P;
const int maxn=1e6+100,maxm=1e5+1000,inf=0x3f3f3f3f,mod=1e9+7;
const ll INF=1e18+7;
vector<ll>G[maxn];
bool used[maxn];
int cx[maxn],cy[maxn];
bool vis[maxn];
int res=0;
bool dfs(int u)
{
    for(int i=0; i<G[u].size(); i++)
    {
        int v=G[u][i];
        if(used[v]) continue;
        used[v]=true;
        if(cy[v]<0||dfs(cy[v]))
        {
            cx[u]=v,cy[v]=u;
            return true;
        }
    }
    return false;
}
int solve(int n)
{
    res+=dfs(n);
    return res;
}
int main()
{
    int n;
    scanf("%d",&n);
    int ans=1;
    memset(cx,-1,sizeof(cx));
    memset(cy,-1,sizeof(cy));
    for(int i=1; i<=2000000; i++)
    {
        int x=(int)(sqrt(i+1));
        while(x*x-i<i)
        {
            if(x*x-i>=1) G[i].push_back(x*x-i);
            x++;
        }
        if(i-solve(i)>n) break;
        ans=i;
    }
    memset(vis,false,sizeof(vis));
    printf("%d\n",ans);
    for(int i=ans; i>=1; i--)
    {
        if(vis[i]) continue;
        printf("%d",i);
        for(int j=i; cx[j]!=-1; j=cx[j],vis[j]=true)
            printf(" %d",cx[j]);
        printf("\n");
    }
}