2020九月&十月暂记

最新推荐文章于 2021-10-12 18:17:52 发布
最新推荐文章于 2021-10-12 18:17:52 发布
阅读量115 收藏
点赞数
分类专栏： 题解 临时暂记 文章标签： 1024程序员节
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/bell041030/article/details/108951366
版权

Sept.&Oct.

新&重&难点：

  1. 状态压缩动态规划：(网格中)逐行状态转移 ;逐格状态转移(轮廓线优化)
  2. 计算几何：扫描线(+线段树)；半平面交；曼哈顿距离与切比雪夫距离转换
  3. 数学：二项式定理与二项式反演；概率期望；逆元求 t × a − 1 m o d    b t\times a^{-1}\mod b t×a−1modb
  4. 图论综合：最短路为主；(LCA、生成树、缩点……)
  5. 数据结构综合：线段树为主
  6. 大赛后希望学习：四边形不等式优化区间dp，斜率优化dp，强连通分量……更多学习计划写在11月暂记里吧。。。

2020.9.4

2020.9.4
T1
找出所有能到达蹦床的位置，将这些位置全部打上标记
如果当前出发点在能到达蹦床的位置，把打上标记的位置全部统计一遍
最后再找一段最大连续上升或最大连续下降序列并统计其中的楼房数便得到最后的答案
如果当前不在能到达蹦床的位置，统计从当前开始的一段最大连续上升或下降序列即为答案

#include<iostream>
#include<string>
#include<cstring>
#include<cstdio>

using namespace std;
int N,K,h[300005],jump[300005],ans,cnt,b[300005],sum[300005],tmp,lt,rt;
string s;
bool start; 

int main()
{
	cin>>N>>K;
	for(int i = 1;i <= N;i ++)
		cin>>h[i];
	cin>>s;
	for(int i = 0;i < s.size();i ++)
		if(s[i]=='T')
		{
			cnt++;
			jump[cnt] = i+1;
		}
	if(cnt != 0)//有蹦床的情况 
	{
		for(int i = 1;i <= cnt;i ++)//给能到蹦床上的区间打上标记 
		{
			if(b[jump[i]] != 0) continue;
			sum[i] = 1;
			b[jump[i]] = i;
			for(int j = jump[i] - 1 ;j >= 1;j --)
			{
				if(h[j]>=h[j+1]&&b[j] == 0)
				{
					b[j] = i,sum[i] ++;
				}
				else break;
			}	
			for(int j = jump[i] + 1;j <= N;j ++)
			{
				if(h[j]>=h[j-1]&&b[j] == 0)
				{
					b[j] = i,sum[i] ++;
				}
				else break;
			}
		}
		if(b[K]!=0)//当前可以到蹦床那么就可以把所有可到达蹦床区间的楼房都跳一遍 
		{
			tmp = 0;int maxn = -1;
			for(int i = 1;i <= cnt;i ++) ans += sum[i];
			if(b[N-1] == 0) tmp = 1;
			for(int i = N-1;i >= 1;i --)
			{
				if(h[i]>=h[i+1] && b[i] == 0) tmp ++;
				else
				{
					maxn = max(maxn,tmp);
					tmp = 0;
					if(b[i] == 0) tmp ++;
				}
			}
			maxn = max(maxn,tmp);
			tmp = 0;
			if(b[N-1] == 0) tmp = 1;
			for(int i = N-1;i >= 1;i --)
			{
				if(h[i]<=h[i+1] && b[i] == 0) tmp ++;
				else
				{
					maxn = max(maxn,tmp);
					tmp = 0;
					if(b[i] == 0) tmp ++;
				}
			}
			maxn = max(maxn,tmp);
			tmp = 0;
			if(b[1] == 0) tmp ++;
			for(int i = 2;i <= N;i ++)
			{
				if(h[i]<=h[i-1] && b[i] == 0) tmp ++;
				else
				{
					maxn = max(maxn,tmp);
					tmp = 0;
					if(b[i] == 0) tmp ++;
				}
			}
			maxn = max(maxn,tmp);
			tmp = 0;
			if(b[1] == 0) tmp ++;
			for(int i = 2;i <= N;i ++)
			{
				if(h[i]>=h[i-1] && b[i] == 0) tmp ++;
				else
				{
					maxn = max(maxn,tmp);
					tmp = 0;
					if(b[i] == 0) tmp ++;
				}
			}
			maxn = max(maxn,tmp);
			ans = maxn + ans;
			cout << ans << endl;
			return 0;
		}
		
	}
	//没有蹦床或是当前无法去到任何一个蹦床
	//统计楼房高度相平横向跳跃的总数 
	lt = rt = K;
	for(int i = K-1;i >= 1;i --)
		if(h[i] != h[i+1])
		{
			lt = i + 1;
			break;
		}
	for(int i = K+1;i <= N;i ++)
		if(h[i] != h[i-1])
		{
			rt = i - 1;
			break;
		}
	//
	//
	tmp = 0;
	for(int i = lt-1;i >= 1;i --)
	{
		if(h[i]<=h[i+1]) tmp ++;
		else
		{
			ans = max(ans,tmp);
			break;
		}
	}
	ans = max(ans,tmp);
	tmp = 0;
	for(int i = rt+1;i <= N;i ++)
	{
		if(h[i]<=h[i-1]) tmp ++;
		else
		{
			ans = max(ans,tmp);
			break;
		}
	}
	ans = max(ans,tmp);
	ans += rt - lt + 1;
	cout << ans << endl;
	return 0;
}
/*
10 1
10 7 3 1 1 9 8 2 4 10
..T..T.T..
*/

2020.9.9

2020.9.9

  • 收录整理

T2

2020.9.30

2020.9.30
T1

  • 收录整理

当所有字符串变换成使得它们的公共前缀最长时，trie树的结点树最少
看到数据范围 1 ≤ n ≤ 16 1\le n\le 16 1≤n≤16，自然先想到状压DP
由于字符串可以任意变换，只需统计字符串中每一种字母出现的个数就可以方便表示字符串，且方便接下来的操作
设 S S S表示被选取加入trie树的字符串的集合， S ′ ⊆ S S'\subseteq S S′⊆S， f ( S ) f(S) f(S)表示在状态 S S S下trie树的最少结点数， p ( S ) p(S) p(S)表示在 S S S状态下，最大公共前缀的长度
则有： f ( S ) = f ( S ′ ) + f ( S − S ′ ) − p ( S ) f(S) = f(S')+f(S-S')-p(S) f(S)=f(S′)+f(S−S′)−p(S)

#include<iostream>
#include<string>
#include<cstring>
#include<algorithm>
#include<fstream>

using namespace std;
const int inf = 0xfffffff;
int n,f[65540],tmp,tot,cnt[65540][30];
string s[25];

int main()
{
//	freopen("test.txt","r",stdin);
	cin >> n;
	tot = (1<<n)-1;
	for(int i = 0; i <= tot;i ++) f[i] = inf;
	for(int i = 0;i <= tot;i ++)
		for(int j = 0;j < 26;j ++)
			cnt[i][j] = inf ;
	for(int i = 1;i <= n;i ++)
	{
		cin >> s[i];
		tmp = 1<<(i-1);
		f[tmp] = s[i].size();
		for(int j = 0;j < 26;j ++) cnt[tmp][j] = 0;
		for(int j = 0;j < s[i].size();j ++) cnt[tmp][s[i][j]-'a'] ++;
	} 
	f[0] = 0;
	for(int i = 1;i <= tot;i ++)
	{
		tmp = 0;
		for(int j = 1;j <= n;j ++)
			if(((1<<(j-1))&i) > 0)
				for(int k = 0;k < 26;k ++)
					cnt[i][k] = min(cnt[i][k],cnt[(1<<(j-1))][k]);
		for(int k = 0;k < 26;k ++)
			tmp += cnt[i][k];
		for(int j = i&(i-1);j;j = (j-1)&i)
		{
			int k = i - j;
			f[i] = min(f[i],f[k]+f[j]-tmp);
		}
	}
	cout << f[tot]+1 << endl;
	return 0;
}

T2
在这里插入图片描述
斜边所在直线 y = − x + b y=-x+b y=−x+b由图中显然可知： y = − x + x ′ + y ′ + r ′ y = -x+x'+y'+r' y=−x+x′+y′+r′
故求两斜边中较下的一条，比较 b b b即 x ′ + y ′ + r ′ x'+y'+r' x′+y′+r′即可

/*
    |
    |\
 x  |__\
 +  | x |\
 r  |   |  \ 
 +  |___|____\
 y  |_y_|__________
 
 */
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cmath>

using namespace std;
const long long inf = 0xfffffffff;
int n,tot,cnt,times[1105];
long long nx,ny,nc,x[1105],y[1105],c[1105];
double ans;

int main()
{
	cin >> n;
	tot = (1<<n)-1;
	for(int i = 1;i <= n;i ++)
		cin >> x[i] >> y[i] >> c[i];
	for(int i = 1;i <= tot;i ++)
	{
		cnt = 0;nc = inf;nx = -inf;ny = -inf;
		for(int j = 1;j <= n;j ++)
			if(((1<<(j-1))&i) > 0)
			{
				nc = min(nc,x[j]+y[j]+c[j]);
				nx = max(nx,x[j]);
				ny = max(ny,y[j]);
				cnt ++;
			}
		ans += pow(-2,cnt-1)*max((long long)0,nc-nx-ny)*max((long long)0,nc-nx-ny)*0.5;
	}
	printf("%.1lf",ans);
	return 0;
}

10.6模拟赛

  • 收录整理
    在这里插入图片描述
    又是状压DP
#include<iostream>
#include<cstdio>

using namespace std;
const int modnum = 1000000007;
int n,m,k,tot,s1;
int u[505],v[505],tmp;
bool map[30][30];
long long f[35][1048580],ans;

int main()
{
	scanf("%d%d%d",&n,&m,&k);
	tot = (1 << n) - 1;
	for(int i = 1;i <= m;i ++)
	{
		scanf("%d%d",&u[i],&v[i]);
		map[u[i]][v[i]] = 1;
	}
	f[0][0] = 1;
	for(int s = 1; s <= tot; s ++)//enumerate all the states
		for(int i = 1;i <= n;i ++)//enumerate the last number in this sequence
		{
			tmp = 0;
			if(((1<<(i-1))&s) > 0)
			{
				s1 = s - (1<<(i-1));
				
				for(int j = 1;j <= n;j ++)
					if(((1<<(j-1))&s1) > 0)
					{
						if(map[i][j]) tmp ++;//if the last number is i , it will cause tmp arguements
						if(tmp > k) break;
					}
				if(tmp > k) continue;
				for(int j = 0;j <= (k-tmp);j ++)//enumerate last status
					f[tmp+j][s] = (f[tmp+j][s]%modnum + f[j][s1]%modnum)%modnum;
			}
		}
	for(int i = 0;i <= k;i ++)//sum up the answer
		ans = (ans%modnum + f[i][tot]%modnum)%modnum;
	printf("%lld",(ans%modnum));
	return 0;
}

恶补状压DP

T1

#include<cstdio>
#include<iostream>
#include<cmath>

using namespace std;
const int modnum = 100000000;
int n,m,map[20][20],f[20][20][32780];
int st[20];
bool vis[20][20][32780];

int dfs(int i,int j,int s)
{
	if(i>n) return 1;//全部放置完毕,算作一种方案 
	if(j>m) return dfs(i+1,1,s&st[i+1])%modnum;//this row has been counted,go to the next grid
	if(vis[i][j][s]) return f[i][j][s]%modnum;//we've known the answer and return it
	vis[i][j][s] = 1;//
	int now = 1<<(j-1);//the status of (i,j)
	if((s&now)==0)//we can't place a plant here
	{
		int s1 = (s|now);//we can't place a plant on (i,j) , but we can place a plant on (i+1,j)
		f[i][j][s] = dfs(i,j+1,s1)%modnum;//go to the next grid
		return f[i][j][s]%modnum;
	}
	int tmp = dfs(i,j+1,s)%modnum;//we don't place a plant here
	int s2 = s - now;int next = 1<<j;//if we place a plant here,we can't place a plant on the next grid
	if(((next&s2)>0)&&(j<m)) s2 -= next;//figure out the status of placing a plant here
	int tmp2 = dfs(i,j+1,s2)%modnum;//go to the next grid
	f[i][j][s] = (tmp+tmp2)%modnum;//sum up the answer
//	cout << i << ',' << j << ':' << s << " plant" << tmp << " empty" << tmp2 << endl;
	return f[i][j][s];
}

int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++) st[i]=0;
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=m;j++)
		{
			scanf("%d",&map[i][j]);
			if(map[i][j])
				st[i]+=pow(map[i][j]*2,j-1);
		}
	}
		
	printf("%d",dfs(1,1,st[1])%modnum);
	
	return 0;
}

T2

#include<iostream>
#include<cstdio>
#include<algorithm>

using namespace std;
int n,map[15][15];
int f[2048],tot;

int main()
{
	while(1)
	{
		scanf("%d",&n);
		if(n == 0) break;
		tot = (1<<n) - 1;
		for(int i = 1;i <= n;i ++)
			for(int j = 1;j <= n;j ++)
				scanf("%d",&map[i][j]);
		for(int i = 0;i <= tot;i ++) f[i] = 0;
		for(int k = 1;k <= tot;k ++)
			for(int i = 0;i < n;i ++)
			{
				int s1 = 1<<i;
				if((s1&k)>0)
					for(int j = i + 1;j < n;j ++)
					{
						int s2 = 1<<j;
						if((s2&k)>0)
						{
							int s3 = k - s1;
							int s4 = k - s2;
							f[k] = max(f[k],f[s3] + map[j+1][i+1]);
							f[k] = max(f[k],f[s4] + map[i+1][j+1]);
						}
					}
			}
		printf("%d\n",f[tot]);
	}
	return 0;
}

T3
T4

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<queue>

using namespace std;
const long long inf = 0xfffffffff;
typedef long long ll;
ll n,m,t,d[55][55],dis[55][55],f[66000][25];
int map[55][55],kx,ky,qx,qy,gx[55],gy[55],cnt,s1,ans,tot,num1;
int mark[55][55],nowx,nowy,nx,ny,dx[5] = {0,0,0,1,-1},dy[5] = {0,1,-1,0,0};
bool vis[55][55];
string s;

int cnt1(int x)
{
	int tmp = 0;
	for(int i = 1;i <= cnt;i ++)
		if((x&(1<<(i-1))) > 0)
			tmp ++;
	return tmp;
}

void bfs(int sx,int sy)
{
	memset(vis,0,sizeof(vis));
	for(int i = 1;i <= n;i ++)
		for(int j = 1;j <= m;j ++)
			d[i][j] = inf;
	queue< pair<int , int > > q;
	vis[sx][sy] = 1;
	d[sx][sy] = 0;
	q.push(make_pair(sx,sy));
	while(!q.empty())
	{
		nowx = q.front().first;
		nowy = q.front().second;
		q.pop();
		vis[nowx][nowy] = 0;
		for(int i = 1;i <= 4;i ++)
		{
			nx = nowx + dx[i];ny = nowy + dy[i];
			if(nx > 0 && nx <= n && ny > 0 && ny <= m )
				if(d[nx][ny] > d[nowx][nowy] + 1 && map[nx][ny] == 1)
				{
					d[nx][ny] = d[nowx][nowy] + 1;
					if(!vis[nx][ny])
					{
						vis[nx][ny] = 1;
						q.push(make_pair(nx,ny));
					}
				}
		}
	}
	for(int i = 0;i <= cnt + 1;i ++)
		dis[mark[sx][sy]][i] = d[gx[i]][gy[i]];
	return ;
}


int main()
{
	memset(mark,128,sizeof(mark));
	//read data
	scanf("%lld%lld%lld",&n,&m,&t);
	for(int i = 1;i <= n;i ++)
	{
		cin >> s;
		for(int j = 0;j < s.size();j ++)
		{
			switch(s[j])
			{
				case '.':
					map[i][j+1] = 1;
					break;
				case '#':
					map[i][j+1] = 0;
					break;
				case 'G':
					cnt ++;
					gx[cnt] = i;gy[cnt] = j + 1;
					mark[i][j+1] = cnt;
					map[i][j+1] = 1;
					break;
				case 'K':
					kx = i;ky = j + 1;
					map[i][j+1] = 1;
 					break;
				case 'Q':
					qx = i;qy = j + 1;
					map[i][j+1] = 1;
					break;
			}
		}
	}
	mark[kx][ky] = 0;mark[qx][qy] = cnt + 1;
	gx[0] = kx;gy[0] = ky;
	gx[cnt + 1] = qx;gy[cnt + 1] = qy;
	for(int i = 0;i <= cnt+1;i ++)
		for(int j = 0;j <= cnt + 1;j ++)
			dis[i][j] = inf;
			
	for(int i = 0;i <= cnt + 1;i ++)
		bfs(gx[i],gy[i]);
	
	tot = 1<<(cnt+1) - 1;
	for(int i = 0;i <= tot;i ++)
		for(int j = 0;j <= cnt+1; j ++)
			f[i][j] = inf;
	for(int i = 1;i <= cnt;i ++)
		f[(1<<(i-1))][i] = dis[0][i];
	for(int s = 0;s <= tot; s ++)
	{
		num1 = cnt1(s);
		for(int i = 1;i <= cnt;i ++)
			if((s&(1<<(i-1))) > 0)
			{
				s1 = s - (1<<(i-1));
				for(int j = 1;j <= cnt;j ++)
					if((s1&(1<<(j-1))) > 0)
						f[s][i] = min(f[s][i],f[s1][j]+dis[j][i]*num1);
			}
	}
	for(int s = 0;s <= tot;s ++)
	{
		num1 = cnt1(s);
		for(int i = 1;i <= cnt;i ++)
			if((s&(1<<(i-1))) > 0)
			{
				ll sum = f[s][i] + dis[i][cnt+1]*(num1+1);
				if(sum <= t)
					ans = max(ans,num1);
			}
	}
	printf("%d",ans);
	return 0;
}

T5

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<queue>

using namespace std;
const int inf = 2139062143;
int n,cp,k,r,u,v,cnt,f[45][1050],tot,tmp,s1,num;
int map[15][45][45],p[15][45],d[50];
bool vis[50];

int cnt1(int x)
{
	int tmp = 0;
	for(int i = 1;i <= cp;i ++)
		if((x&(1<<(i-1))) > 0)
			tmp ++;
	return tmp;
}

void dfs(int c,int now,int fa)
{
	if(vis[now]) return ;
	vis[now] = 1;
	for(int i = 1;i <= n;i ++)
		if(map[c][now][i])
		{
			map[c][fa][i] = 1;
			dfs(c,i,fa);
		}
}

void solve()
{
	//init
	memset(d,127,sizeof(d));
	memset(f,127,sizeof(f));
	memset(map,0,sizeof(map));
	memset(p,0,sizeof(p));
	
	//read
	scanf("%d%d",&n,&cp);
	for(int i = 1;i <= cp;i ++)
	{
		scanf("%d",&k);
		for(int j = 1;j <= k;j ++)
		{
			scanf("%d%d",&u,&v);
			map[i][u][v] = 1;
		}
	}
	for(int i = 0;i < n;i ++)
		for(int j = 1;j <= cp;j ++)
			scanf("%d",&p[j][i]);
	
	//dp
	tot = 1<<(cp+1) - 1;
	for(int i = 0;i <= tot;i ++)
	{
		num = cnt1(i);
		if(num > 3) continue;
		tmp = 0 ;
		for(int j = 1;j <= cp;j ++)
			if((i&(1<<(j-1))) > 0)
				tmp += p[j][0];
		f[0][i] = tmp;
	}
	
	memset(vis,0,sizeof(vis));
	queue<int > q; 
	q.push(0);
	while(!q.empty())
	{
		int i = q.front();
		q.pop();
		vis[i] = 0;
		for(int s = 0;s <= tot;s ++)//enumrate the status of tickets
		{
			num = cnt1(s);
			if(num > 3) continue;
			for(int j = 1;j <= cp;j ++)//We buy tickets on this island
				if((s&(1<<(j-1))) > 0)
				{
					s1 = s - (1<<(j-1));
					f[i][s] = min(f[i][s],f[i][s1]+p[j][i]);
				}
			for(int j = 1;j <= cp;j ++)//We have bought tickets on other island
				if((s&(1<<(j-1))) > 0)
				{
					s1 = s - (1<<(j-1));
					for(int k = 0;k < n;k ++)
						if(map[j][i][k])//Use this ticket to travel from i to k
						{
							d[k] = min(d[k],f[i][s]);//update the minimum distance
							if(f[k][s1] > f[i][s])//update the next status
							{
								f[k][s1] = f[i][s];
								if(!vis[k])
								{
									vis[k] = 1;
									q.push(k);//keep updating
								}
							}
						}
				}
		}
	}

	for(int i = 1;i < n;i ++)
	{
		if(d[i] == inf) printf("-1 ");
		else printf("%d ",d[i]);
	}
	printf("\n");
}

int main()
{
	scanf("%d",&r);
	for(int i = 1;i <= r;i ++)
		solve();

	return 0;
}

T6

#include<iostream>
#include<cstdlib>
#include<string>
#include<cstring>

using namespace std;
int h,w,tot;
long long f[25][25][30000];
bool vis[25][25][30000];

long long dfs(int x,int y,int s)
{
	if(x > h) return 1;//this rectangle is full of bricks 
	if(y > w) return dfs(x+1,1,s);//next row
	if(vis[x][y][s]) return f[x][y][s];
	vis[x][y][s] = 1;
	int s1 = 1<<(y-1);
	if((s1&s)>0) //the current grid is occupied by a vertical brick
	{
		f[x][y][s] = dfs(x,y+1,s-s1);// As for the next row,it won't be affected by this vertical brick 
		return f[x][y][s];
	}
	long long tmp1 = 0;
	if(x < h)//place a vertical brick on this grid
		tmp1 = dfs(x,y+1,s1+s);//the placement of this brick will influence the way of placing bricks on the next row
	int s2 = 1<<y;long long tmp2 = 0;
	if(y < w)//place a transverse brick on this grid
		if((s2&s)==0)//to check whether we can place this transverse brick or not
			tmp2 = dfs(x,y+1,s+s2);
	f[x][y][s] = tmp1 + tmp2;
	return f[x][y][s];
}

int main()
{
	while(1)
	{
		cin>>h>>w;
		if(h == 0 && w == 0) break;
		tot = (1<<(w+1))-1;
		memset(f,0,sizeof(f));
		memset(vis,0,sizeof(vis));
		cout << dfs(1,1,0) << endl;
	}
	return 0;
}

T7

#include<iostream>
#include<cstdio>
#include<string>
#include<cstring>

using namespace std;
const int inf = 2139062143;
int n,m,f[15][1050],tmp,tmp1,ans = inf,tot,b[15];
//1 refers to drawing vertically and 0 refers to drawing transversely
//regard # as 0 
bool map[15][15],flag;
string s;

bool check(int now,int x)//the position of # in the map can't be 1 in the status
{
	if((b[now]&x) > 0) return 0;
	else return 1;
}

int main()
{
	memset(f,127,sizeof(f));
	scanf("%d%d",&n,&m);
	for(int i = 1;i <= n;i ++)
	{
		cin >> s;
		for(int j = 0;j < s.size();j ++)
			if(s[j] == '#')
			{
				map[i][j+1] = 1;
				b[i] += (1<<j);//record the status
			}
	}
	
	tot = 1<<(m+1) - 1;
	for(int s = 0;s <= tot; s++)//initalization
	{
		if(!check(1,s)) continue;
		flag = 0;tmp = 0;
		for(int j = 1;j <= m;j ++)
		{
			if((s&(1<<(j-1))) > 0)//draw vertically
			{
				tmp ++;
				flag = 0;
			}
			else
			{
				if(map[1][j])
				{
					flag = 0;
					continue;
				}//draw transeversely
				if(flag) continue;//draw with one stroke 
				else
				{
					flag = 1;
					tmp ++;
				}
			}
		}
		f[1][s] = tmp;
		if(1 == n) ans = min(ans,f[1][s]);
	}
	
	for(int i = 2;i <= n;i ++)
	{
		for(int s = 0;s <= tot;s ++)//the status of last row
		{
			if(!check(i-1,s)) continue;
			for(int s1 = 0;s1 <= tot;s1 ++)//the status of this row
			{
				if(!check(i,s1)) continue;
				flag = 0;tmp = 0;tmp1 = 0;
				for(int j = 1;j <= m;j ++)
				{
					if((s1&(1<<(j-1))) > 0)
					{
						flag = 0;
						//draw transeversely on the last row
						if((s&(1<<(j-1))) == 0) tmp1 ++;//draw vertically this row
					}
					else
					{
						if(map[i][j])
						{
							flag = 0;
							continue;
						}
						if(flag) continue;
						else flag = 1,tmp ++;
					}
				}
				f[i][s1] = min(f[i][s1],f[i-1][s] + tmp1 + tmp);
//				cout << i << ' ' << s1 << ' ' << f[i][s1] << endl;
				if(i == n) ans = min(ans,f[i][s1]);
			}
		}
	}
	if(ans == inf) printf("0");
	else printf("%d",ans);
	return 0;
}

10.23

10.23模拟赛

T1
通过探索规律可以发现：只要算出2天内团伙的发展情况，此后的发展情况全部都由 t = 0 t=0 t=0， t = 1 t=1 t=1， t = 2 t=2 t=2的发展情况表示
我们可以算出：
t = 0 , f ( 0 ) = B t = 1 , f ( 1 ) = SSB t = 2 , f ( 2 ) = BBSSB \begin{array}{crl} t=0,f(0) = & \text {B}\\ t=1,f(1)= &\text {SSB} \\ t=2, f(2)= &\text {BBSSB} \\ \end{array} t=0,f(0)=t=1,f(1)=t=2,f(2)=​BSSBBBSSB​
再推几个式子发现：
t = 3 , f ( 3 ) = SSBSSBBBSSB t = 4 , f ( 4 ) = BBSSBBBSSBSSBSSBBBSSB ∵ f ( 3 ) = f ( 1 ) + f ( 1 ) + f ( 0 ) + f ( 0 ) + f ( 1 ) f ( 4 ) = f ( 2 ) + f ( 2 ) + f ( 1 ) + f ( 1 ) + f ( 2 ) ∴ 可 推 广 得 : f ( i ) = { B , i = 0 S S B , i = 1 B B S S B , i = 2 f ( i − 2 ) + f ( i − 2 ) + f ( i − 1 ) + f ( i − 1 ) + f ( i − 2 ) , i > 2 \begin{array}{crl} t=3, f(3) = & \text {SSBSSBBBSSB} \\ t=4, f(4) = & \text {BBSSBBBSSBSSBSSBBBSSB}\\ \end{array}\\ \begin{array}{l} \because \\ f(3) = f(1)+f(1)+f(0)+f(0)+f(1)\\ f(4)=f(2)+f(2)+f(1)+f(1)+f(2) \\ \therefore可推广得:\\ f(i) = \begin{cases} B &,i = 0 \\ SSB &,i=1\\ BBSSB&,i=2\\ f(i-2)+f(i-2)+f(i-1)+f(i-1)+f(i-2)&,i>2 \end{cases} \end{array} t=3,f(3)=t=4,f(4)=​SSBSSBBBSSBBBSSBBBSSBSSBSSBBBSSB​∵f(3)=f(1)+f(1)+f(0)+f(0)+f(1)f(4)=f(2)+f(2)+f(1)+f(1)+f(2)∴可推广得:f(i)=⎩⎪⎪⎪⎨⎪⎪⎪⎧​BSSBBBSSBf(i−2)+f(i−2)+f(i−1)+f(i−1)+f(i−2)​,i=0,i=1,i=2,i>2​​
同时容易发现，设发展到第 x x x天 B B B的个数为 g ( x ) g(x) g(x)，全部字符个数为 φ ( x ) \varphi(x) φ(x)，则有：
g ( x ) = { 1 , x = 0 2 g ( x − 1 ) − 1 , x > 0 , x m o d    2 = 1 2 g ( x − 1 ) + 1 , x > 0 , x m o d    2 = 0 φ ( x ) = { 1 , x = 0 2 φ ( x − 1 ) + 1 , x > 0 , x m o d    2 = 1 2 φ ( x − 1 ) − 1 , x > 0 , x m o d    2 = 0 g(x)= \begin{cases} 1&,x = 0 \\ 2g(x-1)-1&,x>0,x\mod 2 = 1 \\ 2g(x-1)+1&,x>0,x\mod 2 = 0\\ \end{cases}\\ \\ \varphi (x) = \begin{cases} 1 &,x=0 \\ 2\varphi(x-1)+1&,x>0,x\mod 2 = 1 \\ 2\varphi(x-1)-1&,x>0,x\mod 2 = 0\\ \end{cases}\\ g(x)=⎩⎪⎨⎪⎧​12g(x−1)−12g(x−1)+1​,x=0,x>0,xmod2=1,x>0,xmod2=0​φ(x)=⎩⎪⎨⎪⎧​12φ(x−1)+12φ(x−1)−1​,x=0,x>0,xmod2=1,x>0,xmod2=0​
然后就能推算出间谍在哪一个区间，递归求其位置即可

#include<iostream>
#include<cstdio>

using namespace std;
int n,k,cnt,tmp;
int len[35],p1,p2,p3,p4;
int num[35],ans;

int l(int x)
{
	if(x%2 == 0) return -1;
	else return 1;
}

void solve(int pre,int sum,int m)//递归
{
	if(m == 0)
	{
		ans = pre+1;
		return ;
	}
	if(m == 1)
	{
		ans = pre+3;
		return ;
	}
	if(m == 2)
	{
		k -= sum;
		if(k == 1) ans = pre+1;
		if(k == 2) ans = pre+2;
		if(k == 3) ans = pre+5;
		return ;
	}
	int tmp1 = num[m-2];int pre1 = len[m-2];
	int tmp2 = tmp1+tmp1;int pre2 = pre1+pre1;
	int tmp3 = tmp2+num[m-3];int pre3 = pre2+len[m-3];
	int tmp4 = tmp3+num[m-3];int pre4 = pre3+len[m-3];
	if(k<=sum+tmp1)
	{
		solve(pre,sum,m-2);
		return ;
	}
	if(k<=sum+tmp2)
	{
		solve(pre+pre1,sum+tmp1,m-2);
		return ;
	}
	if(k<=sum+tmp3)
	{
		solve(pre+pre2,sum+tmp2,m-3);
		return ;
	}
	if(k<=sum+tmp4)
	{
		solve(pre+pre3,sum+tmp3,m-3);
		return ;
	}
	solve(pre+pre4,sum+tmp4,m-2);
	return ;
}

int main()
{
	len[0] = 1;num[0] = 1;num[1] = 1;
	for(int i = 1;i <= 30;i ++)
	{
		len[i] = len[i-1]*2+l(i);//phi(x)
		num[i+1] = len[i];//g(x)
	}
	scanf("%d%d",&n,&k);
	solve(0,0,n);
	printf("%d",ans);
	return 0;
}

T4未改好

#include<iostream>
#include<cstdio>

using namespace std;
const long long modnum = 1000000007;
typedef long long ll;
ll n,a[105],f[105][1005][105],p[105],sum,sum2,ans,ti[1005];
int len[105],tmp,num,maxlen;

int main()
{
	ti[0] = 1;
	for(int i = 1;i <= 1000;i ++)
	{
		ti[i] = ((ti[i-1]%modnum)*10)%modnum;
		ti[i] %= modnum;
	}
	scanf("%lld",&n);
	for(int i = 1;i <= n;i ++)
	{
		scanf("%lld",&a[i]);
		num = a[i];tmp = 0;
		while(num)
		{
			tmp ++;
			num /= 10;
		}
		len[i] = tmp;
		maxlen += tmp;
	}
	p[0] = 1;
	for(int i = 1;i <= n;i ++)
	{
		p[i] = ((p[i-1]%modnum)+(((p[i-1]%modnum)*(i%modnum))%modnum))%modnum;
		p[i] = p[i]%modnum;
	}
	
	for(int i = 1;i <= n;i ++)
		f[i][0][0] = p[n-1];
		
	for(int k = 1;k < n;k ++)
		for(int i = 1;i <= n;i ++)
			for(int j = k;j <= maxlen-len[i];j ++)
				for(int ii = 1;ii <= n;ii ++)
				{
					if(i!=ii&&len[ii] <= j)
					{
						sum = ((p[n-k-1]%modnum)*(f[ii][j-len[ii]][k-1]%modnum))%modnum;
						sum2 = 0;
						f[i][j][k] += sum%modnum;
						f[i][j][k] %= modnum;
					}
				}
					
	
	for(int i = 1;i <= n;i ++)
	{
		for(int k = 0;k <= n;k ++)
		{
			for(int j = k;j <= maxlen;j ++)
			{
				sum = ((((f[i][j][k]%modnum)*(a[i]%modnum))%modnum)*(ti[j]%modnum))%modnum;
				ans += (sum%modnum);
				ans %= modnum;
			}
		}
	}
	
	printf("%lld",ans);
	return 0;
}

10.31

LCS与LIS
最长公共子序列 与 最长上升子序列
P1439 【模板】最长公共子序列
P1233 木棍加工

bell041030
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
吉林省吉林市第一中学校2015-2016学年高一语文9月检测试题
08-07
【标题】和【描述】中提到的是一个高中语文的9月检测试题，这表明我们需要关注的是高中语文教育的相关知识点，特别是针对高一年级的。【标签】为“课件”，暗示了这可能是一份教学资源，包含了一些教学内容。部分...
医疗废物管理制度2018年4月.doc
10-08
9. 专用运输工具和暂存处清洁消毒制度：规定了运输工具和暂存场所的清洁和消毒频率。 10. 医疗废物管理培训制度：对所有相关人员进行培训，提高其处理能力。 11. 应急预案：为应对医疗废物泄漏、扩散等紧急情况制定...
2020中兴捧月算法大赛迪杰斯特拉赛道初赛题解
gzxb1995的博客
06-27 2561
目录摘要1 程序中使用的数据结构1.1 几个基本数据类型1.2 车道(Lane)1.3 道路(Road)1.4 站点(Station)1.5 货物(Goods)1.6 系统资源(SystemResource)1.7 物流系统(LogisticsSystem)2 算法思路2.1 初赛初版：路由表、深度优先搜索、路径惩罚2.1.1 搜索策略2.1.2 路径惩罚策略2.1.3 货物发送策略2.1.4 效果与问题2.1.5 转向新的思路2.2 初赛终版：多条路径的宽度优先搜索、货物组包、路径统合2.2.1 搜索策略
我的人生之路记录（更新2020年9月）
睶
04-14 1539
我的个人信条 生活低姿态 保持终身学习态度 也许5年，10年，15年后，甚至更远，回头一看少年已变。 我的校园招聘 19届应届生，不才，跨专业，普通一本，哈哈，在招聘中真的是茫茫人海中不起眼的一个。自学Java已有三年，一个人耐住了也习惯了独行。毕竟不能强求别人一起面对困难。走错了很多路，但依然坚挺。 小小成果： 华资软件Java开发，SunLine的Java开发........
2020 HW记录-2
qq_22807425的博客
09-11 1150
2020 HW记录-2 0x00 前言 9.11 hw开始了，记录下整个hw期间爆的0day虽然自己有点想法 有所收获，碍于信息比较敏感的缘故，还是hw结束后再发这儿比较好。hw可是个学习的好机会啊，冲鸭大伙儿 0x02 漏洞 一.常规漏洞 1.Apache DolphinScheduler高危漏洞 漏洞编号：CVE-2020-11974、CVE-2020-13922暂无poc 1.1Apache DolphinScheduler权限覆盖漏洞（CVE-2020-13922） 受影响版本 Apac
绿卡日记：2020-10-29
林兴陆
03-21 102
2020年10月29日，2020年11月绿卡排期公布：EB-1A 杰出人才表A前进183天至2018年12月1日，表B依然是2020年9月1日暂无变化，依然开放表B，允许在美国境内的人提交...
Win10 系统和C盘的 新优化记录（2017年12月02日）
hanshou101的博客
12-02 3418
1.清理C盘空间 可选工具： 1.TweakCube（功能强大，但对于程序员的电脑有一定危险性。清理垃圾文件时，容易清除掉有用文件。） 2.CCleaner（大部分CCleaner的功能，Tuneup都有了。） 3.Tuneup（升级为了AVG Tuneup）（实时在后台保证速度。名副其实的瑞士军刀。和Windows Defender是绝配套餐。） 2.设立内存盘 Z...
医疗废物管理制度(2022年4月)定义.pdf
01-20
9. **专用运输工具和暂存处清洁消毒制度**：保证运输工具和暂存场所的卫生，减少废物处理过程中的二次污染。 10. **医疗废物管理培训制度**：定期对工作人员进行培训，提高他们的废物处理知识和技能。 11. **医疗...
医疗废物管理制度(2018年4月).doc
11-08
《医疗废物管理制度(2018年4月)》详细阐述了医疗机构在处理医疗废物时应遵循的各项规章制度，旨在确保医疗废物的安全管理，防止环境污染，保障公众健康。以下是该制度的关键内容： 一、管理职责 医疗机构需明确各...
进球数 概率 组合数学
Bell的博客
03-30 4288
进球数 题目描述 中超继续进行，广州恒大和北京国安马上开赛了，比赛持续90分钟，为了分析方便，我们把5分钟作为一个时间片，那么比赛会进行18个时间片。在每一个时间片内，恒大踢进1球的概率百分比都是A%，国安踢进1球的概率都是B%。当比赛结束后，两支球队当中，至少有一支球队的进球数是质数的概率是多少？ 输入格式 第一行，一个整数R，表示有R组测试数据。1 <= R <= 10。 每组测试...
掷骰子 概率 动态规划
Bell的博客
04-05 3476
掷骰子 题目描述 玩家A和B正在玩骰子游戏。 A骰子有6个面，第i个面的点数是sideA[i]。 B骰子有6个面，第i个面的点数是sideB[i]。 玩家A总共掷X次A骰子，每次掷骰子得到的面都是1/6的概率。 玩家B总共掷Y次B骰子，每次掷骰子得到的面都是1/6的概率。 玩家最终的总得分就是每次掷骰子得到的点数的总和。 计算玩家A赢得游戏的概率，即玩家A总得分高于玩家B的总得分的概率。 输入格式...
找玩具 概率 dfs
Bell的博客
04-12 2665
找玩具 题目描述 在游戏开始之前，游戏大师在房间的某些地方隐藏了N个玩具。玩具编号为1到N。您的任务是尽可能多地找到这些玩具。 你没有任何辅助信息就能找到第i个玩具的概率是p[i]％。您每找到一个玩具后，有可能可以得到一些辅助信息，这些辅助信息是告诉您其他某些玩具所在的位置。如果您已经知道玩具的位置，您一定会找到它。 给出二维数组clue[1…N][1…N]，其中clue[i][j]=‘Y’表示若...
天平问题 动态规划
Bell的博客
06-20 1858
天平问题 题目描述 小C为了试验小X，便为物竞的小X出了一道物理相关的题：现在给出n个质量的砝码，问小X能称出多少种质量的物品，可是总有好事者想要破坏，于是乎，n 达到了500，远远超出了小X能够承受的范围，锲而不舍的他决定寻求你们的帮助. 注意：天平有两边，两边均可放砝码。 输入格式 第一行，一个正整数N； 以下N行，每行一个不超过200的正整数，依次表示每个砝码的质量。 输出格式 一行，一个数...
黑白球 概率 动态规划
Bell的博客
03-30 1844
黑白球 题目描述 一个箱子里面有n个黑球m个白球。你每小时都随机从箱子里面抽出两个小球，然后把这两个球都染成黑球，然后再放回去。问需要多少小时才能把所有小球变成黑色小球？输出期望值。 输入格式 多组测试数据。 第一行，一个整数G，表示有G组测试数据。1&lt;=G&lt;=101&lt;=G&lt;=101<=G<=10 每组测试数据格式如下: 一行，两...
扑克牌游戏 概率 贪心
Bell的博客
03-30 1843
扑克牌游戏 题目描述 有2N张扑克牌，每张扑克牌都印有一个整数，整数的范围是[1,2N]，而且所有扑克牌的整数都不相同。 Alice和Bob在玩这个游戏。Alice选出其中的N张扑克牌，Bob拿剩下的N张扑克牌。接着游戏进行N轮。每一轮游戏，Alice出一张牌，Bob也出一张牌，如果Alice出的牌印着的整数比Bob的大，那么Alice可以得1分，否则Alice得0分。 给出数组b[1…N]，对于...
排列组合简单练习
Bell的博客
04-12 1546
照相 题目描述 某个班级有N个学生，现在从中选择r个学生出来，从左往右排成一行照相，有多少种不同的排列？答案模10007. 输入格式 一行，N和r。1 <= N<=1000, 1<=r<=N 输出格式 一个整数。 输入样例 1000 2 输出样例 8307 解题思路 代码 #include<iostream> #include<fstream&g...
状态压缩动态规划部分习题详解
Bell的博客
10-24 1314
状态压缩动态规划部分习题详解 状压DP部分题目详解状态压缩动态规划部分习题详解简介经典子集类问题原子弹最短路与状压DP结合送礼物P3959宝藏旅游经典网格类铺地砖一笔画其他类型单词队伍统计 简介 此处记录了一些比较经典或巧妙的简单状压dp题，粗略介绍了做题的一些思路。 经典子集类问题 原子弹 最近，火星研究人员发现了N个强大的原子。他们互相都不一样。 这些原子具有一些性质。当这两个原子碰撞时，其中一个原子会消失，产生大量的能量。 研究人员知道每两个原子在碰撞时的能释放的能量。 你要写一个程序，让它们碰撞
敲砖块 动态规划
Bell的博客
05-03 746
敲砖块 题目描述 在一个凹槽中放置了N层砖块，最上面的一层有N块砖，从上到下每层依次减少一块砖。每块砖都有一个分值，敲掉这块砖就能得到相应的分值，如图所示。 如果你想敲掉第i层的第j块砖的话，若i=1，你可以直接敲掉它，若i＞1，则你必须先敲掉第i-1层的第j和第j+1块砖。 你现在可以敲掉最多M块砖，求得分最多能有多少。 输入格式 输入文件第一行有两个正整数N和M； 接下来的N行，描述这N层砖...
P1525 [NOIP2010 提高组] 关押罪犯
Bell的博客
10-12 591
P1525 [NOIP2010 提高组] 关押罪犯 并查集 贪心 二分答案 二分图判定
2023年9月计算机二级C语言题库: 数组处理和整数运算
2023年9月计算机二级C语言题库.doc中包含了三道题目，每道题目都包含了具体的题目描述和要求在C语言中实现相应的功能。第一道题目要求编写一个函数，该函数能够找出给定数组中低于平均分的成绩，并将这些分数放在另...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值