[网络流24题]魔术球问题(简化版)

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
#define inf 0x7fffffff
#define T 10000
using namespace std;
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
struct edge{int to,next,v;}e[200001];
int n,cnt=1,s,ans,head[10001],h[10001];
inline void ins(int u,int v,int w){
<span style="white-space:pre">	</span>e[++cnt]=(edge){v,head[u],w};head[u]=cnt;
<span style="white-space:pre">	</span>e[++cnt]=(edge){u,head[v],0};head[v]=cnt;
}
inline bool bfs(){
<span style="white-space:pre">	</span>int t=0,w=0,now,q[20001];
<span style="white-space:pre">	</span>memset(h,-1,sizeof(h));
<span style="white-space:pre">	</span>h[0]=q[0]=0;
<span style="white-space:pre">	</span>while(t<=w){
<span style="white-space:pre">		</span>now=q[t++];
<span style="white-space:pre">		</span>for(int i=head[now];i;i=e[i].next){
<span style="white-space:pre">			</span>if(e[i].v&&h[e[i].to]==-1){
<span style="white-space:pre">				</span>h[e[i].to]=h[now]+1;
<span style="white-space:pre">				</span>q[++w]=e[i].to;
<span style="white-space:pre">			</span>}
<span style="white-space:pre">		</span>}
<span style="white-space:pre">	</span>}
<span style="white-space:pre">	</span>if(h[T]==-1)return 0;
<span style="white-space:pre">	</span>else return 1;
}
inline int dfs(int x,int f){
<span style="white-space:pre">	</span>if(x==T)return f;
<span style="white-space:pre">	</span>int used=0,rest;
<span style="white-space:pre">	</span>for(int i=head[x];i;i=e[i].next){
<span style="white-space:pre">		</span>if(e[i].v&&h[e[i].to]==h[x]+1){
<span style="white-space:pre">			</span>rest=f-used;
<span style="white-space:pre">			</span>rest=dfs(e[i].to,min(e[i].v,rest));
<span style="white-space:pre">			</span>e[i].v-=rest;
<span style="white-space:pre">			</span>e[i^1].v+=rest;
<span style="white-space:pre">			</span>used+=rest;
<span style="white-space:pre">			</span>if(used==f)return f;
<span style="white-space:pre">		</span>}
<span style="white-space:pre">	</span>}
<span style="white-space:pre">	</span>if(!used)h[x]=-1;
<span style="white-space:pre">	</span>return used;
}
int main(){
<span style="white-space:pre">	</span>freopen("balla.in","r",stdin);
<span style="white-space:pre">	</span>freopen("balla.out","w",stdout);
<span style="white-space:pre">	</span>n=read();
<span style="white-space:pre">	</span>while(1){
<span style="white-space:pre">		</span>ans++,s++;
<span style="white-space:pre">		</span>for(int i=1;i<s;i++)
<span style="white-space:pre">			</span>if(sqrt(i+s)==(int)(sqrt(i+s)))
<span style="white-space:pre">				</span>ins(i,s+5000,1);
<span style="white-space:pre">		</span>ins(0,s,1);ins(s+5000,T,1);
<span style="white-space:pre">		</span>while(bfs())ans-=dfs(0,inf);
<span style="white-space:pre">		</span>if(ans>n)break;
<span style="white-space:pre">	</span>}
<span style="white-space:pre">	</span>printf("%d\n",s-1);
<span style="white-space:pre">	</span>return 0;
}




  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值