10.25 正睿停课训练 Day9
2018.10.25 正睿停课训练 Day9
期望得分:100+60+20
实际得分:100+0+0
Dijkstra模板题爆零了。我还有救吗
A 数独(思路 DP)
先是想到,限制的是1,但其实在数独里1和2,3,...,9别的数没啥特殊的地方啊,可以忽略其它数的种类?
好,上面这句话没什么用。
考虑到每种合法的填1方案(不考虑其它数,先只填1),对数独的解数的贡献是一样的。
所以我们求 \(当前填1合法方案数/无限制填1合法方案数\) 这个比例,再乘以无限制的数独总解数就行了。
当前填1合法方案数可以DP,状压一下,三行三行转移,非常轻松。直接DFS也随便过。
数独总解数大概可以各种剪枝+优化(对称性)在本机跑出来。反正我选择百度。
搜索求解的数量(orz wph):
或者面向样例得到解的数量(orz wzh)
//49ms 504kb
#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
//#define gc() getchar()
#define MAXIN 100000
#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
#define mod 998244353
#define lb(x) (x&-(x))
#define Mod(x) x>=mod&&(x-=mod)
typedef long long LL;
const int N=12,M=(1<<9)+2;
int bit[M];
bool ban[N][N];
char IN[MAXIN],*SS=IN,*TT=IN;
inline int read()
{
int now=0;register char c=gc();
for(;!isdigit(c);c=gc());
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now;
}
inline LL Read()
{
LL now=0;register char c=gc();
for(;!isdigit(c);c=gc());
for(;isdigit(c);now=now*10+c-'0',now%=mod,c=gc());
return now;
}
inline int FP(int x,int k)
{
int t=1;
for(; k; k>>=1,x=1ll*x*x%mod)
if(k&1) t=1ll*t*x%mod;
return t;
}
inline bool Check(int a,int b)
{
if(a<3) return b>=3;
if(a>5) return b<=5;
return b<=2||b>=6;
}
int Calc()
{
static int f[M];
memset(f,0,sizeof f);
f[0]=1;
int all=(1<<9)-1;
for(int s=0; s<all; ++s)
{
int cnt=bit[s];
if(cnt%3) continue;
for(int a=0; a<9; ++a)
if(!(s>>a&1)&&!ban[cnt+1][a+1])
for(int b=0; b<9; ++b)
if(!(s>>b&1)&&!ban[cnt+2][b+1]&&Check(a,b))
for(int c=0; c<9; ++c)
{
if(s>>c&1||ban[cnt+3][c+1]||!Check(a,c)||!Check(b,c)) continue;
int ss=s|(1<<a)|(1<<b)|(1<<c);
f[ss]+=f[s], Mod(f[ss]);
}
}
return f[all];
}
int main()
{
const int Sum=719935075;//Read(); //6670903752021072936960
const int tot=46656,inv=549081465;// int tot=Calc(),inv=FP(tot,mod-2);
for(int i=1; i<M; ++i) bit[i]=bit[i^lb(i)]+1;
for(int T=read(); T--; )
{
for(int i=1; i<=9; ++i)
{
register char c=gc(); for(;c!='0'&&c!='1';c=gc());
ban[i][1]=c=='1';
for(int j=2; j<=9; ++j) ban[i][j]=gc()=='1';
}
int now=Calc();
printf("%d\n",(int)(1ll*Sum*now%mod*inv%mod));
}
return 0;
}
B 红绿灯(最短路Dijkstra)
不管一条边怎么样,我们到达一个点时肯定是越早越好啊(不会更差)。
所以直接Dijkstra就行了。边权就模拟一下。
//2762ms 54532kb
#include <queue>
#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
#define mp std::make_pair
#define pr std::pair<LL,int>
//#define gc() getchar()
#define MAXIN 300000
#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
typedef int LL;//longlong?
const int N=5e5+6,M=2e6+7;
int Enum,H[N],nxt[M],to[M],len[M],K[M],L[M],R[M];
char IN[MAXIN],*SS=IN,*TT=IN;
inline int read()
{
int now=0;register char c=gc();
for(;!isdigit(c);c=gc());
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now;
}
inline void AE(int r,int l,int k,int d,int v,int u)
{
to[++Enum]=v, nxt[Enum]=H[u], H[u]=Enum, len[Enum]=d, K[Enum]=k, L[Enum]=l, R[Enum]=r;
to[++Enum]=u, nxt[Enum]=H[v], H[v]=Enum, len[Enum]=d, K[Enum]=k, L[Enum]=l, R[Enum]=r;
}
inline int Calc(int now,int d,int k,int l,int r)
{
int tm=0;
if(now>r) tm+=k-now+l;
else if(now<l) tm+=l-now;
else if(r-now+1>=d) return d;
else d-=r-now+1, tm+=k-now+l;
now=(d-1)/(r-l+1), tm+=now*k, d-=now*(r-l+1);//(d-1)/len not d/len!
return tm+d;
}
void Dijkstra(int n)
{
static LL dis[N];
static bool vis[N];
static std::priority_queue<pr> q;
memset(dis,0x3f,sizeof dis);
dis[1]=0, q.push(mp(0,1));
while(!q.empty())
{
int x=q.top().second; q.pop();
if(vis[x]) continue;
vis[x]=1; LL ds=dis[x];
for(int i=H[x],v,tm; i; i=nxt[i])
if(tm=Calc(ds%K[i],len[i],K[i],L[i],R[i]),dis[v=to[i]]>ds+tm)
q.push(mp(-(dis[v]=ds+tm),v));
}
for(int i=1; i<=n; ++i) printf("%d\n",dis[i]);
}
int main()
{
Enum=1; int n=read(),m=read();
for(int i=1; i<=m; ++i) AE(read(),read(),read(),read(),read(),read());
Dijkstra(n);
return 0;
}
C 轰炸(计算几何 圆并)
求圆并,然后找覆盖\(k\)次的部分?
好不NOIP啊,咕了。
咕
考试代码
B
迷之全RE懒得调。
#include <queue>
#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
//#define gc() getchar()
#define MAXIN 300000
#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
#define Add(x,y) (x+y>=20?x+y-20:x+y)
typedef long long LL;
const int N=5e5+6,M=2e6+7;
const LL INF=0x3f3f3f3f3f3f3f3f;
int Enum,H[N],nxt[M],to[M],len[M],K[M],L[M],R[M],f[M>>1][20];
char IN[MAXIN],*SS=IN,*TT=IN;
struct Node
{
LL ds; int x,tm;
bool operator <(const Node &a)const
{
return ds>a.ds;
}
};
inline int read()
{
int now=0;register char c=gc();
for(;!isdigit(c);c=gc());
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now;
}
inline void AE(int r,int l,int k,int d,int v,int u)
{
to[++Enum]=v, nxt[Enum]=H[u], H[u]=Enum, len[Enum]=d, K[Enum]=k, L[Enum]=l, R[Enum]=r;
to[++Enum]=u, nxt[Enum]=H[v], H[v]=Enum, len[Enum]=d, K[Enum]=k, L[Enum]=l, R[Enum]=r;
for(int i=0; i<k; ++i)
{
int tmp=i,tm=0,dd=d;
if(tmp>r) tm+=l+k-tmp;
else if(tmp<l) tm+=l-tmp;
else if(tmp!=l)
{
if(r-tmp+1>=dd) {f[Enum>>1][i]=dd; continue;}
else dd-=r-tmp+1, tm+=k-tmp+l;
}
tmp=dd/(r-l+1), tm+=tmp*k, dd-=tmp*(r-l+1);
f[Enum>>1][i]=tm+dd;
}
}
void Dijkstra(int n)
{
static bool vis[N][20];
static LL dis[N][20];
static std::priority_queue<Node> q;
memset(dis,0x3f,sizeof dis);
dis[1][0]=0, q.push((Node){0,1,0});
while(!q.empty())
{
Node tmp=q.top(); q.pop();
LL ds=tmp.ds; int x=tmp.x, tm=tmp.tm;
if(vis[x][tm]) continue;
vis[x][tm]=1;
for(int i=H[x],v,use,tmp; i; i=nxt[i])
{
v=to[i], use=f[i>>1][tm%K[i]], tmp=Add(tm,use);
if(dis[v][tmp]<=ds+use) continue;
q.push((Node){dis[v][tmp]=ds+use,v,tmp});
}
}
for(int i=1; i<=n; ++i)
{
LL ans=INF;
for(int j=0; j<20; ++j) ans=std::min(ans,dis[i][j]);
printf("%lld\n",ans);
}
}
int main()
{
// freopen(".in","r",stdin);
// freopen(".out","w",stdout);
Enum=1; int n=read(),m=read();
for(int i=1; i<=m; ++i) AE(read(),read(),read(),read(),read(),read());
Dijkstra(n);
return 0;
}
C
为啥输出0.0000没分啊 明明有理有据啊
#include <cstdio>
#include <cctype>
#include <algorithm>
#define gc() getchar()
typedef double db;
const int N=105;
int n,K,Lim;
struct Point
{
int x,y,a;
}p[N];
inline int read()
{
int now=0,f=1;register char c=gc();
for(;!isdigit(c);c=='-'&&(f=-1),c=gc());
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now*f;
}
int main()
{
// freopen(".in","r",stdin);
// freopen(".out","w",stdout);
for(int T=read(); T--; )
{
n=read(),K=read(),Lim=read();
for(int i=1; i<=n; ++i) p[i]=(Point){read(),read(),read()};
if(K==1||1) printf("0.0000\n");
// else if(!Lim && K==n) Subtask3::Main();
// else if(!Lim) Subtask1::Main();
}
return 0;
}
10.25 正睿停课训练 Day9的更多相关文章
- 10.31 正睿停课训练 Day13
目录 2018.10.31 正睿停课训练 Day13 A Poker(期望) B Label(高斯消元) C Coin(二分图染色 博弈) 考试代码 A(打表) B 2018.10.31 正睿停课训练 ...
- 10.29 正睿停课训练 Day11
目录 2018.10.29 正睿停课训练 Day11 A 线段树什么的最讨厌了(思路 DFS) B 已经没有什么好害怕的了(差分 前缀和) C 我才不是萝莉控呢(DP 贪心 哈夫曼树) 考试代码 A ...
- 10.30 正睿停课训练 Day12
目录 2018.10.30 正睿停课训练 Day12 A 强军战歌(DP 树状数组 容斥) B 当那一天来临(思路) C 假如战争今天爆发(贪心) 考试代码 B C 2018.10.30 正睿停课训练 ...
- 10.24 正睿停课训练 Day8 AM
目录 2018.10.24 正睿停课训练 Day8 AM A 棒棒糖(组合) B 彩虹糖(思路 博弈) C 泡泡糖(DP) 考试代码 A B C 2018.10.24 正睿停课训练 Day8 AM 期 ...
- 10.23 正睿停课训练 Day7
目录 2018.10.23 正睿停课训练 Day7 A 矩形(组合) B 翻转(思路) C 求和(思路 三元环计数) 考试代码 B1 B2 C 2018.10.23 正睿停课训练 Day7 期望得分: ...
- 11.6 正睿停课训练 Day17
目录 2018.11.6 正睿停课训练 Day17 A chinese(思路 计数) B physics(单调队列/剪枝 DP) C chemistry(期望 DP) 考试代码 A B C 2018. ...
- 11.5 正睿停课训练 Day16
目录 2018.11.5 正睿停课训练 Day16 A 道路规划(思路) B 逻辑判断(枚举 位运算/DP 高维前缀和) C 区间(贪心/树状数组) 考试代码 A B C 2018.11.5 正睿停课 ...
- 11.2 正睿停课训练 Day15
目录 2018.11.2 正睿停课训练 Day15 A 郁闷的小G(二分) B 小G的树(树形DP) C 数的距离(思路) 考试代码 B C 2018.11.2 正睿停课训练 Day15 时间:3.5 ...
- 11.1 正睿停课训练 Day14
目录 2018.11.1 正睿停课训练 Day14 A 字符串 B 取数游戏(贪心) C 魔方(模拟) 考试代码 B C 2018.11.1 正睿停课训练 Day14 时间:3.5h 期望得分:100 ...
随机推荐
- linux系统下安装tomcat及配置
一.下载TOMCAT压缩包 apache-tomcat-6.0.44-client.zip 或 apache-tomcat-7.0.69.tar.gz 点击进入官网:http://tomcat.apa ...
- 量化投资与Python之NumPy
数组计算 NumPy是高性能科学计算和数据分析的基础包.它是pandas等其他各种工具的基础.NumPy的主要功能:ndarray,一个多维数组结构,高效且节省空间无需循环对整组数据进行快速运算的 ...
- 关于nginx报错/usr/share/nginx/html/jiankongshare" failed (2: No such file or directory)的问题解决
nginx的location虚拟目录配置: monitor.conf server { server_name monitor.chinasoft.com; server_ ...
- 域名调整 SEO优化(nginx)
=============================================== 2019/3/31_第1次修改 ccb_warlock == ...
- 10-关于DOM的事件操作
一.JavaScript的组成 JavaScript基础分为三个部分: ECMAScript:JavaScript的语法标准.包括变量.表达式.运算符.函数.if语句.for语句等. DOM:文档对象 ...
- js闭包实例汇总
本文是通过实例来帮助大家深刻理解js闭包,是篇非常不错的文章,这里推荐给大家,有需要的小伙伴可以参考下 Js闭包 闭包前要了解的知识 1. 函数作用域 (1).Js语言特殊之处在于函数内部可以直接读取 ...
- CCF2014123集合竞价(C语言版)
问题描述 某股票交易所请你编写一个程序,根据开盘前客户提交的订单来确定某特定股票的开盘价和开盘成交量. 该程序的输入由很多行构成,每一行为一条记录,记录可能有以下几种: 1. buy p s 表示一个 ...
- CCF2014032窗口(C语言)
问题描述 在某图形操作系统中,有 N 个窗口,每个窗口都是一个两边与坐标轴分别平行的矩形区域.窗口的边界上的点也属于该窗口.窗口之间有层次的区别,在多于一个窗口重叠的区域里,只会显示位于顶层的窗口里的 ...
- AlexNet
AlexNet学习笔记 目录 AlexNet整体结构 CNN 全连接 TensorFlow实现 AlexNet是2012年ImageNet竞赛冠军获得者Hinton和他的学生Alex Krizhevs ...
- hdu3415 单调队列模板题
比较裸的单调队列 先求前缀和,枚举所有结束位置1~n+k即可 #include<iostream> #include<cstdio> #include<cstring&g ...