没贴网络流标签绝对想不到网络流orz。。或许会想到贪心去莽一发?
先说说网络流的。。
其实很容易想到把每个数字作为一个点,对能构成平方数的数对连边,这样一来一根柱子就能看作一个路径,然后就变成了最小覆盖路径了。。。
所以对给定答案是能判断其最大柱子数,所以很容易想到用二分去做。。。
事实上好像用枚举会更快?因为重建图重跑一遍最大流的时间比相比于枚举实在是要大得多了。。不过这种感觉蛮难写的。。重建图思路简单粗暴。。
然后这题竟然还能贪心。。
思路是从前到后寻找柱子,找不到再重开一个。。听起来就十分不靠谱。。
窝大概是这么理解的。。前面的柱子肯定要比后面的柱子长,顶端的数字肯定要比后面的数字要大。。而数字越大就 越难找到匹配的数字构成平方和,因此如果和大数能够构成平方和应该是比较优的方案了。。
然而严格证明并不会。。
然后有人甚至总结出了公式= =!不过不会证明就不贴上来了。。
网络流:
/**
* ┏┓ ┏┓
* ┏┛┗━━━━━━━┛┗━━━┓
* ┃ ┃
* ┃ ━ ┃
* ┃ > < ┃
* ┃ ┃
* ┃... ⌒ ... ┃
* ┃ ┃
* ┗━┓ ┏━┛
* ┃ ┃ Code is far away from bug with the animal protecting
* ┃ ┃ 神兽保佑,代码无bug
* ┃ ┃
* ┃ ┃
* ┃ ┃
* ┃ ┃
* ┃ ┗━━━┓
* ┃ ┣┓
* ┃ ┏┛
* ┗┓┓┏━┳┓┏┛
* ┃┫┫ ┃┫┫
* ┗┻┛ ┗┻┛
*/
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<queue>
#include<cmath>
#include<map>
#include<stack>
#include<set>
#define inc(i,l,r) for(int i=l;i<=r;i++)
#define dec(i,l,r) for(int i=l;i>=r;i--)
#define link(x) for(edge *j=h[x];j;j=j->next)
#define mem(a) memset(a,0,sizeof(a))
#define ll long long
#define eps 1e-12
#define succ(x) (1<<x)
#define lowbit(x) (x&(-x))
#define sqr(x) ((x)*(x))
#define mid (x+y>>1)
#define NM 20005
#define nm 100005
#define pi 3.1415926535897931
using namespace std;
const int inf=1000000005;
ll read(){
ll x=0,f=1;char ch=getchar();
while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
return f*x;
}
struct edge{int t,v;edge*next,*rev;}e[nm],*h[NM],*o=e,*tmp[NM],*p[NM];
void _add(int x,int y,int v){o->t=y;o->v=v;o->next=h[x];h[x]=o++;}
void add(int x,int y,int v){_add(x,y,v);_add(y,x,0);h[x]->rev=h[y];h[y]->rev=h[x];}
int n,m,tot,d[NM],cnt[NM],ans;
bool v[NM];
int maxflow(){
int flow=0;edge*j;
inc(i,0,n)tmp[i]=h[i];
cnt[0]=tot=n+1;
for(int x=0,s=inf;d[x]<tot;){
for(j=tmp[x];j;j=j->next)if(j->v&&d[j->t]+1==d[x])break;
if(j){
s=min(s,j->v);p[j->t]=tmp[x]=j;
if((x=j->t)==n){
for(;p[x];x=p[x]->rev->t)p[x]->v-=s,p[x]->rev->v+=s;
flow+=s;s=inf;
}
}else{
if(!--cnt[d[x]])break;d[x]=tot;
link(x)if(j->v&&d[j->t]+1<d[x])d[x]=d[j->t]+1,tmp[x]=j;
cnt[d[x]]++;
if(p[x])x=p[x]->rev->t;
}
}
return flow;
}
bool check(int t){
n=t*2+1;
mem(e);mem(h);o=e;mem(d);mem(cnt);mem(tmp);mem(p);
inc(i,1,t)inc(j,i+1,t)if(sqr((int)sqrt(i+j))==i+j)add(i,j+t,1);
inc(i,1,t)add(0,i,1),add(i+t,n,1);
int s=maxflow();
// printf("%d %d\n",t,s);
return t-s<=m;
}
int main(){
m=read();
for(int x=1,y=2000;x<=y;)if(check(mid)){ans=mid;x=mid+1;}else y=mid-1;
printf("%d\n",ans);
check(ans);
inc(i,1,ans)if(!v[i]){
bool f=true;
for(int t=i;f;v[t]++){
printf("%d ",t);f=false;
link(t)if(j->t>ans&&!j->v){t=j->t-ans;f=true;break;}
}
putchar('\n');
}
return 0;
}
贪心:
/**
* ┏┓ ┏┓
* ┏┛┗━━━━━━━┛┗━━━┓
* ┃ ┃
* ┃ ━ ┃
* ┃ > < ┃
* ┃ ┃
* ┃... ⌒ ... ┃
* ┃ ┃
* ┗━┓ ┏━┛
* ┃ ┃ Code is far away from bug with the animal protecting
* ┃ ┃ 神兽保佑,代码无bug
* ┃ ┃
* ┃ ┃
* ┃ ┃
* ┃ ┃
* ┃ ┗━━━┓
* ┃ ┣┓
* ┃ ┏┛
* ┗┓┓┏━┳┓┏┛
* ┃┫┫ ┃┫┫
* ┗┻┛ ┗┻┛
*/
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<queue>
#include<cmath>
#include<map>
#include<stack>
#include<set>
#define inc(i,l,r) for(int i=l;i<=r;i++)
#define dec(i,l,r) for(int i=l;i>=r;i--)
#define link(x) for(edge *j=h[x];j;j=j->next)
#define mem(a) memset(a,0,sizeof(a))
#define ll long long
#define eps 1e-12
#define succ(x) (1<<x)
#define lowbit(x) (x&(-x))
#define sqr(x) ((x)*(x))
#define mid (x+y>>1)
#define NM 20005
#define nm 100005
#define pi 3.1415926535897931
using namespace std;
const int inf=1000000005;
ll read(){
ll x=0,f=1;char ch=getchar();
while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
return f*x;
}
int n,cnt;
vector<int>a[NM];
bool check(int t){return sqr((int)sqrt(t))==t;}
int main(){
n=read();
inc(k,1,1e6){
bool f=false;
inc(i,1,cnt)
if(check(k+a[i][a[i].size()-1])){a[i].push_back(k);f++;break;}
if(!f)a[++cnt].push_back(k);
if(cnt>n){printf("%d\n",k-1);break;}
}
inc(i,1,n){
int t=a[i].size()-1;
inc(j,0,t)printf("%d ",a[i][j]);
putchar('\n');
}
}
P2765 魔术球问题
题目描述
«问题描述:
假设有n根柱子,现要按下述规则在这n根柱子中依次放入编号为1,2,3,...的球。
(1)每次只能在某根柱子的最上面放球。
(2)在同一根柱子中,任何2个相邻球的编号之和为完全平方数。
试设计一个算法,计算出在n根柱子上最多能放多少个球。例如,在4 根柱子上最多可放11 个球。
«编程任务:
对于给定的n,计算在n根柱子上最多能放多少个球。
输入输出格式
输入格式:第1 行有1个正整数n,表示柱子数。
输出格式:程序运行结束时,将n 根柱子上最多能放的球数以及相应的放置方案输出。文件的第一行是球数。接下来的n行,每行是一根柱子上的球的编号。
输入输出样例
说明
感谢 @PhoenixEclipse 提供spj
4<=n<=55