6.14 YZBOI模拟赛solution
\(6.14\ YZBOI\)模拟赛\(solution\)
本来不想写题解来着。。。毕竟是自己找的题还是写一写吧
上午为了整活,就把赛制改成\(IOI\)赛制了,于是乎拯救了大家的\(70pts\)代码。
三道巨水的题。。。\(bz\)迅速就\(AK\)了。。。
\(T1\ Delta\)
比较直观的想,肯定是找一个少的的完全图划分了,又发现平方具有传递性,就变成维护联通块了
#include
#define int long long
#define MAXN 5005
using namespace std;
int n,a[MAXN],fa[MAXN],ans[MAXN];
bool vis[MAXN];
int Find(int x)
{
if(fa[x]==x) return x;
return fa[x]=Find(fa[x]);
}
signed main()
{
scanf("%lld",&n);
for(int i=1;i<=n;i++)
{
fa[i]=i;
scanf("%lld",&a[i]);
}
for(int i=1;i<=n;i++)
{
for(int j=i+1;j<=n;j++)
{
if(a[i]*a[j]>0)
{
if((int)sqrt(a[i]*a[j])*(int)sqrt(a[i]*a[j])!=a[i]*a[j]) continue;
int fi=Find(i);
int fj=Find(j);
if(fi==fj) continue;
fa[fi]=fj;
}
}
}
for(int i=1;i<=n;i++)
{
memset(vis,0,sizeof(vis));
int num=0;
for(int j=i;j<=n;j++)
{
if(a[j]!=0&&vis[Find(j)]==false)
{
vis[Find(j)]=true;
num++;
}
if(num==0) ans[num+1]++;
else ans[num]++;
}
}
for(int i=1;i<=n;i++)
{
cout<
\(T2 Add\)
线性规划会吧,模拟会吧,单纯形法会吧,转对偶会吧,那么这就是一道水题了
#include
using namespace std;
#define MAXN 100005
const double eps=1e-9;
const int INF=1000000000;
double a[1005][MAXN],b[MAXN],v,c[MAXN];
int m,n;
void solve(int l,int e)
{
b[l]/=a[l][e];
for(int i=1;i<=n;++i) if(i!=e) a[l][i]/=a[l][e];
a[l][e]=1/a[l][e];
for(int i=1;i<=m;++i)
{
if(i!=l&&abs(a[i][e])>eps)
{
b[i]-=b[l]*a[i][e];
for(int j=1;j<=n;++j)
{
if(j!=e) a[i][j]-=a[l][j]*a[i][e];
}
a[i][e]=-a[i][e]*a[l][e];
}
}
v+=c[e]*b[l];
for(int i=1;i<=n;++i)
{
if(i!=e) c[i]-=c[e]*a[l][i];
}
c[e]=-c[e]*a[l][e];
}
void simplex()
{
int l,e;
while(true)
{
for(e=1;e<=n;++e)
{
if(c[e]>eps) break;
}
if(e>n) break;
double t=INF;
for(int i=1;i<=m;++i)
{
if(a[i][e]>eps&&b[i]/a[i][e]>m>>n;
for(int i=1;i<=m;++i) scanf("%lf",&b[i]);
for(int i=1,l,r;i<=n;++i)
{
scanf("%d%d",&l,&r);
for(int j=l;j<=r;++j) ++a[j][i];
scanf("%lf",&c[i]);
}
simplex();
printf("%.0lf\n",v);
return 0;
}
\(T3\ Times\)
我记得这道题是我省选前一天为了帮忙\(debug\)写的
考虑最优解必然是在最小生成树上的,然后考虑从大往小删除边,我们需要保证删边之后每个联通块至少有一个守卫,那么大概就是我们对每个守卫匹配一个边去割掉,而且不存在交叉,那么就可以很直观的用费用流来表示
#include
#define MIN 0xcfcfcfcfcfcfcfcf
#define INF 2147483647
#define int long long
#define MAXN 1000005
using namespace std;
struct node
{
int u,v,w;
}lu[MAXN];
int head[MAXN],cost[MAXN],val[MAXN],nxt[MAXN],to[MAXN],tot=1;
int dis[MAXN],fa[MAXN],Sum,s,t,n,m,k;
bool vis[MAXN],in[MAXN];
vector >rd[MAXN];
int Find(int x)
{
if(fa[x]==x) return x;
return fa[x]=Find(fa[x]);
}
int R(int x)
{
return x+n+n;
}
int ru(int x)
{
return x;
}
int chu(int x)
{
return x+n;
}
void add(int u,int v,int w,int cts)
{
tot++; to[tot]=v; val[tot]=w; cost[tot]=cts; nxt[tot]=head[u]; head[u]=tot;
swap(u,v),w=0,cts*=-1;
tot++; to[tot]=v; val[tot]=w; cost[tot]=cts; nxt[tot]=head[u]; head[u]=tot;
}
bool bfs(int s,int t)
{
memset(dis,0xcf,sizeof(dis));
queueq;
dis[s]=0;
in[s]=true;
q.push(s);
while(q.size())
{
int now=q.front();
q.pop();
in[now]=false;
for(int i=head[now];i;i=nxt[i])
{
int y=to[i];
if(val[i]&&dis[y]=INF)
{
cout<<-1;
return 0;
}
cout<