2019.11.11【NOIP提高組】模擬 B 組

JZOJ 3888 正確答案

題目

試卷中共有$m$道判斷題，$n$個神犇中有$p$個考了滿分，$q$個考了零分，其他神犇不爲滿分或零分。還原出標準答案嗎，如有多解則輸出字典序最小的那個。無解輸出-1。

---

分析

分類討論，若$p$不爲0，枚舉正確答案，用$\text{STL::map}$維護，否則若$q$不爲0，枚舉錯誤答案，同樣用該數據結構維護，最後若$p,q$都爲0直接深搜找答案，時間複雜度不會太高的呀，也就是$O(nm)$左右吧

---

代碼

#include <iostream>
#include <string>
#include <map> 
#define rr register
using namespace std;
map<string,int>uk; int n,m,p,q,flag;
map<string,int>::iterator it;
inline string rev(string p){
	rr int len=p.size();
	for (rr int i=0;i<len;++i)
	    p[i]=p[i]=='N'?'Y':'N';
	return p;
}
inline void dfs(int dep,string str){
    if (flag) return;
	if (dep>m){
		if (uk.find(str)==uk.end()&&uk.find(rev(str))==uk.end()) flag=1,cout<<str;
		return;
	}
	dfs(dep+1,str+'N');
	dfs(dep+1,str+'Y');
} 
signed main(){
	ios::sync_with_stdio(0); string str;
	cin.tie(0),cout.tie(0),cin>>n>>m>>p>>q;
	for (rr int i=1;i<=n;++i) cin>>str,++uk[str];
	if (p){
		for (it=uk.begin();it!=uk.end();++it)
		if ((it->second)==p){
			rr string s1=rev(it->first);
			if (q&&(uk.find(s1)==uk.end()||uk[s1]!=q)) continue;
			cout<<it->first; flag=1;
	        break;
		}
	}else if (q){
		it=uk.end(); 
		do{
			--it;
			if (it->second==q){
				rr string s1=rev(it->first);
				if (uk.find(s1)!=uk.end()) continue;
				cout<<s1; flag=1;
				break;
			}
		}while (it!=uk.begin()); 
	}else dfs(1,str="");
	if (!flag) cout<<"-1";
	return 0;
} 

---

JZOJ 3889 序列問題

題目

在一個長度爲$n$的數列中，找出一段子序列$t1$，再在其後面找到一段子序列$t2$，使$t1$的異或和與$t2$的按位與相等，問有多少種找出$t1,t2$的方案

---

分析

由於它們相等，所以總異或和爲0，考慮dp（這裏的數列要倒置）
設$dp[i][j][k]$表示前$i$個數進行的結果爲$j$，$k=0$時表示絲毫未取，$k=1$時表示取按位與，$k=2$時表示取異或和，最後輸出$dp[n][0][2]$
$dp[i][j][0]=dp[i-1][j][0]$
$dp[i][j\: and\: a[i]][1]=dp[i-1][j][0]+dp[i-1][j][1]+dp[i-1][j\: and\: a[i]][1]$
$dp[i][j]\: xor\: a[i]][2]=dp[i-1][j][1]+dp[i-1][j][2]+dp[i-1][j\: xor\: a[i]][2]$
然後要滾動數組+高精度

---

代碼

#include <cstdio>
#include <cctype>
#include <cstring>
#define rr register
using namespace std;
const int N=1031; const long long mod=1e17; int a[N],n;
struct Bignum{int len; long long p[21];}dp[2][N][3];
inline signed iut(){
	rr int ans=0; rr char c=getchar();
    while (!isdigit(c)) c=getchar();
    while (isdigit(c)) ans=(ans<<3)+(ans<<1)+(c^48),c=getchar();
    return ans;
}
inline void print(long long ans){
	if (ans>9) print(ans/10);
	putchar(ans%10+48);
}
inline void plus(Bignum &t1,const Bignum t2){
	rr int len=t1.len>t2.len?t1.len:t2.len,g=0;
	for (rr int i=1;i<=len;++i){
	    rr long long s=t1.p[i]+t2.p[i]+g;
	    t1.p[i]=s>=mod?s-mod:s,g=s>=mod;
	}
	if (g) t1.p[++len]=g; t1.len=len;
}
signed main(){
	rr Bignum *t=&dp[0][1023][0];
	n=iut(),(*t).p[(*t).len=1]=1;
	for (rr int i=n;i>=1;--i) a[i]=iut();
	for (rr int i=1;i<=n;++i){
		memcpy(dp[i&1],dp[(i&1)^1],sizeof(dp[(i&1)^1]));
		for (rr int j=0;j<1024;++j){
			t=&dp[i&1][j&a[i]][1];
			plus(*t,dp[(i&1)^1][j][0]);
			plus(*t,dp[(i&1)^1][j][1]);
			t=&dp[i&1][j^a[i]][2];
			plus(*t,dp[(i&1)^1][j][1]);
			plus(*t,dp[(i&1)^1][j][2]);
		}
	}
	t=&dp[n&1][0][2];
	print((*t).p[(*t).len]);
	for (rr int i=(*t).len-1;i>=1;--i){
		rr long long now=(*t).p[i];
		while (now*10<mod) now*=10,putchar(48);
		print((*t).p[i]);
	}
	return 0;
}

---

JZOJ 3890 長途旅行

弱化版，此題邊權較大，不適宜矩陣乘法

---

分析

考慮答案應爲鏈+環的總和，而環是可以增加的，環必須經過0號點，所以可以設0號點中的一條邊爲循環節，設$dis[x][now]$表示到第$x$個點時模循環節的答案爲$now$時的路徑長度，顯然只要$dis[n-1][T\bmod m]\leq T$那麼必須有解

---

代碼

#include <cstdio>
#include <cctype>
#include <cstring>
#include <queue>
#define rr register
using namespace std;
struct node{int y,w,next;}e[211];
typedef long long ll; ll dis[51][10011],T;
int n,m,k,ls[51],mn; bool v[51][10011];
inline ll iut(){
	rr ll ans=0; rr char c=getchar();
	while (!isdigit(c)) c=getchar();
	while (isdigit(c)) ans=(ans<<3)+(ans<<1)+(c^48),c=getchar();
	return ans;
}
signed main(){
	for (rr int t=iut();t;--t){
		n=iut(),m=iut(),T=iut(),k=1;
		memset(ls,0,sizeof(ls)),mn=20001;
		for (rr int i=1;i<=m;++i){
			rr int x=iut()+1,y=iut()+1,w=iut();
			e[++k]=(node){y,w,ls[x]},ls[x]=k,
			e[++k]=(node){x,w,ls[y]},ls[y]=k;
			if (x==1||y==1) mn=mn<w?mn:w;
		}
		mn<<=1; memset(dis,42,sizeof(dis));
		rr queue<pair<int,int> >q; dis[1][0]=0;
		q.push(make_pair(1,0)); v[1][0]=1;
		while (q.size()){
			rr int x=q.front().first,now=q.front().second; q.pop();
			for (rr int i=ls[x];i;i=e[i].next){
			    rr int stat=(now+e[i].w)%mn;
				if (dis[e[i].y][stat]>dis[x][now]+e[i].w){
			        dis[e[i].y][stat]=dis[x][now]+e[i].w;
			        if (!v[e[i].y][stat]){
			    	    v[e[i].y][stat]=1;
			    	    q.push(make_pair(e[i].y,stat));
				    }
			    }
		    }
			v[x][now]=0;
		}
		if (dis[n][T%mn]<=T) printf("P");
		    else printf("Imp");
		printf("ossible\n");
	}
	return 0;
}