蓝桥杯刷题日记
文章目录
DAY1
1、递归实现指数型枚举
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-Lx8KAlYW-1653371059604)(C:\Users\12587\AppData\Roaming\Typora\typora-user-images\image-20220110145030958.png)]
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-5bUBXPIl-1653371059605)(C:\Users\12587\AppData\Roaming\Typora\typora-user-images\image-20220110145332484.png)]
图抄网上的,思路很简单(但是竟然不会做,只能说数据结构全部还给王丽萍了),每一个数选择是一种情况,不选择是一种情况,所以只需要每个数标记选择和不选择进行递归就可以了。
#include <iostream>
#include <algorithm>
using namespace std;
int n;
bool flag[16];
void rec(int v)
{
if(v>n)
{
for(int i=1;i<=n;i++)
if(flag[i]) cout<<i<<" ";
cout<<endl;
return;
}
flag[v]=true;
rec(v+1);
flag[v]=false;
rec(v+1);
}
int main(){
cin>>n;
rec(1);
return 0;
}
2、递归实现组合型枚举
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-BLaf4zSd-1653371059606)(C:\Users\12587\AppData\Roaming\Typora\typora-user-images\image-20220110152643697.png)]
还是递归题,参考上题的思路,我们只需要flag里有m个状态为1的时候也返回就可以了(flag里有m个返回要在v>n之前返回,我真是笨蛋)
然后递归多加1个参数来记录有多少个已经被标记为1了就可以了。
这个是自己写的代码:
#include <iostream>
#include <algorithm>
using namespace std;
int n,m;
bool flag[16];
void rec(int v,int num)
{
if(num==m)
{
for(int i=1;i<=n;i++)
if(flag[i]) cout<<i<<" ";
cout<<endl;
return;
}
if(v>n) return;
flag[v]=true;
rec(v+1,num+1);
flag[v]=false;
rec(v+1,num);
}
int main(){
cin>>n>>m;
rec(1,0);
return 0;
}
别人写的代码,我先看看思路:
#include <bits/stdc++.h>
using namespace std;
const int N = 30;
int n, m;
int path[N];
bool st[N];
void dfs(int u, int start)
{
if (u > m)
{
for (int i = 1; i <= m; i++)
cout << path[i] << " ";
cout << endl;
return;
}
for (int i = start; i <= n; i++)
{
if(!st[i])
{
st[i] = true;
path[u] = i;
dfs(u + 1, i);
st[i] = false;
}
}
}
int main()
{
cin >> n >> m;
dfs(1, 1);
return 0;
}
其实思路差不多,只不过他的第二个参数标记成开始的位置,所以其实也差不了多少,不过确实<bits/stdc++.h>这个头文件还有代码的格式更加工整一些。
3、递归实现排列型枚举
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-C2anXm4v-1653371059606)(C:\Users\12587\AppData\Roaming\Typora\typora-user-images\image-20220111091137678.png)]
王丽萍讲过的题,每一次选取一个数,然后这个数不能再次选择,并且记录到数组里去,到最后将所有的数都选择完了,输出这个数组,并且将状态回溯到之前没选择这个数的状态,选择下一个数。
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-Xgn91pEK-1653371059607)(C:\Users\12587\AppData\Roaming\Typora\typora-user-images\image-20220111091335113.png)]
#include <iostream>
#include <algorithm>
using namespace std;
int n;
bool flag[20];
int record[20];
void rec(int v)
{
if(v>n)
{
for(int i=1;i<=n;i++)
cout<<record[i]<<" ";
cout<<endl;
return;
}
for(int i=1;i<=n;i++)
{
if(!flag[i])
{
flag[i]=1;
record[v]=i;
rec(v+1);
flag[i]=0;
}
}
}
int main(){
cin>>n;
rec(1);
return 0;
}
Day2
1、八皇后问题
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-qGWImP01-1653371059607)(C:\Users\12587\AppData\Roaming\Typora\typora-user-images\image-20220111091515395.png)]
王丽萍绝对讲过原题了,做的时候又尬住了,满脑子想的就是当时好像说要用三个bool数组牺牲空间复杂度来简化时间复杂度,三个数组分别标记列、右上斜线和右下斜线。
还是把回溯法 :先放然后递归然后还原原状态的这个关键点忘记了,只要想起来回溯大概是怎么回溯的实际上很快就做的出来。
#include <iostream>
#include <algorithm>
using namespace std;
int n;
char map[10][10];
bool col[10];
bool rg_down[19];
bool rg_up[19];
void rec(int v)
{
if(v==n)
{
for(int i=1;i<=n;i++)
{
for(int j=1;j<=n;j++)
{
cout<<map[i][j];
}
cout<<endl;
}
cout<<endl;
return;
}
for(int i=1;i<=n;i++)
{
if(!col[i]&&!rg_down[v-i+n+1]&&!rg_up[i+v+1])
{
col[i]=rg_down[v-i+n+1]=rg_up[i+v+1]=true;
map[v+1][i]='Q';
rec(v+1);
col[i]=rg_down[v-i+n+1]=rg_up[i+v+1]=false;
map[v+1][i]='.';
}
}
}
int main(){
cin>>n;
for(int i=0;i<10;i++)
{
for(int j=0;j<10;j++)
map[i][j]='.';
}
rec(0);
return 0;
}
*2、费解的开关
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-9PyElJ6m-1653371059608)(C:\Users\12587\AppData\Roaming\Typora\typora-user-images\image-20220111094102026.png)]
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-JdyeiKK9-1653371059608)(C:\Users\12587\AppData\Roaming\Typora\typora-user-images\image-20220111094139993.png)]
这题其实从昨天开始就开始想了,看了答案一开始还没看懂,当时说"枚举第一行开灯的所有情况",我真的是大大的?,后来去搜了份题解才知道,要使第一行不开的灯打开,必须要让正下方的第二行的灯改变状态,直到最后一行查看是否全为绿色,最后得到最小的次数。
确实还是太菜了,尝试一下自己写写看看有什么问题没有。
很奇怪真的太奇怪了,为什么先把灯关了递归和先开了递归得到的结果不一样,真的是匪夷所思,之后可能再看看才能知道。
思路基本上就是先把第一行的所有情况遍历一遍,然后通过check方法判断这个情况是否可行。
check方法其实就是第一行为关灯的话那么把第二行正下方的灯转变状态,然后上面的灯就开。之后第三行在变化让第二行全1………一直到第4行,把第四行全变完直接查看第五行的状态就行了,每一次都计数,这样的话就可以得到最后的次数。
自己写的代码:
#include <iostream>
#include <algorithm>
#include<cstring>
#include<cstdio>
#define INF 0x3f3f3f3f
using namespace std;
int map[7][7];
int back[7][7];
int dx[5]={0,-1,0,1,0};
int dy[5]={0,0,1,0,-1};
int ans;
int check(int k)
{
memcpy(back,map,sizeof(map));
for(int i=1;i<=4;i++)
for(int j=1;j<=5;j++)
{
if(!back[i][j])
{
k++;
for(int m=0;m<5;m++)
back[i+1+dx[m]][j+dy[m]]=!back[i+1+dx[m]][j+dy[m]];
}
}
for(int i=1;i<=5;i++)
if(!back[5][i]) return INF;
return k;
}
void dfs(int cnt,int k)//枚举第一行所有情况
{
if(cnt>5)
{
ans=min(ans,check(k));
return;
}
for(int i=0;i<5;i++)
map[1+dx[i]][cnt+dy[i]]=!map[1+dx[i]][cnt+dy[i]];
dfs(cnt+1,k+1);//按这个开关
for(int i=0;i<5;i++)
map[1+dx[i]][cnt+dy[i]]=!map[1+dx[i]][cnt+dy[i]];
dfs(cnt+1,k);//不按这个开关
}
int main()
{
int n;
cin>>n;
for(int m=0;m<n;m++)
{
memset(map,0,sizeof(map));
for(int i=1;i<=5;i++)
for(int j=1;j<=5;j++)
scanf("%1d",&map[i][j]);
ans=INF;
dfs(1,0);
if(ans<7)
cout<<ans<<endl;
else
cout<<-1<<endl;
}
return 0;
}
别人的代码:
#include <bits/stdc++.h>
using namespace std;
const int N = 6;
const int INF = 0x3f3f3f3f;
int dx[N] = { -1,0,1,0,0 }, dy[N] = { 0,1,0,-1,0 };
char g[N][N], backup[N][N];
int n;
void turn(int x, int y)//将g[x][y]的状态进行改变
{
for (int i = 0; i < 5; i++)
{
int a = x + dx[i], b = y + dy[i];
if (a < 0 || a >= 5 || b < 0 || b >= 5)continue;
g[a][b] ^= 1;
}
}
int main()
{
cin >> n;
while (n--)
{
for (int i = 0; i < 5; i++)cin >> g[i];//直接读一行用位运算
int res = INF;
for (int op = 0; op < 1 << 5; op++)
{
memcpy(backup, g, sizeof g);
int step = 0;
for (int i = 0; i < 5; i++)//第一行变化
{
if (op >> i & 1)
{
step++;
turn(0, i);
}
}
for (int i = 0; i < 4; i++)//check函数
{
for (int j = 0; j < 5; j++)
{
if (g[i][j] == '0')
{
step++;
turn(i + 1, j);
}
}
}
bool is_dark = false;
for (int j = 0; j < 5; j++)
{
if (g[4][j] == '0')
{
is_dark = true;
break;
}
}
if (!is_dark)res = min(res, step);
memcpy(g, backup, sizeof backup);
}
if (res > 6)res = -1;
cout << res << endl;
}
return 0;
}
3、带分数
这个题应该是比较难的,反正我看答案了
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-dNlEhsz3-1653371059609)(C:\Users\12587\AppData\Roaming\Typora\typora-user-images\image-20220112160017161.png)]
但是仔细一想其实假设我知道了要用递归做这题肯定很简单啊,因为其实假设要用递归做了,那么就只能递归构造数字,这个题有三个数字要构造,那么先递归构造第一个,再第一个基础上递归构造两个,然后最后拿两个算第三个,第三个如果满足要求的话就说明是对的,就算做一种。麻烦的点在于要写两个递归函数。
#include <iostream>
#include <algorithm>
#include<cstring>
using namespace std;
long long n,len,ans;
bool flag[10],back[10];
bool check(long long a,long long c)
{
long long b= (n-a)*c;
if(!b||!a||!c) return false;
memcpy(back,flag,sizeof(flag));
for(;b;b/=10)
{
int tmp=b%10;
if(!tmp||back[tmp]) return false;
back[tmp]=true;
}
for(int i=1;i<=9;i++)
if(!back[i]) return false;
return true;
}
void dfs_c(int u,int a,int c)
{
if(u>=10) return;
if(check(a,c)) ans++;
for(int i=1;i<=9;i++)
{
if(!flag[i])
{
flag[i]=true;
dfs_c(u+1,a,c*10+i);
flag[i]=false;
}
}
}
void dfs_a(int u,int a)
{
if(a>=n)return;
if(u>=10) return;
if(a) dfs_c(u,a,0);
for(int i=1;i<=9;i++)
{
if(!flag[i])
{
flag[i]=true;
dfs_a(u+1,a*10+i);
flag[i]=false;
}
}
}
int main()
{
cin>>n;
len=0;
for(long long tmp=n;tmp;tmp/=10)
len++;
dfs_a(0,0);
cout<<ans<<endl;
return 0;
}
Day3
1、飞行员兄弟
这两天去练车了真的挺累的,没怎么写,明天估计还得忙双创,测核酸,哎,这题看了就不会,想了半天了,最后还是抄的答案,
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-NoYBe9Zg-1653371059609)(C:\Users\12587\AppData\Roaming\Typora\typora-user-images\image-20220114233112578.png)]
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-ymI1IsBf-1653371059609)(C:\Users\12587\AppData\Roaming\Typora\typora-user-images\image-20220114233123512.png)]
一直想着怎么用DFS把这个题搜索出来,想着和费解的开关差不多一样来做,但是确实没有找到递推的关系,也没想到最后答案用暴力的做法直接把答案全部枚举出来2^16次方根本每办法递归,所以就用了很巧的位运算来记录每一个开关需要变化的关联开关,然后全部遍历一遍,得出答案。(因为这个题最大的变化开关数为16次,如果相同的开关变化两次,那么就相当于没有任何变化)
/*
基本抄的
*/
#include <iostream>
#include <algorithm>
#include<cstring>
#include<utility>
#include<vector>
#define N 4
#define INF 100
using namespace std;
typedef pair<int,int> PII;
int change[N][N];
int get(int x,int y)
{
return x*N+y;
}
int main()
{
for(int i=0;i<N;i++){
for(int j=0;j<N;j++)
{
for(int k=0;k<N;k++)
change[i][j]+=(1<<get(i,k))+(1<<get(k,j));
change[i][j]-=1<<get(i,j);
}
}
int state=0;
for(int i=0;i<N;i++)
{
string line;
cin>>line;
for(int j=0;j<N;j++)
{
if(line[j]=='+')
state+=1<<get(i,j);
}
}
vector<PII>path;
for(int i=0;i<1<<16;i++)
{
int now=state;
vector<PII> tmp;
for(int j=0;j<16;j++)
{
if(i>>j&1)
{
int x=j/4;
int y=j%4;
now^=change[x][y];
tmp.push_back({x,y});
}
if(!now&&(path.empty()||path.size()>tmp.size()))
path=tmp;
}
}
cout<<path.size()<<endl;
for(int i=0;i<path.size();i++)
{
cout<<path[i].first+1<<" "<<path[i].second+1<<endl;
}
return 0;
}
2、翻硬币
很晚了QAQ,12:00想不出来就睡了,
//第二天
果然休息好做题也快,这题没什么难度,不知道为什么把它列在搜索里,就遍历一遍字符串比比他们哪里不一样,用滑动窗口一下就做出来了
自己写的代码:
#include<iostream>
#include<algorithm>
#include<cstring>
#include<utility>
#include<vector>
using namespace std;
string a;
string b;
int main()
{
cin>>a>>b;
int st=0,end=0,ans=0;
int len=a.size();
for(;st<len;st++)
{
if(a[st]!=b[st])
{
end=st+1;
while(end<len&&a[end]==b[end])
end++;
ans+=end-st;
st=end;
}
}
cout<<ans<<endl;
return 0;
}
别人的代码也差不多
DAY4(二分和前缀和)
1、数的范围
啊这题,耍赖皮拿upper_bound和lower_bound()做出来了,其实这种升序的题一看就知道应该要拿二分什么之类的做
自己写的依靠库函数的代码:
#include<iostream>
#include<algorithm>
#include<cstring>
#include<utility>
#include<vector>
using namespace std;
int n,m;
vector<int> v;
int main()
{
cin>>n>>m;
for(int i=0;i<n;i++)
{
int tmp;
cin>>tmp;
v.push_back(tmp);
}
for(int i=0;i<m;i++)
{
int tmp;
cin>>tmp;
vector<int>::iterator start=lower_bound(v.begin(),v.end(),tmp);
vector<int>::iterator end=upper_bound(v.begin(),v.end(),tmp);
int x=-1,y=-1;
if(start!=v.end()&&*start==tmp)
x=start-v.begin();
if(x!=-1)
y=end-v.begin()-1;
cout<<x<<" "<<y<<endl;
}
return 0;
}
别人自己写的二分:
#include <bits/stdc++.h>
using namespace std;
const int N = 1e5 + 10;
int a[N];
int n, m;
int main()
{
cin >> n >> m;
for (int i = 0; i < n; i++)cin >> a[i];
while (m--)
{
int x;
cin >> x;
int l = 0, r = n - 1;
while (l < r)//第一次二分确定开端
{
int mid = l + r >> 1;//(除2写的好花哨)
if (a[mid] >= x)r = mid;
else l = mid + 1;
}
if (a[l] != x)
cout << -1 << ' ' << -1 << endl;
else//如果找到的话就第二次二分搜索,
{
cout << l << ' ';
l = 0, r = n - 1;
while(l<r)
{
int mid = l + r + 1 >> 1;
if (a[mid] <= x)l = mid;
else r = mid - 1;
}
cout << r << endl;
}
}
return 0;
}
2、数的三次方根
没什么好说的,范围给出了,所以可以2分搜索是否满足条件
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<utility>
#include<vector>
using namespace std;
double get(double x)
{
return x*x*x;
}
int main()
{
double st,end,n,mid;
cin>>n;
st=-10000.0;
end=10000.0;
while(end-st>=1e-7)
{
mid=(st+end)/2.0;
if(get(mid)>=n)
end=mid;
else
st=mid;
}
printf("%.6lf",st);
return 0;
}
DAY5
1、机器人跳跃问题
总的来说和原来程序设计的轻功水上漂差不多就是二分搜索判断它是不是可以
写的过程中没有注意到一个问题就是,假设我的mid很大,带到check函数里,我每次都要吧能量值翻倍,这样的话mid就会越界,所以当能量大于1e5的时候就直接返回true,因为这时候所有的能量值都比h[i]要大。
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<utility>
#include<vector>
using namespace std;
vector<int> v;
vector<int>::iterator it;
int n,l,r,mid;
bool check(int e)
{
for(it=v.begin();it!=v.end();it++)
{
if(e<*it)
e-=*it-e;
else
e+=e-*it;
if(e<0)
return false;
if(e>=1e5)
return true;
}
return true;
}
int main()
{
cin>>n;
for(int i=0;i<n;i++)
{
int tmp;
cin>>tmp;
v.push_back(tmp);
}
l=1;
r=100000;
while(l<r)
{
mid=(l+r)/2;
if(check(mid))
r=mid;
else
l=mid+1;
}
cout<<l<<endl;
return 0;
}
2、四平方和
这题真的是,本来应该是1/17(DAY5)的题了,结果拖到1/18才搞定
这题理论上不应该是二分的,我一开始的想法肯定是暴力枚举,但是枚举肯定不行,它题目要求是联合主键的排序第一个,我的想法是除去前三个以外最后一个二分排序O(n3logn),他最大的5*106,超了
刷题单的答案给了2分搜索,他的思路我不是很能理解,把所有的c和d枚举,然后因为a、b的组合和c和d组合是对称的,所有枚举出来的c和d也包含了a和b的取值,那么就定好a和b以后2分搜索c和d就感觉也是O(n^3logn)不懂,可能顺序搜a和b比搜c和d快,所以就可以做出来
二分的代码:
#include <bits/stdc++.h>
using namespace std;
const int N = 2500010;
struct sum
{
int s, c, d;
bool operator<(const sum& t)const
{
if (s != t.s)return s < t.s;
if (c != t.c)return c < t.c;
return d < t.d;
}
}sum[N];
int n, m;
int main()
{
cin >> n;
for (int c = 0; c * c <= n; c++)
for (int d = c; c * c + d * d <= n; d++)
sum[m++] = { c * c + d * d,c,d };
sort(sum, sum + m);
for(int a = 0; a * a <= n; a++)
for (int b = 0; a * a + b * b <= n; b++)
{
int t = n - a * a - b * b;
int l = 0, r = m - 1;
while (l < r)
{
int mid = l + r >> 1;
if (sum[mid].s >= t)r = mid;
else l = mid + 1;
}
if (sum[l].s == t)
{
cout << a << ' ' << b << ' ' << sum[l].c << ' ' << sum[l].d << endl;
return 0;
}
}
return 0;
}
网页上的解析就是枚举+哈希:
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-sBW1rbCC-1653371059610)(C:\Users\12587\AppData\Roaming\Typora\typora-user-images\image-20220118162311106.png)]
说实话比二分肯定快,估计空间也差不了太多但是想不到(我觉得绝对是比二分好理解的)
代码:
#include<iostream>
#include<cmath>
using namespace std;
const int N = 1e8 + 10;
int h[N];
int main() {
int n;
cin >> n;
//打表,找出1 - n,所有完全平方数两两之和,如果存在只记第一次出现(题目要求找出字典序小的)
for (int i = 0; i * i * 2<= n; i++) {
for (int j = i; j * j + i * i <= n; j++) {
if (!h[i * i + j * j])
h[i * i + j * j] = i + 1;//防止i = 0时在后面判断查找跳过 i = 0的情况
}
}
//0<= a <= b <= c <= d,可以得出a^2 <= n / 4, a^2 + b^ 2 <= n / 2;
for (int i = 0; i * i * 4 <= n; i++) {
for (int j = i; j * j + i * i <= n / 2; j++) {
int t = n - i * i - j * j;
if (h[t]) {
int c = h[t] - 1;
//防止开根号后因为精度关系,向下取整,例:25 开根号得到4.99999向下取整为4;
int d = (sqrt(t - c * c) + 1e-4);
printf("%d %d %d %d", i, j, c, d);
return 0;
}
}
}
return 0;
}
DAY6
1、分巧克力
这题一看就有思路了,明显的二分搜索,在题目给的范围内判断边长能不能符合条件
就是假设每一条巧克力长a[i]宽b[i],分成k*k的正方形,那么总的个数就是(a[i]/k)**(b[i]/k)的总和,那只要二分搜索k就可以了
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<utility>
#include<vector>
using namespace std;
typedef pair<int,int>PII;
const int N=1e5;
int n,m,l,r;
vector<PII> v;
bool check(int k)
{
int sum=0;
for(int i=0;i<v.size();i++)
sum+=(v[i].first/k)*(v[i].second/k);
if(sum>=m) return true;
return false;
}
int main()
{
cin>>n>>m;
for(int i=0;i<n;i++)
{
int a,b;
cin>>a>>b;
v.push_back({a,b});
}
l=1,r=N;
while(l<r)
{
int mid=(l+r)/2;
if(!check(mid)) r=mid;
else l=mid+1;
}
while(!check(l)) l--;
cout<<l;
return 0;
}
Day7
1、激光炸弹(前缀和)
前缀和就是尺取法
就比如
a[1],a[2],a[3]……………………………………a[n];
那么s[1],s[2]…………………………s[n]就表示从a[i]加到a[i],可以o[1]的算出a[i]加到a[j] (s[j]-s[i]);
二维的也是一样,
这个题就是典型的二维前缀和问题:
代码如下:
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<utility>
#include<vector>
using namespace std;
const int N=5010;
int g[N][N];
int n,r;
int main()
{
cin>>n>>r;
r=min(r,5001);
int s=r,m=r;
for(int i=0;i<n;i++)
{
int x,y,w;
cin>>x>>y>>w;
x++;y++;
s=max(s,x);
m=max(m,y);
g[x][y]+=w;
}
for(int i=1;i<=s;i++)
for(int j=1;j<=m;j++)
g[i][j]+=g[i-1][j]+g[i][j-1]-g[i-1][j-1];
int res=0;
for(int i=r;i<=s;i++)
for(int j=r;j<=m;j++)
res=max(res,g[i][j]-g[i-r][j]-g[i][j-r]+g[i-r][j-r]);
cout<<res;
return 0;
}
2、K倍区间
这道题……真的无语
我看着觉得挺简单的,结果还是没做出来;
附图:
我知道要用前缀和,但是假设除了前缀和其他什么都不用,就可以想想,我要定左边界i和右边界j利用s[j]-s[i-1]%k==0才能判断,那么这样的话i和j都要循环就变成了O(n^2)了,题目要求数据最大1e5,所以不行。
要提高效率只能牺牲空间,那这时候应该想到用哈希。问题是哈希怎么做?
抄答案:
这样的话我只需要再开一个数组cnt[N]即可,cnt[x]表示有多少个s[i]%k==x,这样的话时间复杂度就提高到了
想到这还是会寄,因为cnt[0]起始值是1:
哭了。所以最后还是抄完了答案:
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<utility>
#include<vector>
using namespace std;
typedef long long ll;
const int N=100010;
ll a[N];
int cnt[N];
ll n,k;
int main()
{
cin>>n>>k;
for(int i=1;i<=n;i++)
{
int w;
cin>>w;
a[i]=a[i-1]+w;
}
ll res=0;
cnt[0]++;
for(int i=1;i<=n;i++)
{
res+=cnt[a[i]%k];
cnt[a[i]%k]++;
}
cout<<res;
return 0;
}
3、买不到的数目
蓝桥赛A组原题,冲!
这题满脑子想的都是dp,就硬算……,因为dp[i]表示i能不能满足条件,那么dp[i]=dp[i-n]+dp[i-m];在注意一下特殊的情况然后遍历一遍肯定是能过的,代码如下:
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<utility>
#include<vector>
using namespace std;
const int N=1000010;
int a[N];
int n,m;
int main()
{
cin>>n>>m;
if(n>m)
{
int tmp=n;
n=m;
m=tmp;
}
a[0]=1;
a[n]=1;
int ans=n-1;
for(int i=n+1;i<N;i++)
{
if(i<m)
a[i]=a[i-n];
else
a[i]=a[i-m]|a[i-n];
if(!a[i])
ans=i;
}
cout<<ans;
return 0;
}
但是,额,标答用数学方法直接算出来了,我大无语……
(36条消息) 数论系列 求证:设自然数a,b互质,则不能表示成ax+by(x,y为非负整数)的最大整数是ab-a-b_未来AI-CSDN博客
代码:
#include <bits/stdc++.h>
using namespace std;
int main()
{
int p, q;
cin >> p >> q;
cout << (p - 1) * (q - 1) - 1 << endl;
return 0;
}
DAY8
1、蚂蚁感冒
这道题其实就是看蚂蚁相遇了多少次
1、最简单的情况:第一只向左,没有蚂蚁向右,则没有蚂蚁再感染,同理第一只蚂蚁向右,没有蚂蚁向左,也没有蚂蚁感染
2、复杂的情况就是,如果第一只蚂蚁向左,有蚂蚁在左边向右在左边向右都会感染
所以……应该可以了?
这题在算法数上看到过,不是很难,仔细想想就可以出来,但还是看了答案,哎,分清楚就好了
代码:
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<utility>
#include<vector>
using namespace std;
const int N=51;
int n;
int main()
{
cin>>n;
int f,left=0,right=0;
cin>>f;
for(int i=1;i<n;i++)
{
int x;
cin>>x;
if(x>0&&abs(x)<abs(f)) left++;
else if(x<0&&abs(x)>abs(f)) right++;
}
if((f<0&&left==0)||(f>0&&right==0)) cout<<1;
else cout<<(left+right+1);
return 0;
}
2、饮料换购
我的评价是签到题,不过竟然是蓝桥杯A组的题,开始有自信了
代码:
#include<iostream>
#include<algorithm>
using namespace std;
int main()
{
int n;
cin>>n;
int ans=n,reminder=0;
while(n>=3)
{
ans+=n/3;
reminder=n%3;
n=n/3+reminder;
}
cout<<ans;
return 0;
}
3、 01背包
啊这,基本功不扎实了
试一下先;
还行还记得点,但是还是寄了
很简单啊
dp[i] [j]表示选到第i个物品,最大的价值
那么状态转移方程也不难写
因为对于每件物品都有选和不选两种选择,分别对应两种情况,这时候取最大值就可以了;
dp[i] [j]= max(dp[i-1] [j], dp[i-1] [j-v[i]]+w[i]);(如果就j>=v[i])
dp[i] [j]= dp[i-1] [j] (如果j<v[i],放不下了)老是忘
代码:
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<utility>
#include<vector>
using namespace std;
const int M=1010;
int N,V;
int v[M],w[M];
int dp[M][M];
int main()
{
cin>>N>>V;
for(int i=1;i<=N;i++)
cin>>v[i]>>w[i];
for(int i=1;i<=N;i++)
for(int j=0;j<=V;j++)
if(j < v[i])
dp[i][j] = dp[i - 1][j];
else dp[i][j]=max(dp[i-1][j],dp[i-1][j-v[i]]+w[i]);
cout<<dp[N][V];
return 0;
}
4、摘花生
简单的2维dp,不多说,上代码
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<utility>
#include<vector>
using namespace std;
const int N=101;
int g[N][N];
int dp[N][N];
int main()
{
int t;
cin>>t;
while(t--)
{
int c,r;
cin>>r>>c;
memset(dp,0,sizeof(dp));
for(int i=1;i<=r;i++)
for(int j=1;j<=c;j++)
cin>>g[i][j];
for(int i=1;i<=r;i++)
for(int j=1;j<=c;j++)
dp[i][j]=max(dp[i-1][j],dp[i][j-1])+g[i][j];
cout<<dp[r][c]<<endl;
}
return 0;
}
DAY9
1、错误票据
额,只能说输入流没搞明白,
我反正觉得这题应该就是简单的,只要哈希一下就可以了
但是就是这个输入流用的很烂
读到文件末尾可以用(cin>>x).eof()来处理
下次记得了。
代码:
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<utility>
#include<vector>
using namespace std;
const int N=1e6+10;
int a[N];
int l=N,r=0;
int main()
{
int n;
cin>>n;
int ans1,ans2;
int tmp;
while(!(cin>>tmp).eof())
{
a[tmp]++;
l=min(l,tmp);
r=max(r,tmp);
}
for(int i=l;i<=r;i++)
{
if(a[i]==0) ans1=i;
else if(a[i]==2) ans2=i;
}
cout<<ans1<<" "<<ans2;
return 0;
}
我真的是傻瓜,为什么不while(cin>>tmp)呢?
2、剪格子
这就是省赛的压轴题吗
恐怖如斯
我看看能不能做
首先我一开始想的是bfs,试一下:
发现不对
哎,看了答案就是dfs然后就剪枝,把和不是2倍的剪了,把s超过sum/2的剪了
打了遍代码,发现我的代码风格还是挺规范的(傲娇)
蓝桥杯好像和acwing的题不一样,不管了。
代码:
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<utility>
#include<vector>
#define INF 0x3f3f3f3f
using namespace std;
const int N=12;
int n,m;
int sum;
int ans=INF;
int dx[4]={0,1,0,-1};
int dy[4]={1,0,-1,0};
bool flag[N][N];
int g[N][N];
void dfs(int x,int y,int cnt,int s)
{
if(s==sum/2)
{
ans=min(ans,cnt);
return;
}
if(s>sum/2){
return;
}
for(int i=0;i<4;i++)
{
int nx=x+dx[i];
int ny=y+dy[i];
if(nx<=n&&nx>0&&ny<=m&&ny>0&&!flag[nx][ny])
{
flag[nx][ny]=1;
dfs(nx,ny,cnt+1,s+g[nx][ny]);
flag[nx][ny]=0;
}
}
}
int main()
{
cin>>n>>m;
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
{
cin>>g[i][j];
sum+=g[i][j];
}
if(sum%2==0)
{
flag[0][0]=1;
dfs(1,1,1,g[1][1]);
}
ans=(ans==INF?0:ans);
cout<<ans;
return 0;
}
DAY10
昨天颓了
*1、大臣的旅费(树形dp)
算每个城市之间的最大距离就类似那个数据结构里学的,数据量1e5感觉可以解决,试一下。
我的评价是直接寄,这个题遇到了一个没有见过的知识点:树的直径;
在每条边都有权值的情况下,树上两个最远点的距离就是树的直径
求法就是用树形dp,看上去也比较厉害
有点无语,真的好难
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<utility>
#include<vector>
#define INF 0x3f3f3f3f
using namespace std;
typedef long long ll;
const int N=1e5+10;
ll head[N*2],edge[N*2],nxt[N*2],ver[N*2];
ll dis[N],vis[N];
ll total=0,ans=0;
ll n,p,q,d;
void add(int u,int v,int w)
{
ver[++total]=v;
edge[total]=w;
nxt[total]=head[u];
head[u]=total;
}
void dp(int s)
{
vis[s]=1;
for(int i=head[s];i;i=nxt[i])
{
int y=ver[i];
if(vis[y]) continue;
dp(y);
ans=max(ans,dis[s]+dis[y]+edge[i]);
dis[s]=max(dis[s],dis[y]+edge[i]);
}
}
int main()
{
cin>>n;
for(int i=1;i<n;i++)
{
cin>>p>>q>>d;
add(p,q,d);
add(q,p,d);
}
dp(1);
cout<<ans*10+ans*(ans+1)/2;
return 0;
}
我觉得我的图论可能要重新学起。
因为这个有n-1个操作所以才会想到用树存,树的很好用的存储方式就是
tot代表点的数量(最后一个点的标号)
用四个数组来存,ver是存儿子节点的编号的,edge是存边的权重的,next是存这个点的下一个兄弟节点,head是存这个节点的第一个儿子节点
树形dp就是求树的直径,也就是两个最远距离的点,那么最远距离的两个点就是对于公共祖宗节点,这个祖宗节点的最长链和次长链。用dp遍历一遍树,dp存最大的链,如果ans比dp[x]+dp[y]+edge[i]小就说明dp[x]<dp[y]+edge[i],说明dp【x】是次长链,dp[y]+edge[i]是最长链;然后最后求出ans就行了
注意这个因为ans还要计算,所以拿long long存
痛苦
DAY11
1、最长上升子序列
字符串dp,先试一下
袁野:菜,我说不出话了,真的菜,
这个题应该是做过的啊,不知道为什么就是那么的脆弱,有一些可惜。
n^2的复杂度
代码:
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<utility>
#include<vector>
#define INF 0x3f3f3f3f
using namespace std;
const int N=1010;
int a[N];
int dp[N];
int n;
int ans=0;
int main()
{
cin>>n;
for(int i=0;i<n;i++)
cin>>a[i];
for(int i=0;i<n;i++)
{
dp[i]=1;
for(int j=0;j<i;j++)
if(a[j]<a[i]) dp[i]=max(dp[i],dp[j]+1);
ans=max(dp[i],ans);
}
cout<<ans;
}
2、连号区间数
题目描述半天实际就是枚举每个区间的最大值最小值,很简单
代码:
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<utility>
#include<vector>
#define INF 0x3f3f3f3f
using namespace std;
const int N=1e4+10;
int a[N];
int n;
int ans=0;
int main()
{
cin>>n;
for(int i=1;i<=n;i++)
cin>>a[i];
for(int i=1;i<=n;i++){
int l=INF,r=0;
for(int j=i;j<=n;j++)
{
l=min(l,a[j]);
r=max(r,a[j]);
if(r-l==j-i) ans++;
}
}
cout<<ans;
return 0;
}
*3、地宫取宝
感觉就是二维dp,试一下
抄了答案,我真没想到dp数组超二维,再加上两个坐标记录拿了多少物品和最大的价值
这样的话状态转移方程就可以写了,注意边界条件真的还是挺难的
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<utility>
#include<vector>
#define INF 0x3f3f3f3f
using namespace std;
const int N=55;
const int M=15;
const int rem=1000000007;
int n,m,k;
int g[N][N];
int dp[N][N][M][M];
int main()
{
cin>>n>>m>>k;
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++){
cin>>g[i][j];
g[i][j]++;
}
dp[1][1][0][0]=1;
dp[1][1][1][g[1][1]]=1;
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
for(int c=0;c<=k;c++)
for(int d=0;d<=13;d++)
{
dp[i][j][c][d]=(dp[i][j][c][d]+dp[i-1][j][c][d])%rem;
dp[i][j][c][d]=(dp[i][j][c][d]+dp[i][j-1][c][d])%rem;
if(g[i][j]==d&&c>0)
for(int e=0;e<d;e++)
{
dp[i][j][c][d]=(dp[i][j][c][d]+dp[i-1][j][c-1][e])%rem;
dp[i][j][c][d]=(dp[i][j][c][d]+dp[i][j-1][c-1][e])%rem;
}
}
int res=0;
for(int i=0;i<=13;i++)
res=(res+dp[n][m][k][i]);
cout<<res;
return 0;
}
还要注意边界检查,还有0价值的物品,反正挺离谱的
DAY12
昨天又摆烂了
1、快速幂
这个题就是利用任何一个数都可以表示为2的幂之和,再利用a(2i)=(a(2(i-1)))^2得到
这样的话就可以用log^2b的复杂度得到结果
#include<iostream>
#include<algorithm>
using namespace std;
int a,b,p;
int power(int a, int b,int p)
{
int ans=1%p;
for(;b;b>>=1){
if(b&1) ans=(long long)ans*a%p;
a=(long long)a*a%p;
}
return ans;
}
int main()
{
cin>>a>>b>>p;
int ans=power(a,b,p);
cout<<ans;
return 0;
}
2、64位快速乘法
这个题也是一样的,就是也是利用b可以表示成2的幂次的和,然后a*2^i=2a *2^(i-1)
代码:
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
ll a,b,p;
ll mul(ll a,ll b,ll p)
{
ll ans=0;
for(;b;b>>=1)
{
if(b&1) ans=(ans+a)%p;
a=2*a%p;
}
return ans;
}
int main()
{
cin>>a>>b>>p;
ll ans=mul(a,b,p);
cout<<ans;
return 0;
}
这两题都是《算法竞赛进阶指南》的例题,详细的解析可以看到数上,其实我大概也明白了
3、波动数列
看上去直到是dp因为是组合问题,但是实在是想不明白到底是怎么样去dp
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-n2L1Gw8v-1653371059611)(C:\Users\12587\AppData\Roaming\Typora\typora-user-images\image-20220130104129476.png)]
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-yVrFtnKL-1653371059611)(C:\Users\12587\AppData\Roaming\Typora\typora-user-images\image-20220130104136855.png)]
思路大致如上,抄的代码,主要是组合问题一般是2维dp但是j具体选择什么量是需要好好思考的
学一波,打代码:
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<utility>
#include<vector>
#include<limits.h>
using namespace std;
const int MOD=100000007 ;
const int N=1010;
int n,s,a,b;
int dp[N][N];
int mod(int a,int b)
{
return (a%b+b)%b;
}
int main()
{
dp[0][0]=1;
cin>>n>>s>>a>>b;
for(int i=1;i<n;i++)
{
for(int j=0;j<n;j++)
dp[i][j]=(dp[i-1][mod(j-(n-i)*a,n)]+dp[i-1][mod(j+(n-i)*b,n)])%MOD;
}
cout<<dp[n-1][mod(s,n)];
return 0;
}
DAY13
*1、小朋友排队
一开始想用一些比较……暴力的方法来做的(因为看不出来应该用什么算法)但是很明显就是最小换位次数应该是左边有多少个比他大的数和右边有多少个比他小的数,那么就是求逆序数了。
看了看答案,确实不太可能会,因为别人用的就是树状数组,可以很快的求出逆序数,算出逆序数后就可以用
树状数组和线段树有些像(线段树是啥我也不太知道),它的原理又和前缀和比较像,也就是我的a[i]是我要存的数组,然后把它编程线段树以后
例如c[2]就包括a[1]+a[2];
这是基本概念
然后就出现了几个比较离谱的操作:
1、位运算:
int lowbit(x){return x&(-x);}返回最低位的1代表的数值,
2、单点更新:
void update(int x,int y,int n){
for(int i=x;i<=n;i+=lowbit(i)) //x为更新的位置,y为更新后的数,n为数组最大值
c[i] += y;
}
3、区间查询:
因为就是假设要得到sum(5);
c[4]=a[1]+ a[2] +a[3] +a[4];
c[5]=a[5];
所以就是将每个2进制上为1的位的树状数组项进行叠加;
int getsum(int x){
int ans = 0;
for(int i=x;i;i-=lowbit(i))
ans += c[i];
return ans;
}
树状数组可以很快的算出逆序对个数,也就是数组下标作为这个数,而数组记录每个数的个数,这样就可以算出逆序对的个数,也就是这题的思路
代码如下:
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<utility>
#include<vector>
#include<limits.h>
using namespace std;
const int N=100010;
const int H=1000010;
int c[H],b[H];
int n;
int a[N];
int lowbit(int i)
{
return i&(-i);
}
void add(int x,int y)
{
for(int i=x;i<H;i+=lowbit(i))
c[i]+=y;
}
int query(int x)
{
int ans=0;
for(int i=x;i>0;i-=lowbit(i))
ans+=c[i];
return ans;
}
int main()
{
cin>>n;
for(int i=0;i<n;i++){
cin>>a[i];
add(a[i]+1,1);
b[i]=i-query(a[i]);
b[i]-=query(a[i]+1)-query(a[i])-1;
}
memset(c,0,sizeof(c));
for(int i=n-1;i>=0;i--)
{
add(a[i]+1,1);
b[i]+=query(a[i]);
}
long long ans=0;
for(int i=0;i<n;i++)
{
long long total=(long long) b[i]*(b[i]+1)/2;
ans+=total;
}
cout<<ans;
return 0;
}
2、垒骰子
这个题看上去是dp,每个骰子从上往下看应该是有24种状态的,那么也是一种排列组合问题,先试一下
额,不是dp,这个题考的是快速矩阵乘法,
思路大概是每个色子从上往下看的那个面,每个面对应着4种方法,然后f[6]数组记录着最后一个骰子第i面的所有方法,第a[6] [6]记录着往上叠一个骰子的总共的方法,然后用矩阵的快速幂(刚学过)算矩阵的N次方,最后再乘f向量即可。
代码:
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<utility>
#include<vector>
#include<limits.h>
using namespace std;
const int MOD=1e9+7;
const int N=6;
int f[N];
int a[N][N];
int n,m;
void mul1(int *a,int *b, int c[N][N])
{
int tmp[N]={0};
for(int i=0;i<N;i++)
for(int j=0;j<N;j++)
{
tmp[i]=(tmp[i]+(long long)b[j]*c[j][i])%MOD;
}
memcpy(a,tmp,sizeof(tmp));
}
void mul2(int a[][N],int b[][N],int c[][N])
{
int tmp[N][N]={0};
for(int i=0;i<N;i++)
for(int j=0;j<N;j++)
for(int k=0;k<N;k++)
{
tmp[i][j]=(tmp[i][j]+(long long)b[i][k]*c[k][j])%MOD;
}
memcpy(a,tmp,sizeof(tmp));
}
int main()
{
cin>>n>>m;
for(int i=0;i<N;i++)
for(int j=0;j<N;j++)
a[i][j]=4;
for(int i=0;i<m;i++)
{
int x,y,c,d;
cin>>x>>y;
x--;y--;
c=(x+3)%N;
d=(y+3)%N;
a[x][d]=0;
a[y][c]=0;
}
for(int i=0;i<N;i++) f[i]=4;
n--;
for(;n;n>>=1)
{
if(n&1) mul1(f,f,a);
mul2(a,a,a);
}
int res=0;
for(int i=0;i<N;i++)
res=(res+f[i])%MOD;
cout<<res;
return 0;
}
数组越界真的服了
DAY14
1、生命之树
这个题我学会树形dp有思路了
试一下。
啊这,其实思路已经很对了,但是树形dp的编程过程仍然需要熟悉熟悉;
下次应该救会了,这个题世纪就是看子节点的dp结果和0相比的大小,因为最后只能选择最大的dp结果为答案,
其实就是连通域上的点的最大和。
代码:
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<utility>
#include<vector>
#include<limits.h>
using namespace std;
typedef long long ll;
const int N=100010;
int n,tot=0;
ll ans=-1e6-1;
int ver[N*2],head[N*2],Next[N*2],val[N];
ll dp[N];
bool flag[N];
void add(int x,int y)
{
ver[++tot]=y;
Next[tot]=head[x];
head[x]=tot;
}
void dfs(int x){
flag[x]=1;
dp[x]=val[x];
int i;
for(i=head[x];i;i=Next[i])
{
int y=ver[i];
if(flag[y]) continue;
dfs(y);
dp[x]+=max(0LL,dp[y]);
}
ans=max(ans,dp[x]);
}
int main()
{
cin>>n;
for(int i=1;i<=n;i++)
cin>>val[i];
for(int i=1;i<n;i++)
{
int x,y;
cin>>x>>y;
add(x,y);
add(y,x);
}
dfs(1);
cout<<ans;
return 0;
}
2、密码脱落
这个题用dp,很像王丽萍讲过的最长子串匹配
其实如果有思路的话就是找出字符串和他的逆序字符串的最长子串
然后用总长度减去就好
但是想不到啊…………
代码:
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<utility>
#include<vector>
#include<limits.h>
#include<stack>
using namespace std;
const int N=1010;
char a[N],b[N];
int dp[N][N];
int main()
{
cin>>a+1;
int len=strlen(a+1);
for(int i=1;i<=len;i++)
b[i]=a[len-i+1];
for(int i=1;i<=len;i++)
for(int j=1;j<=len;j++)
{
if(a[i]==b[j])
dp[i][j]=dp[i-1][j-1]+1;
else
dp[i][j]=max(dp[i][j-1],dp[i-1][j]);
}
cout<<len-dp[len][len];
return 0;
}
DAY15(枚举)
大年初一开始学北大的算法课,现在开始做做北大算法课的作业
1、特殊密码锁
啊……枚举题都不会做,我真的服了,这个题因为假设要递归,那么他的2^30肯定超时了,所以只能枚举,并且枚举
代码:
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<utility>
#include<vector>
#include<limits.h>
using namespace std;
const int N=32;
char a[N],v[N],t[N];
int s[N];
int n;
void change(int i)
{
v[i]^=1;
if(i>0) v[i-1]^=1;
v[i+1]^=1;
}
int main()
{
cin>>a>>t;
int n=strlen(a);
for(int k=0;k<2;k++)
{
memcpy(v,a,sizeof(a));
memset(s,0,sizeof(s));
if(k==0)
{
for(int i=0;i<n-1;i++)
{
if(v[i]!=t[i]){
change(i+1);
s[i+1]=1;
}
}
if(v[n-1]==t[n-1]) break;
}
else
{
change(0);
s[0]=1;
for(int i=0;i<n-1;i++)
{
if(v[i]!=t[i]){
change(i+1);
s[i+1]=1;
}
}
if(v[n-1]==t[n-1]) break;
}
}
if(v[n-1]==t[n-1])
{
int sum=0;
for(int i=0;i<n;++i)
{
sum+=s[i];
}
cout<<sum;
}
else cout<<"impossible";
return 0;
}
其实和费解的开关一样,因为只要我枚举2种开头开关的情况,如果第i个开关和目标不一样,那么就改变第i+1个开关,直到第n-1个开关,这时候查看第n个开关的状态判断是否满足条件,这个时候不会超时
2进制的巧用和清晰的思路是解决这类题的关键
2、拨钟问题
这个题也是让我大受震撼,首先,有九种枚举的方式,假设把这九种方式都枚举出来,那么算出来应该是9重循环,然后最后是49次的结果也就是218次方,肯定没超时,绝对是可以做的吧
暴力枚举肯定可以做的
然后我去找答案发现确实又很厉害的方式
网上的答案代码,算的更快:
//拨钟问题
#include"stdafx.h"
#include<iostream>
#include<vector>
using namespace std;
//钟表输入
int clock[3][3]={3, 3 ,0,
2, 2, 2,
2, 1, 2};
int mainRoot[9]={3,3,3,3,3,3,3,3,3};//当前的最短路径
void isOk(const int i1,const int i2 ,const int i3);//如果验证成功 那么返回1 否则返回0
void enumerate()//将前3行的拨钟数目进行枚举 返回枚举方法 循环写在这里面
{
int i1=0;int i2=0;int i3=0;//初始时 前三种方法都使用0次
for(int i=0;i<4;i++)
{
i2=0;
for(int j=0;j<4;j++)
{
i1=0;
for(int l=0;l<4;l++)//I1次数
{
isOk(i1,i2,i3);
i1++;
}
i2++;
}
i3++;
}
return ;//结果存在mainRoot中
}
/*
1 ABDE
2 ABC
3 BCEF
4 ADG
5 BDEFH
6 CFI
7 DEGH
8 GHI
9 EFHI
*/
void isOk(const int i1,const int i2 ,const int i3)//如果验证成功 那么返回1 否则返回0
{
int i4=(4-(clock[0][0]+i1+i2)%4)%4;//确定a
int i5=(4-(clock[0][1]+i1+i2+i3)%4)%4;//确定B
int i6=(4-(clock[0][2]+i2+i3)%4)%4; //确定C
int i7=(4-(clock[1][0]+i1+i4+i5)%4)%4;//确定D
int i8=(4-(clock[2][0]+i4+i7)%4)%4;//确定G
int i9=(4-(clock[1][1]+i1+i3+i5+i7)%4)%4;//确定E
//f h i
if((clock[1][2]+i3+i5+i9)%4==0 && (clock[2][1]+i5+i7+i8+i9)%4==0 && (clock[2][2]+i6+i8+i9)%4==0)//满足条件
{
int resultTemp[9]={i1,i2,i3,i4,i5,i6,i7,i8,i9};
int sum=(mainRoot[1]+mainRoot[2]+mainRoot[3]+mainRoot[4]+mainRoot[5]+mainRoot[6]+mainRoot[7]+mainRoot[8]+mainRoot[0]);
if((i1+i2+i3+i4+i5+i6+i7+i8+i9)<sum)//找到了更优秀的解
{
for(int i=0;i<9 ;i++)
{
mainRoot[i]=resultTemp[i];//更新解决方案
}
return;
}
else//不符合要求 直接return继续寻找
{
return;
}
}
return ;
}
int main()
{
for(int i=0;i<3;i++)
{
for(int j=0;j<3;j++)
{
cout<<clock[i][j]<<" ";
}
cout<<endl;
}
enumerate();//将前3行的拨钟数目进行枚举 返回枚举方法
for(int i=0;i<9;i++)
{
cout<<"step "<<i+1<<" "<<mainRoot[i]<<endl;
}
return 0;
}
DAY16(递归)
1、全排列问题
感觉比较轻松,这种题现在也比较熟练了
代码
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<utility>
#include<vector>
#include<limits.h>
using namespace std;
const int N=27;
string a;
char b[N];
vector<string> vs;
int len;
bool flag[N];
void rec(int u)
{
if(u==len) {
cout<<b<<endl;
return ;
}
for(int i=0;i<len;i++)
{
if(!flag[i])
{
flag[i]=1;
b[u]=a[i];
rec(u+1);
flag[i]=0;
}
}
}
int main()
{
cin>>a;
sort(a.begin(),a.end());
len=a.size();
rec(0);
}
2、2的幂次方表示
很常见的递归题,我用了栈,其实可以不用栈的,但是其实都一样,就是每次得到的数如果比2大那就在递归的打印出来直到比2小结束递归就可以了
我的代码:
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<utility>
#include<vector>
#include<limits.h>
#include<stack>
using namespace std;
int n;
void rec(int m)
{
stack<int> s;
int tmp=0;
while(m>0)
{
int mod=m%2;
if(mod) s.push(tmp);
++tmp;
m/=2;
}
while(s.size())
{
int tmp=s.top();
s.pop();
if(tmp==0) cout<<"2(0)";
else if(tmp==1) cout<<"2";
else
{
cout<<"2(";
rec(tmp);
cout<<")";
}
if(s.size()>0) cout<<"+";
}
}
int main()
{
cin>>n;
rec(n);
return 0;
}
别人的代码思路就是减掉2进制最大位数的1,也很厉害很清楚
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<math.h>
void work(int n)
{
if(n==1)//初始判断条件,如果n为1或2则直接输出
{
printf("2(0)");
return;
}
else if(n==2)
{
printf("2");
return;
}
else
{
int j=1,i=0;//j每次乘2,如果大于了n就分解结束,i为当前次数
do
{
j*=2;
if(j>n)
{
j/=2;
if(i==1)//这步非常重要,确定是否需要继续 2()
printf("2");
else
{
printf("2(");
work(i);
printf(")");
}
if(n-j!=0)//如果n分解之后还有剩余的数,那么继续分解
{
printf("+");
work(n-j);
}
return;
}
else
i++;
}while(1);
}
}
int main()
{
int n;
scanf("%d",&n);
work(n);
}
Day17(递归)
1、表达式问题
OpenJudge - 005:Boolean Expressions
服了,这一个递归题竟然在IO上弄了半天
哎getline(cin,string)我真的是用一次忘一次
问题的思路很简单,就是分成两种,一种是表达式,表达式通过|和&连接,需要事先用项去算,另一种是项,项又分三种,一种是()里包含着表达式,通过表达式的值去算,第二种就是通过!直接获得取反的值,第三种就是直接获得该值。
思路不难,代码如下:
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<utility>
#include<vector>
#include<limits.h>
#include<sstream>
#include<stack>
using namespace std;
int n;//记录有多少个表达式
bool factor_value();
bool expression_value();
string expression;
int i;
void trim(string &s)
{
int index = 0;
if(!s.empty())
{
while( (index = s.find(' ',index)) != string::npos)
{
s.erase(index,1);
}
}
}
int main()
{
while(getline(cin,expression))
{
n++;
i=0;
trim(expression);
char out=expression_value()?'V':'F';
cout<<"Expression "<<n<<": "<<out<<endl;
}
return 0;
}
bool expression_value()
{
bool result=factor_value();
bool more=true;
while(more)
{
char op=expression[i];
if(op=='|'||op=='&')
{
i++;
bool value=factor_value();
if(op=='|') result|=value;
else result&=value;
}
else more=false;
}
return result;
}
bool factor_value()
{
bool result=0;
char c=expression[i];
if(c=='(')
{
i++;
result=expression_value();
i++;
}
else
{
i++;
if(c=='!') result=!factor_value();
if(c=='V') result=1;
}
return result;
}
2、简单数的划分
这题决定用dp做锻炼一下
我真的无语,这个题简直就是最简单的dp稍微推出来做一下就可以了
我是直接暴力拿数学推出的dp公式,代码如下
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<utility>
#include<vector>
#include<limits.h>
#include<sstream>
#include<stack>
using namespace std;
const int N=51;
int dp[N][N];
int n;
int main()
{
while(cin>>n)
{
memset(dp,0,sizeof(dp));
for(int i=1;i<=n;i++){
dp[i][1]=1;
dp[1][i]=1;
dp[i][0]=1;
}
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
{
if(i>j) dp[i][j]=dp[i-1][j];
else dp[i][j]=dp[i-1][j]+dp[i][j-i];
}
cout<<dp[n][n]<<endl;
}
return 0;
}
DAY18(动态规划1)
1、拦截导弹
这题简单1维线性dp怎么都会做
状态转移方程就是:dp[i]=max{ dp [j] + 1 }(其中0<=j<i , a[J] < a[I])
要注意的是dp[i]其实为1
代码:
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<utility>
#include<vector>
#include<limits.h>
#include<sstream>
#include<stack>
#include<queue>
using namespace std;
const int N=26;
int n;
int a[N],dp[N];
int main()
{
cin>>n;
for(int i=0;i<n;i++){
cin>>a[i];
dp[i]=1;
}
int ans=0;
for(int i=0;i<n;i++){
for(int j=0;j<i;j++)
{
if(a[i]<=a[j])
dp[i]=max(dp[i],dp[j]+1);
}
ans=max(dp[i],ans);
}
cout<<ans;
return 0;
}
2、Zipper
思路来源于网上(其实也不是想不到,真的就是想复杂了):(64条消息) Zipper动态规划_hejnhong的博客-CSDN博客
String 3能由前两个字符串中的字符混合而成,显然,它的子字符串也必然能由前两个字符串混合成。反之,如果String 3的子串不能由String 1和String 2混合而成,那么该问题的答案就是 no。基于此,我们首先从String 3的第一个字符开始判断,依次到前两个,前三个……如果直到最后的全部字符都满足条件则答案是yes。
定义String 3的前 i+j 个字符,能否由String 1的前 i 个字符和String 2 的前 j 个字符组成的问题,为Yes[i][j]。(注意,这里的 i、j 不是下标,本题的下标都从0开始。)共有两种情况:
(1)String 1 的第 i 个字符是String 3 的第 i+j 个字符,也即String 1[i-1]=String 3[i+j-1],此时子问题为String 3 的前 i+j-1个字符能否由String 1的前 i-1 个字符和String 2的前 j 个字符组成,即 Yes[i-1][j];
(2)String 2 的第 j 个字符是String 3 的第 i+j 个字符,也即String 2[j-1]=String 3[i+j-1],此时子问题为String 3 的前 i+j-1个字符能否由String 1的前 i 个字符和String 2的前 j-1 个字符组成,即 Yes[i][j-1]。
代码:
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<utility>
#include<vector>
#include<limits.h>
#include<sstream>
#include<stack>
#include<queue>
using namespace std;
const int N=201;
char a[N],b[N],c[2*N];
bool dp[N][N];
int n;
int main()
{
cin>>n;
for(int m=1;m<=n;m++)
{
cin>>a>>b>>c;
int i=0,j=0,k=0;
int len1=strlen(a),len2=strlen(b);
memset(dp,0,sizeof(dp));
for(int i=1;i<=len1;i++)
dp[i][0]=(a[i-1]==c[i-1]);
for(int j=1;j<len2;j++)
dp[0][j]=(b[j-1]==c[j-1]);
for(int i=1;i<=len1;i++)
for(int j=1;j<=len2;j++)
{
dp[i][j]=((a[i-1]==c[i+j-1]&&dp[i-1][j])||(b[j-1]==c[i+j-1]&&dp[i][j-1]));
}
cout<<"Data set "<<m<<": ";
if(dp[len1][len2])
cout<<"yes"<<endl;
else
cout<<"no"<<endl;
}
return 0;
}
DAY19
1、最佳加法表达式(高精,动归)
这个题真的也挺难的,最后debug弄了好久,但是思路就是老师给出的那样,
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-l5ghrVUP-1653371059612)(C:\Users\12587\AppData\Roaming\Typora\typora-user-images\image-20220211085215466.png)]
但是不太方便的就是要使用高精度计算真的还是挺麻烦的
代码:
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<string>
#include<utility>
#include<vector>
#include<limits.h>
#include<sstream>
#include<stack>
#include<queue>
using namespace std;
const int N=51;
int m;
string op;
string dp[N][N];
void init()
{
for(int i=0;i<N;i++)
for(int j=0;j<N;j++)
dp[i][j]="";
}
string mins(string a,string b)
{
if(a.empty()) return b;
if(b.empty()) return a;
if(a.length()!=b.length())
return a.length()<b.length()?a:b;
int len=a.length();
for(int i=0;i<len;i++)
if(a[i]!=b[i])
return a[i]<b[i]?a:b;
return a;
}
string add(string a,string b)
{
//cout<<a<<endl<<b<<endl;
string res;
reverse(a.begin(),a.end());
reverse(b.begin(),b.end());
int i=0, cnt=0,l1=a.length(),l2=b.length();
for(i=0;i<l1&&i<l2;i++)
{
char tmp=a[i]+b[i]-'0'+cnt;
if(tmp>'9'){
tmp-=10;
cnt=1;
}
else cnt=0;
res.push_back(tmp);
}
for(;i<l1;i++)
{
char tmp=a[i]+cnt;
if(tmp>'9'){
tmp-=10;
cnt=1;
}
else cnt=0;
res.push_back(tmp);
}
for(;i<l2;i++)
{
char tmp=b[i]+cnt;
if(tmp>'9'){
tmp-=10;
cnt=1;
}
else cnt=0;
res.push_back(tmp);
}
if(cnt) res.push_back('1');
reverse(res.begin(),res.end());
return res;
}
int main()
{
while(cin>>m>>op)
{
init();
int len=op.length();
op=" "+op;
for(int i=0;i<=len;i++){
dp[0][i]=op.substr(1,i);
//cout<<dp[0][i]<<endl;
}
for(int i=1;i<=m;i++)
for(int j=i+1;j<=len;j++)
for(int k=i;k<j;k++)
{
dp[i][j]=mins(dp[i][j],add(dp[i-1][k],op.substr(k+1,j-k)));
//cout<<dp[i][j]<<endl;
}
cout<<dp[m][len]<<endl;
}
return 0;
}
2、复杂的整数划分问题
这个题把动态规划的例题都融合在一起了算是
第1个方式是分盘子问题,每个状态需要考虑有1的情况和没1的情况
若j个数中不包含1,则种数为dp[i-j] [j]
若j个数中包含1,则种数为dp[i-1] [j-1]
因此有dp[i] [j]=dp[i-j] [j]+dp[i-1] [j-1]
第2种方式就是0-1背包不分析了
第3种方式就是普通的完全背包也不分析了
代码如下:
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<string>
#include<utility>
#include<vector>
#include<limits.h>
#include<sstream>
#include<stack>
#include<queue>
using namespace std;
const int N=55;
int n,k;
int dp1[N][N];
int dp2[N][N];
int dp3[N][N];
int main()
{
while(cin>>n>>k)
{
memset(dp1,0,sizeof(dp1));
memset(dp2,0,sizeof(dp2));
memset(dp3,0,sizeof(dp3));
dp1[1][1]=1;
for(int i=2;i<=n;i++)
{
for(int j=1;j<=k;j++)
{
dp1[i][j]=dp1[i-j][j]+ dp1[i-1][j-1];
}
}
cout<<dp1[n][k]<<endl;
for(int i=1;i<=n;i++) dp2[i][0]=1;
dp2[1][1]=1;
for(int i=2;i<=n;i++)
{
for(int j=1;j<=n;j++)
{
if(j>=i) dp2[i][j]=dp2[i-1][j]+dp2[i-1][j-i];
else dp2[i][j]=dp2[j][j];
}
}
cout<<dp2[n][n]<<endl;
dp3[1][1]=1;
for(int i=1;i<=n;i++)
if(i&1)
dp3[i][0]=1;
for(int i=1;i<=n;i++) dp3[1][i]=1;
for(int i=3;i<=n;i+=2)
{
for(int j=1;j<=n;j++)
{
if(j>=i) dp3[i][j]=dp3[i-2][j]+dp3[i][j-i];
else dp3[i][j]=dp3[i-2][j];
}
}
cout<<dp3[n-((n+1)&1)][n]<<endl;
}
return 0;
}
3、Charm Bracelet
OpenJudge - 016:Charm Bracelet
状态压缩的0-1背包,挺简单的
不分析了
代码:
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<string>
#include<utility>
#include<vector>
#include<limits.h>
#include<sstream>
#include<stack>
#include<queue>
using namespace std;
const int N=3403;
const int M=12881;
int dp[2][M];
int n,m;
int w[N];
int v[N];
int main()
{
cin>>n>>m;
for(int i=1;i<=n;i++)
cin>>w[i]>>v[i];
for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++)
{
if(j-w[i]>=0)
dp[i%2][j]=max(dp[(i+1)%2][j],dp[(i+1)%2][j-w[i]]+v[i]);
else dp[i%2][j]=dp[(i+1)%2][j];
}
}
cout<<dp[n%2][m];
return 0;
}
4、分蛋糕
思路参考郭老师的讲法
dp[i] [j] [k]记录着这块蛋糕最大的子蛋糕的面积,i表示蛋糕长度,j表示蛋糕宽度,k表示砍了多少刀,
那么这个时候横着切一刀就变成两块子蛋糕,竖着切1刀也变成两块子蛋糕,问题规模就减小了,
那么选择r作为左边子蛋糕的边长,那么i-r就是右边子蛋糕的边长,记左边子蛋糕切kk刀,那么因为总共切k刀,所以右边要切k-kk-1刀,这样的话两块蛋糕中大的面积为最大子蛋糕,每次只要取最小的就可以了,
竖着切也是同理
代码:
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<string>
#include<utility>
#include<vector>
#include<limits.h>
#include<sstream>
#include<stack>
#include<queue>
#define INF 0x3f3f3f3f
using namespace std;
const int N=21;
int w,h,m;
int dp[N][N][N];
int main()
{
while(cin>>w>>h>>m&&(w||h||m))
{
m--;
memset(dp,INF,sizeof(dp));
for(int a=1;a<=w;a++)
{
for(int b=1;b<=h;b++)
{
for(int c=0;c<=m;c++)
{
if(c==0) dp[a][b][c]=a*b;
else if(a*b<c+1) dp[a][b][c]=INF;
else{
for(int r=1;r<a;r++)
for(int k=0;k<c;k++)
dp[a][b][c]=min(dp[a][b][c],max(dp[r][b][k],dp[a-r][b][c-k-1]));
for(int l=1;l<b;l++)
for(int k=0;k<c;k++)
dp[a][b][c]=min(dp[a][b][c],max(dp[a][l][k],dp[a][b-l][c-k-1]));
}
}
}
}
cout<<dp[w][h][m]<<endl;
}
return 0;
}
新学期开始啦
Day1(最短路和最小生成树问题)
1、3255 – Roadblocks (poj.org)
求次短路参考dijkstra算法,
一个数组记录最短路,一个数组记录次短路,
最短路的计算方式和标准的dijstra算法仍然一致,
次短路其实就是要么是最短路加一条较长边,要么就是次短路加一条短边,比较一下即可
#include<bits/stdc++.h>
#define INF 0x3f3f3f3f
using namespace std;
typedef long long ll;
typedef pair<int,int> PII;
const int N=5050;
struct edge{
int cost;
int to;
edge(int a,int b)
{
to=a;
cost=b;
}
};
int n,r;
vector<edge> g[N];
int d1[N],d2[N];
int main()
{
cin>>n>>r;
for(int i=0;i<r;i++)
{
int a,b,c;
cin>>a>>b>>c;
a--;
b--;
g[a].push_back(edge(b,c));
g[b].push_back(edge(a,c));
}
memset(d1,INF,sizeof(d1));
memset(d2,INF,sizeof(d2));
priority_queue<PII,vector<PII>,greater<PII> > que;
d1[0]=0;
que.push(PII(0,0));
while(que.size())
{
PII p=que.top();
que.pop();
int v=p.second;
int d=p.first;
if(d2[v]<d) continue;
for(int i=0;i<g[v].size();i++)
{
edge& e=g[v][i];
int dist=d+e.cost;
if(d1[e.to]>dist)
{
swap(d1[e.to],dist);
que.push(PII(d1[e.to],e.to));
}
if(d2[e.to]>dist&&d1[e.to]<dist)
{
d2[e.to]=dist;
que.push(PII(d2[e.to],e.to));
}
}
}
cout<<d2[n-1];
return 0;
}
2、3723 – Conscription (poj.org)
思考:最小生成树,这里用的是kruskal算法(poj的runtimeerror太奇怪了,希望下次不要这么阴间了)
#include<iostream>
#include<algorithm>
#include<utility>
#include<vector>
#include<cstring>
#include<cstdio>
#define INF 0x3f3f3f3f
using namespace std;
typedef long long ll;
typedef pair<int,int> PII;
const int N=100010;
const int R=100010;
struct edge{
int from;
int cost;
int to;
bool operator <(const edge& t)const
{
return cost<t.cost;
}
};
int f[N],r[N];
void init(int n)
{
for(int i=0;i<n;i++)
{
f[i]=i;
r[i]=0;
}
}
int find(int x)
{
return x==f[x]?x:find(f[x]);
}
void unite(int x,int y)
{
x=find(x);
y=find(y);
if(r[x]<r[y])
{
f[x]=y;
}
else{
f[y]=x;
if(r[x]==r[y])r[x]++;
}
}
bool same(int x,int y)
{
return find(x)==find(y);
}
int n,m,rd;
edge e[R];
void solve()
{
cin>>n>>m>>rd;
init(n+m);
for(int i=0;i<rd;i++)
{
scanf("%d%d%d",&e[i].from,&e[i].to,&e[i].cost);
e[i].cost=-e[i].cost;
}
sort(e,e+rd);
int res=0;
for(int i=0;i<rd;i++)
{
edge e1=e[i];
if(!same(e1.from,e1.to+n))
{
unite(e1.from,e1.to+n);
res+=e1.cost;
}
}
cout<<(n+m)*10000+res<<endl;
}
int main()
{
int t;
cin>>t;
while(t--)
{
solve();
}
return 0;
}
DAY2(图论day2)
1、3169 – Layout (poj.org)(Bellman-Ford算法)
因为对于这个题目我们需要求从1到N号牛的最大距离,所有的距离都满足:
1、没有关系:d[i]<=d[i+1]
2、有亲密关系:d[AL]+DL>=d[BL];
3、有不好关系:d[AD]+DD<=d[BD];
这个形式统一可以化成求最短路: d [ i ] = m i n ( d [ j ] + e ( j , i ) ) ∣ e ( i , j ) ∈ E ) d[i]=min(d[j]+e(j,i))|e(i,j)\in E) d[i]=min(d[j]+e(j,i))∣e(i,j)∈E)
这个时候我们需要求最大距离d其实也就是求d[N]-d[1];
最后代码如下:
#include<iostream>
#include<algorithm>
#include<utility>
#include<vector>
#include<cstring>
#include<cstdio>
#define INF 0x3f3f3f3f
using namespace std;
typedef long long ll;
typedef pair<int,int> PII;
struct edge{
int from;
int cost;
int to;
};
const int N=20010;
int n,ml,md;
edge eg[2*N];
int d[N];
bool updated;
void bellmen()
{
for(int k=0;k<=n;k++)
{
updated=false;
for(int i=0;i+1<n;i++)
{
if(d[i+1]<INF&&d[i]>d[i+1])
{
d[i]=d[i+1];
updated=true;
}
}
for(int i=0;i<ml;i++)
{
edge e=eg[i];
int x=e.from;
int y=e.to;
int cost=e.cost;
if(d[x]<INF&&d[y]>d[x]+cost){
updated=true;
d[y]=d[x]+cost;
}
}
for(int i=ml;i<ml+md;i++)
{
edge e=eg[i];
int x=e.from;
int y=e.to;
int cost=e.cost;
if(d[y]<INF&&d[x]>d[y]+cost){
updated=true;
d[x]=d[y]+cost;
}
}
}
}
int main()
{
cin>>n>>ml>>md;
for(int i=0;i<ml;i++)
{
int al,bl,len;
cin>>al>>bl>>len;
eg[i].from=al-1;
eg[i].to=bl-1;
eg[i].cost=len;
}
for(int i=ml;i<ml+md;i++){
int al,bl,len;
cin>>al>>bl>>len;
eg[i].from=al-1;
eg[i].to=bl-1;
eg[i].cost=-len;
}
bellmen();
if(updated){
cout<<-1<<endl;
return 0;
}
memset(d,INF,sizeof(d));
d[0]=0;
bellmen();
int res=d[n-1];
if(res==INF) res=-2;
cout<<res<<endl;
return 0;
}
DAY3
1、P1119 灾后重建 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)(Floyd算法)
这个题最重要的就是他的点非常少,但是边很多。
Bellmen算法:O(VE)复杂度较高,但可以作用于有负边的图
Dijstra算法:O(Elog(V))复杂度相对均衡高效的算法,但是不能作用于有负边的图
Floyd算法:O( V 3 V^3 V3)复杂度较高,但是实现简单,在密集矩阵的情况下用是非常有效的
确定使用Floyd算法之后就非常顺利了,因为重建时间是递增的,所以每一次只要按照时间戳来进行最短距离的更新就可以很快的解决每一次询问的问题;
需要注意的是两种特殊情况下输出-1,第一种是还没到时间,目标城镇还没有重建,第二种是目标城镇距离无限大(d==INF)
#include<bits/stdc++.h>=
#define INF 0x3f3f3f3f
using namespace std;
const int N=205;
const int M=20000;
typedef long long ll;
typedef pair<int,int> PII;
int n,m,q;
int edge[N][N];
int t[N];
void update(int k)
{
for(int i=0;i<n;i++)
{
for(int j=0;j<n;j++)
{
if(edge[i][j]>edge[i][k]+edge[k][j])
{
edge[i][j]=edge[j][i]=edge[i][k]+edge[k][j];
}
}
}
}
int main()
{
cin>>n>>m;
for(int i=0;i<n;i++)
cin>>t[i];
for(int i=0;i<n;i++)
for(int j=0;j<n;j++)
edge[i][j]=INF;
for(int i=0;i<m;i++)
{
int a,b,c;
cin>>a>>b>>c;
edge[a][b]=c;
edge[b][a]=c;
}
cin>>q;
int now=0;
for(int i=0;i<q;i++)
{
int a,b,c;
cin>>a>>b>>c;
while(t[now]<=c&&now<n)
{
update(now);
now++;
}
if(t[a]>c||t[b]>c)
cout<<-1;
else{
if(edge[a][b]==INF) cout<<-1;
else cout<<edge[a][b];
}
cout<<endl;
}
}
DAY4
1、(23条消息) 蓝桥杯 RSA解密_林十六要努力的博客-CSDN博客_蓝桥杯rsa解密(扩展欧几里得+快速幂)
思路博客里有了,代码:
#include<bits/stdc++.h>=
#define INF 0x3f3f3f3f
using namespace std;
const int N=205;
const int M=20000;
typedef long long ll;
ll exctgcd(ll a,ll b,ll &x,ll &y)
{
ll d=a;
if(b!=0)
{
d=exctgcd(b,a%b,y,x);
y-=(a/b)*x;
}
else{
x=1;y=0;
}
return d;
}
ll quick_multi(ll a,ll b,ll n)
{
ll res=0;
for(;b;b>>=1)
{
if(b&1)
res=(res+a)%n;
a=(a+a)%n;
}
return res;
}
ll quick_power(ll a,ll b,ll n)
{
ll res=1;
for(;b;b>>=1)
{
if(b&1)
res=quick_multi(res,a,n);
a=quick_multi(a,a,n);
}
return res;
}
int main()
{
/*
ll n= 1001733993063167141;
for(ll i=2;i*i<n;i++)
{
if(n%i==0)
{
cout<<i<<" "<<n/i<<endl;
}
}
*/
/*
ll p=891234941;
ll q=1123984201;
ll a=(p-1)*(q-1);
ll d= 212353;
ll x=0,y=0;
ll res=exctgcd(a,d,x,y);
cout<<res<<" "<<x<<" "<<y<<endl;
y=(y+a)%a;
cout<<y;
*/
ll p=891234941;
ll q=1123984201;
ll e=823816093931522017;
ll n=1001733993063167141;
ll c=20190324;
ll res=quick_power(c,e,n);
cout<<res%n;
return 0;
}
DAY5
(昨天摆了)
1、[P1551 - 蓝桥杯2021初赛] 直线 - New Online Judge (ecustacm.cn)
这个题是一个模拟题,真不知道怎么简化,大概就是想着把所有直线表示成 a x + b y = c ax+by=c ax+by=c的样式,这样的话每一次将 ( a , b , c ) ∈ S (a,b,c)\in S (a,b,c)∈S求出S的size就可以了。
但是打代码的速度好慢,修bug修的也好慢。
(下次或许想清楚了再写,符号搞来搞去真的太慢了)
源代码:
#include<iostream>
#include<vector>
#include<algorithm>
#include<cstring>
#include<queue>
#define INF 0x3f3f3f3f
using namespace std;
const int N=205;
const int M=20000;
typedef long long ll;
typedef pair<int,int> P;
struct edge{
int to;
int cost;
edge(int a,int b){
to=a;
cost=b;
}
};
vector<edge> g[2022];
int gcd(int a,int b)
{
return b?gcd(b,a%b):a;
}
void init()
{
for(int i=1;i<=2021;i++)
{
for(int j=1;j<=21&&j<=2021-i;j++)
{
int a=i,b=i+j;
int c=gcd(a,b);
c=a*b/c;
g[a].push_back(edge(b,c));
g[b].push_back(edge(a,c));
}
}
}
int d[2022];
void dijstra()
{
priority_queue<P,vector<P>,greater<P> > q;
memset(d,INF,sizeof(d));
d[1]=0;
q.push(P(0,1));
while(q.size())
{
P p=q.top();
q.pop();
int f=p.second;
int c=p.first;
//cout<<f<<" "<<c<<endl;
if(d[f]<c) continue;
for(int i=0;i<g[f].size();i++)
{
edge ed=g[f][i];
if(d[ed.to]>d[f]+ed.cost)
{
d[ed.to]=d[f]+ed.cost;
q.push(P(d[ed.to],ed.to));
}
}
}
return;
}
int main()
{
init();
dijstra();
cout<<d[2021]<<endl;
return 0;
}
2、[P1553 - 蓝桥杯2021初赛] 路径 - New Online Judge (ecustacm.cn)
这个题就是标准的Dijstra算法,但是我模板背不好,初始化还搞错
真的服了,气死了,去写数据库了
#include<iostream>
#include<vector>
#include<algorithm>
#include<cstring>
#include<queue>
#define INF 0x3f3f3f3f
using namespace std;
const int N=205;
const int M=20000;
typedef long long ll;
typedef pair<int,int> P;
struct edge{
int to;
int cost;
edge(int a,int b){
to=a;
cost=b;
}
};
vector<edge> g[2022];
int gcd(int a,int b)
{
return b?gcd(b,a%b):a;
}
void init()
{
for(int i=1;i<=2021;i++)
{
for(int j=1;j<=21&&j<=2021-i;j++)
{
int a=i,b=i+j;
int c=gcd(a,b);
c=a*b/c;
g[a].push_back(edge(b,c));
g[b].push_back(edge(a,c));
}
}
}
int d[2022];
void dijstra()
{
priority_queue<P,vector<P>,greater<P> > q;
memset(d,INF,sizeof(d));
d[1]=0;
q.push(P(0,1));
while(q.size())
{
P p=q.top();
q.pop();
int f=p.second;
int c=p.first;
//cout<<f<<" "<<c<<endl;
if(d[f]<c) continue;
for(int i=0;i<g[f].size();i++)
{
edge ed=g[f][i];
if(d[ed.to]>d[f]+ed.cost)
{
d[ed.to]=d[f]+ed.cost;
q.push(P(d[ed.to],ed.to));
}
}
}
return;
}
int main()
{
init();
dijstra();
cout<<d[2021]<<endl;
return 0;
}
DAY6
*1、[P1559 - 蓝桥杯2021初赛] 异或数列 - New Online Judge (ecustacm.cn)
第一题直接抄答案了(哭哭)
看的出来是dp但是实在是太弱了,但实际上用递归已经把状态转移方程写出来了,只是想的太慢了;
dp[i] [j]表示第i轮转换能否称取j重量的物品,最后再全部取和就可以算出最后的结果、
代码:
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<set>
using namespace std;
const int N=110;
const int M=100010;
int n,sum,ans;
int a[N];
int dp[N][M];
int main()
{
cin>>n;
for(int i=1;i<=n;i++){
cin>>a[i];
sum+=a[i];
}
dp[1][a[1]]=1;
for(int i=2;i<=n;i++)
{
int x=a[i];
for(int j=1;j<=sum;j++)
{
dp[i][j]=dp[i-1][j];
}
dp[i][x]=1;
for(int j=1;j<=sum;j++)
{
if(dp[i-1][j]==1)
{
dp[i][j+x]=1;
dp[i][abs(j-x)]=1;
}
}
}
for(int i=1;i<=sum;i++)
{
if(dp[n][i]==1)
ans++;
}
cout<<ans;
return 0;
}
蓝桥杯结束了
因为疫情,事情好多,做什么事情都断断续续的,虽说失败不应该找借口,但是不得不承认,这次蓝桥杯结果确实是差强人意,我觉得算法还是非常重要的,明年应该巩固算法的知识,等到回到家里面应该是有时间对这次蓝桥杯好好复盘的。
关于CSDN博客
虽然比较失败,但是我觉得我的笔记里应该还是包含很多经典的题,希望能帮到大家。也希望明年蓝桥杯大家能取得好成绩吧!