T1:math

题目链接:

http://zhengruioi.com/contest/156/problem/471

题解:

先讲讲我的乱搞做法。对于前面70%,我跑了背包。因为背包有后效性...我做了两次,也就是迭代了一下

剩下的30%随机化了一波。就是先把每个数的20以内的倍数暴力的算出来对k取模然后丢到一个大小为k的桶里面去。因为题目就是让你给每个数一个系数,于是我就每次随机两个位置相加判断在模k的意义下是否出现过,如果没有出现过就加入答案中,咳咳重复1e7次即可A掉本题

下面说说题解做法:

$ax+by=z$存在整数解,当且仅当$gcd(a, b)∣z$。

那么,若z可以被凑出,即 $\sum_{i=1}^{n} x_ia_i = z$,当且仅当 $gcd(a_1, a_2,⋯, a_n)∣z$。

因此,答案只能是gcd的整数倍。

但是,这样考虑x有可能是负数,但是在mod k的条件下,我们可以把x调为非负整数。

时间复杂度$O((n + k)logv)$。

乱搞代码

#include<algorithm>
#include<cstring>
#include<cstdio>
#include<iostream>
#include<vector>
#include<time.h>
using namespace std; const int N=1e6+;
int n,k;
int a[N],f[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;
}
namespace task1
{
void main()
{
for (int i=;i<=n;i++)
for (int j=;!f[1ll*j*a[i]%k];j++) f[1ll*j*a[i]%k]=;
for (int i=;i<=n;i++)
{
for (int j=;j<=k;j++) f[j]|=f[((j-a[i])%k+k)%k];
for (int j=;j<=k;j++) f[j]|=f[((j-a[i])%k+k)%k];
}
int s=;
for (int j=;j<k;j++) s+=f[j];
printf("%d\n",s);
for (int j=;j<k;j++) if (f[j]) printf("%d ",j);
}
}
int gcd(int a,int b) {if (!b) return a;else return gcd(b,a%b);}
namespace task2
{
void main()
{
printf("%d\n",k);
for (int i=;i<k;i++) printf("%d ",i);
}
}
int main()
{
n=read();k=read();
for (int i=;i<=n;i++) a[i]=read()%k;
for (int i=;i<=n;i++) if (a[i]==||gcd(a[i],k)==) {task2::main();return ;}
if (n<=) {task1::main();return ;}
srand(time());
vector <int> p;
for (int i=;i<=n;i++)
{
for (int j=;j<=;j++)
{
int q=1ll*j*a[i]%k;
if (!f[q])
{
p.push_back(q);
f[q]=;
}
}
}
for (int i=;i<=1e7;i++)
{
int si=p.size();
int l=rand()%si,r=rand()%si;
if (!f[(p[l]+p[r])%k])
{
p.push_back((p[l]+p[r])%k);
f[(p[l]+p[r])%k]=;
}
}
int s=;
for (int j=;j<k;j++) s+=f[j];
printf("%d\n",s);
for (int j=;j<k;j++) if (f[j]) printf("%d ",j);
return ;
}

T2:biology

题目链接:

http://zhengruioi.com/contest/156/problem/472

题解:

我们显然可以把元素按$a$排序,然后宽搜转移。

$f_{x,y}$ 表示当前路径的结尾在$(x,y)$位置的最大吸引度之和。 $f_{x,y} = b_{x,y}+max (f_{z,k} + ∣x −z ∣ + ∣y −k∣, a_{x,y} >a_{z,k} )$

暴力转移时间复杂度最差为$O(n ^2 m^2)$,考虑优化。

坐标转化$(x,y)->(x+y,x-y)$

这样原来的曼哈顿距离$|x-z|+|y-k|$就变成了切比雪夫距离$max(|x-z|,|y-k|)=max(x-z,z-x,y-k,k-y)$

这样坐标是最大值,转移也是最大值,因此可以用4个变量分别记录最大的$f_{z,k}-z,f_{z,k}+z,f_{z,k}-k,f_{z,k}+k$

时间复杂度为排序复杂度$O(nmlognm)$

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<queue>
#include<cmath>
using namespace std;
typedef long long ll; const int N=2e3+;
const ll inf=1e9;
int n,m,tot;
ll A,B,C,D;
ll b[N][N],dp[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;
}
struct node{
int x,y;
int d;
}s[N*N];
bool operator < (node x,node y) {return x.d<y.d;}
void chkmx(ll &a,ll b) {if (b>a) a=b;}
void calc()
{
int l=;ll mx=;
A=-inf;B=-inf;C=-inf;D=-inf;
while (l<=tot)
{
int r=l;
while (r<tot&&s[r].d==s[r+].d) ++r;
for (int k=l;k<=r;k++)
{
int x=s[k].x,y=s[k].y;
int i=x+y,j=x-y;
dp[x][y]=b[x][y];
chkmx(dp[x][y],b[x][y]+A+i);
chkmx(dp[x][y],b[x][y]+B-i);
chkmx(dp[x][y],b[x][y]+C+j);
chkmx(dp[x][y],b[x][y]+D-j);
chkmx(mx,dp[x][y]);
}
for (int k=l;k<=r;k++)
{
int x=s[k].x,y=s[k].y;
int i=x+y,j=x-y;
chkmx(A,dp[x][y]-i);
chkmx(B,dp[x][y]+i);
chkmx(C,dp[x][y]-j);
chkmx(D,dp[x][y]+j);
}
l=r+;
}
printf("%lld\n",mx);
}
int main()
{
n=read();m=read();
for (int i=;i<=n;i++)
for (int j=,x;j<=m;j++)
{
x=read();
if (x) s[++tot]=(node){i,j,x};
}
for (int i=;i<=n;i++)
for (int j=;j<=m;j++) b[i][j]=read();
sort(s+,s++tot);
calc();
return ;
}

正睿NOIP赠送附加赛1的更多相关文章

  1. 10.25 正睿停课训练 Day9

    目录 2018.10.25 正睿停课训练 Day9 A 数独(思路 DP) B 红绿灯(最短路Dijkstra) C 轰炸(计算几何 圆并) 考试代码 B C 2018.10.25 正睿停课训练 Da ...

  2. 11.6 正睿停课训练 Day17

    目录 2018.11.6 正睿停课训练 Day17 A chinese(思路 计数) B physics(单调队列/剪枝 DP) C chemistry(期望 DP) 考试代码 A B C 2018. ...

  3. 8.10 正睿暑期集训营 Day7

    目录 2018.8.10 正睿暑期集训营 Day7 总结 A 花园(思路) B 归来(Tarjan 拓扑) C 机场(凸函数 点分治) 考试代码 A B C 2018.8.10 正睿暑期集训营 Day ...

  4. 10.31 正睿停课训练 Day13

    目录 2018.10.31 正睿停课训练 Day13 A Poker(期望) B Label(高斯消元) C Coin(二分图染色 博弈) 考试代码 A(打表) B 2018.10.31 正睿停课训练 ...

  5. 11.5 正睿停课训练 Day16

    目录 2018.11.5 正睿停课训练 Day16 A 道路规划(思路) B 逻辑判断(枚举 位运算/DP 高维前缀和) C 区间(贪心/树状数组) 考试代码 A B C 2018.11.5 正睿停课 ...

  6. 11.2 正睿停课训练 Day15

    目录 2018.11.2 正睿停课训练 Day15 A 郁闷的小G(二分) B 小G的树(树形DP) C 数的距离(思路) 考试代码 B C 2018.11.2 正睿停课训练 Day15 时间:3.5 ...

  7. 11.1 正睿停课训练 Day14

    目录 2018.11.1 正睿停课训练 Day14 A 字符串 B 取数游戏(贪心) C 魔方(模拟) 考试代码 B C 2018.11.1 正睿停课训练 Day14 时间:3.5h 期望得分:100 ...

  8. 10.29 正睿停课训练 Day11

    目录 2018.10.29 正睿停课训练 Day11 A 线段树什么的最讨厌了(思路 DFS) B 已经没有什么好害怕的了(差分 前缀和) C 我才不是萝莉控呢(DP 贪心 哈夫曼树) 考试代码 A ...

  9. 10.30 正睿停课训练 Day12

    目录 2018.10.30 正睿停课训练 Day12 A 强军战歌(DP 树状数组 容斥) B 当那一天来临(思路) C 假如战争今天爆发(贪心) 考试代码 B C 2018.10.30 正睿停课训练 ...

随机推荐

  1. Bootstrap表格内容居中

    1.<th style='text-align: center;'>host</th> 水平居中 2.<td rowspan=$row_host1 style='vert ...

  2. Java Servlet 配置

    图片太大,可以右键另存再查看,也可以鼠标点击拖置一下,用浏览器单独承载放大查看.

  3. LeetCode(11)Container With Most Water

    题目如下: 题目的意思是求容器能装的最大的水量,当时我按梯形的面积来算,一直不对,后来才发现要按矩形的面积来算 Python代码如下: def maxArea(self, height): " ...

  4. C# 正则表达式大全(转载)

    文章导读 正则表达式的本质是使用一系列特殊字符模式,来表示某一类字符串.正则表达式无疑是处理文本最有力的工具,而.NET提供的Regex类实现了验证正则表达式的方法.Regex 类表示不可变(只读)的 ...

  5. BZOJ 3282 Link Cut Tree (LCT)

    题目大意:维护一个森林,支持边的断,连,修改某个点的权值,求树链所有点点权的异或和 洛谷P3690传送门 搞了一个下午终于明白了LCT的原理 #include <cstdio> #incl ...

  6. Vue的数据依赖实现原理简析

    首先让我们从最简单的一个实例Vue入手: const app = new Vue({ // options 传入一个选项obj.这个obj即对于这个vue实例的初始化 }) 通过查阅文档,我们可以知道 ...

  7. VUE:渐进式JavaScript框架(小白自学)

    VUE:渐进式JavaScript框架 一.官网 英文 https://vuejs.org/ 中文 https://cn.vuejs.org/ 二:渐进式 即有一个核心库,在需要的时候再逐渐添加插件的 ...

  8. IDEA使用快捷键

     sout+TAB键---->System.out.println();你可以按ctrl+j里面各种快捷键模板都可以看到. Intellij Idea get/set方法快捷键:Alt+Inse ...

  9. UVALIVE 4256 Salesmen

    Salesmen Time Limit: 3000ms Memory Limit: 131072KB This problem will be judged on UVALive. Original ...

  10. 楼宇自控-BA系统流程总图

    总结一下过程中的节点和技能,希望能对其他人有所帮助