[雅礼NOIP2018集训 day3]
考试的时候刚了T1两个小时线段树写了三个子任务结果发现看错了题目,于是接下来一个半小时我自闭了
result=历史新低
这告诉我们,打暴力要端正态度,尤其是在发现自己之前出锅的情况下要保持心态的平和,不能和今天的比赛一样后面差不多直接弃疗
T1:
题意就是我们要做多次倒三角的区间加,最后统计全部的异或和。不幸的是当我看到空间限制512MB的时候就直接暴力上线段树了,凉心出题人
正解是很巧妙的二维前缀和做法
考虑我们暴力怎么做--对倒三角的每一行差分,最后统计一次,这样的复杂度是$O(nq)$的
这个时候可以发现每一次倒三角我们改变的差分序列是可以二维差分优化的。其实就是对三角的竖着的直角边和那条斜边在维护差分数组,最后再统计答案就好了
其他的做法不会啊,果然还是太弱
#include<algorithm>
#include<cstdio>
#include<iostream>
#include<cstring>
#include<time.h>
using namespace std;
typedef long long ll; const int N=3e3+;
ll n,q;
ll a[N][N],b[N][N],c[N][N],d[N][N],e[N][N];
inline ll read()
{
char ch=getchar();
ll s=,f=;
while (ch<''||ch>'') {if (ch=='-') f=-;ch=getchar();}
while (ch>=''&&ch<='') {s=(s<<)+(s<<)+ch-'';ch=getchar();}
return s*f;
}
int main()
{
freopen("u.in","r",stdin);
freopen("u.out","w",stdout);
//double st=clock();
n=read();q=read();
while (q--)
{
ll r=read(),c=read(),l=read(),s=read();
a[r][c]+=s;a[r+l][c]-=s;
b[r-c+n-][r]-=s;b[r-c+n-][r+l]+=s;
}
for (int i=;i<=n;i++)
for (int j=;j<=n;j++)
{
c[i][j]=c[i-][j]+a[i][j];
d[i][j]=d[i-][j-]+b[i-j+n][i];
e[i][j]=e[i][j-]+c[i][j]+d[i][j];
}
/*for (int i=1;i<=n;i++)
{
for (int j=1;j<=n;j++)
printf("%d ",c[i][j]+d[i][j]);
printf("\n");
}*/
//printf("\n");
ll ans=;
for (int i=;i<=n;i++)
for (int j=;j<=n;j++) ans^=e[i][j];
/*for (int i=1;i<=n;i++)
{
for (int j=1;j<=n;j++) printf("%d ",e[i][j]);
printf("\n");
}*/
printf("%lld\n",ans);
//double ed=clock();
//printf("%lf\n",ed-st);
}
T2:
数据范围好像就是给你状压DP的,状态分别是当前还剩下的球(一个0/1序列)和剩下球的个数(这里有个坑点,不能直接通过0/1序列记忆化,因为最高位可能是0,这样可能出现长度不同但是0/1序列相同的两个状态)
我们数组记忆化肯定是不行的,那么大的就只好开map了。鉴于上述的坑点,看代码注释了解如何避免吧
所谓最优策略,其实就是状态转移的时候取max就好
还有就是关于那个erase函数,删掉一个球,相当于在一个二进制数里面去掉一位。设去掉的二进制位为k,笔者的思路就是把0-k-1先取出来(预处理111...这样的数字&一下就好了),再把原来的数后面变成0(注意0的个数要比原来的位数少1,因为有一位被去掉了)
然后这题好像就A了,考场上像个傻叉一样的写搜索,关键是之前的1分暴力写完了没打return 0,输出两个答案的我彻底凉凉了,22分暴力都没有拿到
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<iostream>
#include<map>
using namespace std; const int N=;
int n,k;
int pre[N];
namespace calc
{
const int M=;
double a[<<M+];
map <int,double> p[N];
void init()
{
for (int i=;i<<<M+;i++) a[i]=-;
}
bool count(int bit,int len)
{
if (len<=M) return a[<<len|bit]!=-;//坑点处理在这儿,看到那个|没有?std太伟大了
else return p[len].count(bit);
}
double &find(int bit,int len)
{
if (len<=M) return a[<<len|bit];
else return p[len][bit];
}
}
int erase(int bit,int i)
{
return (bit&pre[i-])|((bit>>(i+))<<i);
}
double dfs(int bit,int len)
{
if (len<=k) return ;
if (calc::count(bit,len)) return calc::find(bit,len);//记忆化
double &res=calc::find(bit,len);
res=;
for (int i=,j=len-;i<=j;++i,--j)
{
if (i<j) res+=max(dfs(erase(bit,i),len-)+(bit>>i&),dfs(erase(bit,j),len-)+(bit>>j&))*;//*2是因为正的第i个,反的第len-i+1个也是这个位置
else res+=dfs(erase(bit,i),len-)+(bit>>i&);
}
return res/=len;//别忘了/len
}
int main()
{
//freopen("v.in","r",stdin);
//freopen("v.ans","w",stdout);
scanf("%d%d",&n,&k);
char ch[N];
scanf("%s",ch);
calc::init();
pre[]=;
for (int i=;i<N;i++) pre[i]=pre[i-]|(<<i);
int bit=;
k=n-k;
for (int i=;i<n;i++)
{
bit|=(ch[i]=='W')<<i;
}
printf("%.10lf\n",dfs(bit,n));
return ;
}
T3:
有两个需要明确的性质
1.不要变的边我们不变,因为变了我们还要至少花费一次代价把它变回来
2.一棵树里变了的边的条数就是在这些边覆盖的点集中奇数度数点的个数/2(度数是翻转的边带来的)。这个好像挺显然的,因为一条翻转的边只会给两个端点带来奇数的度数,中间的都是偶数的度数
考虑dp,dp[x][0/1]表示节点x与父亲的边是否翻转,最少的奇数点的个数和翻转的总长度(我们的dp数组存的这样的结构体)
发现目标是在最少的奇数点的基础上最小化翻转的总长度,也就是前者优先,于是我们用pair
考虑如何转移
设x与父亲的边的类型是type(if (d==2) type=2 else type=c^d)
如果type==1||type==2的话,dp[x][0]显然可以设成inf了,我们只能转移dp[x][1]。设tmp0表示仅仅考虑x的子树,当前点x不是路径端点的最小代价(注意x可能是被一条路径经过);tmp1表示仅仅考虑x的子树,当前点是路径端点的最小代价。这两个如何计算参见代码。dp[x][1]=max((tmp0.fi+1,tmp0.se+1),(tmp1.fi,tmp1.se+1))
还有就是type==0||type==2,这种情况也差不多,不理解参考代码吧
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<vector>
#define pii pair<int,int>
using namespace std; const int N=1e5+;
const int inf=1e9+;
int n;
vector <pii> mp[N];
pii dp[N][];
inline int read()
{
char ch=getchar();
int s=,f=;
while (ch<''||ch>'') {if (ch=='-') f=-;ch=getchar();}
while (ch>=''&&ch<='') {s=(s<<)+(s<<)+ch-'';ch=getchar();}
return s*f;
}
pii operator + (pii a,pii b){return make_pair(a.first+b.first,a.second+b.second);}
void dfs(int x,int fa,int type)
{
pii tmp0(,),tmp1(inf,inf);//tmp0是不以它为断电答案,tmp1是以它为端点的最小答案
for (int i=;i<mp[x].size();i++)
{
int y=mp[x][i].first;
if (y==fa) continue;
dfs(y,x,mp[x][i].second);
pii nxt0,nxt1;
nxt0=min(tmp0+dp[y][],tmp1+dp[y][]);
nxt1=min(tmp0+dp[y][],tmp1+dp[y][]);
tmp0=nxt0;tmp1=nxt1;
}
if (type==||type==)
{
dp[x][]=min(tmp0,make_pair(tmp1.first+,tmp1.second));
}
else dp[x][]=make_pair(inf,inf);
if (type==||type==)
{
dp[x][]=min(make_pair(tmp1.first,tmp1.second+),make_pair(tmp0.first+,tmp0.second+));
}
else dp[x][]=make_pair(inf,inf);
}
int main()
{
n=read();
for (int i=;i<n;i++)
{
int a=read(),b=read(),c=read(),d=read();
if (d!=) d=c^d;
mp[a].push_back(make_pair(b,d));
mp[b].push_back(make_pair(a,d));
}
dfs(,-,);
printf("%d %d\n",dp[][].first/,dp[][].second);
return ;
}
最新文章
- (翻译)FIFO In Hardware
- java基础复习 - 自动装箱
- HTTP,TCP/IP协议
- Redis入门(一)系统安装
- Tracert 转
- LightOJ1125 Divisible Group Sums(DP)
- oracle启动关闭命令
- form表单验证2
- SpriteKitCommonUse
- WeX5是主要进行app开发吗?能开发微信App吗?
- 【C语言】01-函数
- Delphi实现无标题有边框的窗体
- SQL排序 空值的后面
- SVN 代码下载,上传
- 微信小程序之----弹框组件modal
- 我的Spring学习记录(四)
- pandas groupby生成新的dataframe
- The iOS Simulator deployment target is set to 6.0
- CI框架 -- 密码哈希
- Spring Boot 基础配置
热门文章
- win32 Service memory leak
- 均匀分布(uniform distribution)期望的最大似然估计(maximum likelihood estimation)
- BZOJ 3674/BZOJ 3673 主席树
- 什么是CNN--Convolutional Neural Networks
- php 生成不重复的随机字符串
- vue-cli webpack配置中 如何启动less-loader sass-loader
- Book---强连通分量
- RocketMQ学习笔记(10)----RocketMQ的Producer 事务消息使用
- Unity 设置2台摄像机的叠加
- win10 MongoDB安装