其他分享
首页 > 其他分享> > noip模拟72(待补)

noip模拟72(待补)

作者:互联网

T1 出了个大阴间题

状压 \(dp\),设 \(f_{s,k}\) 表示 \(s\) 状态下,当前的 \(max_a -\) 最大值等于 \(j\) 的所有代价.

有一个不会对分数有什么影响的规律,就是 \(k\) 永远 \(\le 1\).

A_code
#include<bits/stdc++.h>
using namespace std;
namespace BSS {
	#define ll long long 
	#define ull unsigned ll
	#define lf long double
	#define lbt(x) (x&(-x))
	#define mp(x,y) make_pair(x,y)
	#define lb lower_bound 
	#define ub upper_bound
	#define Fill(x,y) memset(x,y,sizeof x)
	#define Copy(x,y) memcpy(x,y,sizeof x)
	#define File(x) freopen(#x".in","r",stdin),freopen(#x".out","w",stdout)
	inline ll read() {
		ll w=0; bool cit=1; char ch;
		while(!isdigit(ch=getchar())) if(ch=='-') cit=0; 
		while(isdigit(ch)) w=(w<<3)+(w<<1)+(ch^48),ch=getchar();
		return cit?w:-w;
	}
} using namespace BSS;

const ll mod=1e9+7;

ll m,n,U,maxval;
ll val[25],b[25];
ll num[1<<18];
map<ll,ll> f[1<<18],g[1<<18];
signed main(){
	File(repair);
	n=read(),m=read()%mod,U=(1<<n)-1;
	for(ll i=1;i<=n;i++) val[i]=read();
	for(ll i=2;i<=n;i++) b[i]=b[i-1]<<1|1;
	for(ll i=1;i<=n;i++) f[1<<i-1][val[i]]=0,g[1<<i-1][val[i]]++;
	for(ll i=1;i<U;i++){
		num[i]=num[i&(i-1)]+1;
		for(ll j=1,k;j<=n;j++){
			if((i>>j-1)&1) continue; k=i|(1<<j-1);
			for(auto c : f[i]){
				ll to=(c.first==val[j] ? val[j]+1 : max(c.first,val[j]));
				f[k][to]=(f[k][to]+f[i][c.first]+(to*m%mod+b[num[i]])%mod*g[i][c.first]%mod)%mod;
				g[k][to]=(g[k][to]+g[i][c.first])%mod;
			}
		}
	}
	for(auto i : f[U]) maxval=max(maxval,i.first);
	printf("%lld %lld\n",maxval,f[U][maxval]%mod),exit(0);
}

T2 最简单辣快来做

自以为是打的是正解,其实连别人的暴力都不如.

自以为把简单的想复杂就可能是正解,像个笑话.

正解是利用分块思想,考虑把每个横纵坐标都有当成一条线,于是变成了 \(n^2\) 个矩形.

然后前缀和即可.

学习了光速幂.

B_code
#include<bits/stdc++.h>
using namespace std;
namespace BSS {
	#define ll long long 
	#define ull unsigned ll
	#define lf long double
	#define lbt(x) (x&(-x))
	#define mp(x,y) make_pair(x,y)
	#define lb lower_bound 
	#define ub upper_bound
	#define Fill(x,y) memset(x,y,sizeof x)
	#define Copy(x,y) memcpy(x,y,sizeof x)
	#define File(x) freopen(#x".in","r",stdin),freopen(#x".out","w",stdout)
	inline ll read() {
		ll w=0; bool cit=1; char ch;
		while(!isdigit(ch=getchar())) if(ch=='-') cit=0; 
		while(isdigit(ch)) w=(w<<3)+(w<<1)+(ch^48),ch=getchar();
		return cit?w:-w;
	}
} using namespace BSS;

const ll N=2e3+21,W=1e5;

ll m,n,s,ops,mod,cntx,cnty;
ll lshx[N],lshy[N];
ll rk[N][N];
ll pre[5][N][N];
struct I { ll x,y,h,id; } p[N*100];
struct Flash_Pow{
	ll val; ll sm[W+21],bg[W+21],po[N];
	inline void Pre(){
		sm[0]=1; for(ll i=1;i<=W;i++) sm[i]=sm[i-1]*val%mod;
		bg[0]=1; for(ll i=1;i<=W;i++) bg[i]=bg[i-1]*sm[W]%mod;
	}
	inline ll qpow(ll x){ return sm[x%W]*bg[x/W]%mod; }
} A,B;
signed main(){
	File(satellite);
	s=read(),ops=read(),n=read(),m=read(); ll x,y,res;
	mod=read(),A.val=read(),B.val=read(),A.Pre(),B.Pre();
	for(ll i=1;i<=s;i++){
		p[i].h=read()%mod,p[i].id=i;
		lshx[++cntx]=(p[i].x=read()),lshy[++cnty]=(p[i].y=read());
	}
	sort(lshx+1,lshx+1+cntx),cntx=unique(lshx+1,lshx+1+cntx)-lshx-1;
	sort(lshy+1,lshy+1+cnty),cnty=unique(lshy+1,lshy+1+cnty)-lshy-1;
	for(ll i=1;i<=cntx;i++) A.po[i]=A.qpow(lshx[i]-lshx[i-1]);
	for(ll i=1;i<=cnty;i++) B.po[i]=B.qpow(lshy[i]-lshy[i-1]);
	for(ll i=1;i<=s;i++){
		p[i].x=lb(lshx+1,lshx+1+cntx,p[i].x)-lshx,p[i].y=lb(lshy+1,lshy+1+cnty,p[i].y)-lshy;
		rk[p[i].x][p[i].y]=p[i].id;
		pre[1][p[i].x][p[i].y]+=p[i].h,pre[2][p[i].x][p[i].y]+=p[i].h;
		pre[3][p[i].x][p[i].y]+=p[i].h,pre[4][p[i].x][p[i].y]+=p[i].h;
	}
	for(ll i=1;i<=cntx;i++){
		for(ll j=1;j<=cnty;j++){
			y=pre[1][i][j]%mod;
			y=(y+pre[1][i-1][j]*A.po[i]%mod+pre[1][i][j-1]*B.po[j]%mod)%mod;
			y=(y-pre[1][i-1][j-1]*A.po[i]%mod*B.po[j]%mod+mod)%mod,pre[1][i][j]=y;
		}
		for(ll j=cnty;j>=1;j--){
			y=pre[2][i][j]%mod;
			y=(y+pre[2][i-1][j]*A.po[i]%mod+pre[2][i][j+1]*B.po[j+1]%mod)%mod;
			y=(y-pre[2][i-1][j+1]*A.po[i]%mod*B.po[j+1]%mod+mod)%mod,pre[2][i][j]=y;
		}
	}
	for(ll i=cntx;i>=1;i--){
		for(ll j=1;j<=cnty;j++){
			y=pre[3][i][j]%mod;
			y=(y+pre[3][i+1][j]*A.po[i+1]%mod+pre[3][i][j-1]*B.po[j]%mod)%mod;
			y=(y-pre[3][i+1][j-1]*A.po[i+1]%mod*B.po[j]%mod+mod)%mod,pre[3][i][j]=y;
		}
		for(ll j=cnty;j>=1;j--){
			y=pre[4][i][j]%mod;
			y=(y+pre[4][i+1][j]*A.po[i+1]%mod+pre[4][i][j+1]*B.po[j+1]%mod)%mod;
			y=(y-pre[4][i+1][j+1]*A.po[i+1]%mod*B.po[j+1]%mod+mod)%mod,pre[4][i][j]=y;
		}
	}
	while(ops--){
		x=read(),y=read(),res=0;
		ll i=ub(lshx+1,lshx+1+cntx,x)-lshx-1,j=ub(lshy+1,lshy+1+cnty,y)-lshy-1;
		res=(res+pre[1][i][j]*A.qpow(abs(x-lshx[i]))%mod*B.qpow(abs(y-lshy[j]))%mod)%mod;
		res=(res+pre[2][i][j+1]*A.qpow(abs(x-lshx[i]))%mod*B.qpow(abs(y-lshy[j+1]))%mod)%mod;
		res=(res+pre[3][i+1][j]*A.qpow(abs(x-lshx[i+1]))%mod*B.qpow(abs(y-lshy[j]))%mod)%mod;
		res=(res+pre[4][i+1][j+1]*A.qpow(abs(x-lshx[i+1]))%mod*B.qpow(abs(y-lshy[j+1]))%mod)%mod;
		printf("%lld\n",res);
	}
	exit(0);
}

T3 是我的你不要抢

\(Hash\) 能水过.

正解还没打.

C_code
#include<bits/stdc++.h>
using namespace std;
namespace BSS {
	#define ll int
	#define ull unsigned ll
	#define lf long double
	#define lbt(x) (x&(-x))
	#define mp(x,y) make_pair(x,y)
	#define lb lower_bound 
	#define ub upper_bound
	#define Fill(x,y) memset(x,y,sizeof x)
	#define Copy(x,y) memcpy(x,y,sizeof x)
	#define File(x) freopen(#x".in","r",stdin),freopen(#x".out","w",stdout)
	inline ll read() {
		ll w=0; bool cit=1; char ch;
		while(!isdigit(ch=getchar())) if(ch=='-') cit=0; 
		while(isdigit(ch)) w=(w<<3)+(w<<1)+(ch^48),ch=getchar();
		return cit?w:-w;
	}
} using namespace BSS;

const ll N=6e5+21;
const ull P=131;

ll m,n,ops,maxlen,ans;
ll len[N];
ull po[N];
vector<ull> c[N],pre[N];
unordered_map<ll,ll> mp1[N];
inline ull getsh(ll x,ll l,ll r){ return pre[x][r]-pre[x][l-1]*po[r-l+1]; }
signed main(){
	File(string);
	n=read(),ops=read(); char s[N]; ll x,y;
	for(ll i=1;i<=n;i++){
		scanf("%s",s+1); c[i].push_back(0),pre[i].push_back(0);
		len[i]=strlen(s+1),maxlen=max(maxlen,len[i]);
		for(ll j=1;j<=len[i];j++){
			c[i].push_back(s[j]-'a'+1);
			pre[i].push_back(pre[i][j-1]*P+c[i][j]);
		}
	}
	po[0]=1; for(ll i=1;i<=maxlen;i++) po[i]=po[i-1]*P;
	while(ops--){
		y=read(),x=read(); maxlen=min(len[x],len[y]);
		if(mp1[x].find(y)!=mp1[x].end()){ printf("%lld\n",mp1[x][y]); continue; }
		for(ll i=maxlen;i>=0;i--){
			if(getsh(x,1,i)==getsh(y,len[y]-i+1,len[y])){
				ans=i; break;
			}
		}
		mp1[x][y]=ans;
		printf("%d\n",ans);
	}
	exit(0);
}

T4 显然也是我整的

还没会,先鸽.

这几场考得都不是很让自己满意.

每次考试的策略都很差.

这次考试上来先入手码了码量最长的题,甚至连 \(T1\) 和 \(T3\) 的暴力都没打.

\(T1\) 绝对在自己的能力范围内,自己知道是状压,可能推个规律就出来了.

可是自己偏偏就是看着 \(T2\) 去码较难的.

\(T4\) 连部分分都没拿全,这样的策略怎么可能行..

一定要规划好.

标签:pre,ch,noip,待补,res,ll,72,define,mod
来源: https://www.cnblogs.com/AaMuXiiiiii/p/15390750.html