7.29 NOI模拟赛 题答 npc问题 三染色 随机 贪心

LINK:03colors

avatar

这道题虽然绝大多数的人都获得了满分 可是我却没有.

老师讲题的时候讲到了做题答的几个技巧 这里总结一下.

数据强度大概为n=5000,m=60000的随机数据.

老师说:一部分题答是要观察数据的性质 投入时间越多得分越多.

一部分题答没有什么特点考验选手的算法的优秀程度.

不管对于前者还是后者 都可以写一个通解来先使得自己活得50~80分.

对于后者是这样的 对于前者 观察性质10几分钟 写代码十几分钟 30分钟获得10分确实不太值得.

他交了一种随机方面比较优的方法:类似模拟退火 应该算是爬山算法。

不断调整法 先随机出解 每次调整使得结果不会比原来更差。

然后有多种选择 随机选择一种。类似爬上 比爬上多了一点点随机性.

而这种调整算法的特点尽可能要保证 调整后不变劣 可以一直调整下去/尽可能多调.

实际上对于这道题 两个策略.

调整某个点 使得答案减小尽可能多 选择一样就随机选择。

当然调整不动的时候可以找到两个颜色的联通块 然后吧颜色翻转.

这样调整也是如此 不变劣且带来更多可能.

综上使用这种方法可以快速得到问题的解.

code
//#include<bits/stdc++.h>
#include<iostream>
#include<cstdio>
#include<ctime>
#include<cctype>
#include<queue>
#include<deque>
#include<stack>
#include<iostream>
#include<iomanip>
#include<cstdio>
#include<cstring>
#include<string>
#include<ctime>
#include<cmath>
#include<cctype>
#include<cstdlib>
#include<queue>
#include<deque>
#include<stack>
#include<vector>
#include<algorithm>
#include<utility>
#include<bitset>
#include<set>
#include<map>
#define ll long long
#define db double
#define INF 1000000000
#define inf 100000000000000000ll
#define ldb long double
#define pb push_back
#define put_(x) printf("%d ",x);
#define get(x) x=read()
#define gt(x) scanf("%d",&x)
#define gi(x) scanf("%lf",&x)
#define put(x) printf("%d
",x)
#define putl(x) printf("%lld
",x)
#define rep(p,n,i) for(RE int i=p;i<=n;++i)
#define go(x) for(int i=lin[x],tn=ver[i];i;tn=ver[i=nex[i]])
#define fep(n,p,i) for(RE int i=n;i>=p;--i)
#define vep(p,n,i) for(RE int i=p;i<n;++i)
#define pii pair<int,int>
#define mk make_pair
#define RE register
#define P 1000000007ll
#define gf(x) scanf("%lf",&x)
#define pf(x) ((x)*(x))
#define uint unsigned long long
#define ui unsigned
#define EPS 1e-10
#define sq sqrt
#define S second
#define F first
#define mod 998244353
#define id(i,j) ((i-1)*m+j)
#define max(x,y) ((x)<(y)?y:x)
using namespace std;
char *fs,*ft,buf[1<<15];
inline char gc()
{
	return (fs==ft&&(ft=(fs=buf)+fread(buf,1,1<<15,stdin),fs==ft))?0:*fs++;
}
inline int read()
{
	RE int x=0,f=1;RE char ch=gc();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=gc();}
	return x*f;
}
const int MAXN=200010,G=3;
int n,m,len;
int c[MAXN];
int lin[MAXN],ver[MAXN<<1],nex[MAXN<<1];
inline int check()
{
	rep(1,n,j)
	{
		go(j)
		{
			if(c[tn]==c[j])return 0;
		}
	}
	return 1;
}
inline void add(int x,int y)
{
	ver[++len]=y;nex[len]=lin[x];lin[x]=len;
	ver[++len]=x;nex[len]=lin[y];lin[y]=len;
}
	int b[4];
int main()
{
	freopen("1.in","r",stdin);
	freopen("1.out","w",stdout);
	srand(time(0));
	get(n);get(m);
	rep(1,m,i)add(read(),read());
	rep(1,n,i)c[i]=rand()%3+1;
	//put(n);	

	while(!check())
	{
		int mx=0,p,cc;
		rep(1,n,j)
		{
			int ww=0;
			b[1]=b[2]=b[3]=0;
			go(j)++b[c[tn]];
			ww=b[c[j]];
			rep(1,3,k)
			{
				int vv=0;
				if(c[j]==k)continue;
				vv=ww-b[k];
				if(vv==mx&&rand()%2==0)p=j,cc=k;
				if(vv>mx)mx=vv,p=j,cc=k;
			}
		}
		c[p]=cc;
		//if(!mx)cout<<mx<<' '<<check()<<' '<<p<<' '<<cc<<endl;
		//cout<<mx<<endl;
	}
	//rep(1,n,i)get(c[i]);
	rep(1,n,i)cout<<c[i]<<' ';
	//cout<<check();
	return 0;
}
原文地址:https://www.cnblogs.com/chdy/p/13401764.html