分析:做这个题目的时候是在图论的题目里面看到的,到时读了题目推了一下,发现好像有点规律,试了一下果然过了。
后来看了一下数据,才50,那么试了一下模拟,也过了。
好像zoj有一道题目卡模拟,模拟的时候必须贪心一下才能过
这道题出题人的意图在于考大家的:二分图最小路径覆盖。
把每一个球看做一个点,然后如果两个和为平方数的话就给这两个点之间连接一条边,然后用一个特殊的匹配算法,类似于匈牙利算法,但是每次找匹配的时候加入一条边上连接的有匹配的话就不能匹配,最后求一个最大匹配就好了。这个题目50个点,当然要预处理成离线的。
如果每次都建图跑的话时间复杂度是非常高的,所以我们可以用标记的方法,如果匹配过的就不用再匹配了。这样能降低很多时间,但是还是很长。
找规律代码:
#include <cstdio>
int ans[55];
void isit()
{
ans[1]=1,ans[2]=3;
for(int i=3;i<51;i++)
ans[i]=ans[i-1]+(i+1)/2 * 2;
}
int main()
{
int T,n;isit();
scanf("%d",&T);
for(;T--&&scanf("%d",&n);printf("%d\n",ans[n]));
}
#include <cstdio>
#include <stack>
#include <cmath>
using namespace std;
stack<int> st[60];
int ans[60];
void isit()
{
int tmp=1,i=1;
while(i<55)
{
int ok=0;
for(i=1;!st[i].empty();i++)
{
int zhi = st[i].top()+tmp;
int sq=sqrt(zhi);
if(sq * sq == zhi)
{
st[i].push(tmp);
ok=1;
break;
}
}
if(!ok){
st[i].push(tmp);
ans[i]=tmp-1;
}
tmp++;
}
}
int main()
{
int T,n;isit();
scanf("%d",&T);
for(;T--&&scanf("%d",&n);printf("%d\n",ans[n+1])){}
}
二分图匹配代码:
#include <cstdio>
#include <stack>
#include <cmath>
#include <cstring>
using namespace std;
const int N = 60;
const int M = 1500;
int ans[N];
bool ok[M*2+10];
bool mp[M][M];
int mx[M],my[M]; //mx保存正向边 my保存反向边
bool used[M]; //
bool dfs(int v) {
for (int u = 1; u < v; ++u) //优化
if (mp[v][u] && !used[u]) {
used[u] = true;
if (my[u] == -1 || dfs(my[u])) { //一直向下找
my[u] = v;
mx[v] = u;
return true;
}
}
return false;
}
int Max_Pi(int ans,int n) //最大匹配
{
for(int i=1;i<=n;i++)
{
if(mx[i]<0){
memset(used,false,sizeof(used));
if(dfs(i))
ans++;
}
}
return ans;
}
void Min_Road() //最小路径覆盖
{
memset(mx,-1,sizeof(mx)); //优化
memset(my,-1,sizeof(my));
for(int i=1;i*i<=2*M;i++)
ok[i*i]=true;
ans[1]=1;
int tmp=0;
for(int i=2;i<=M;i++)
{
for(int j=1;j<i;j++)
if(ok[i+j])
mp[i][j]=true;
tmp=Max_Pi(tmp,i);
//printf("xx%d %d",i,tmp);
ans[i-tmp]=i;
if(i-tmp>55)
break;
}
// for(int i=1;i<=50;i++)
// printf("%d ",ans[i]);
}
int main()
{
int T,n;Min_Road();
scanf("%d",&T);
for(;T--&&scanf("%d",&n);printf("%d\n",ans[n])){}
}