【奢侈的旅行】
老博客关了之后就再也没有写过科学的dijstra+heap,稍微mark一下
#include<cstdio>
#include<cmath>
#include<cstring>
#include<queue>
#include<algorithm>
#include<iostream>
#define ll long long
#define N 404040
using namespace std;
const double eps=1e-13;
int head[N],vet[N],jie[N],edgenum,inq[N],n,cas,m;
ll pri[N],dis[N],cons[N];
struct Node{
int id;
ll dis;
bool operator <(const Node &a)const{
return dis>a.dis;
}
}tt;
priority_queue<Node>q;
void add(int u,int v,int a,int b){
edgenum++;vet[edgenum]=v;jie[edgenum]=head[u];
head[u]=edgenum;pri[edgenum]=a;cons[edgenum]=b;
}
int ok(int u,int v,ll al,ll b,ll a){
if(b==0)return 1;
b=1ll<<b;b--;
if(b<=a/al)return 1;else return 0;
}
void spfa(){
for(int i=0;i<=n;i++)dis[i]=1ll<<61;
while(!q.empty())q.pop();
dis[1]=1;
tt.id=1;
tt.dis=1;
q.push(tt);
while(!q.empty()){
Node node=q.top();
q.pop();
int u=node.id,e=head[u];
if(inq[u]) continue;
inq[u]=1;
while(e>0){
int v=vet[e];
if(dis[u]+pri[e]<dis[v] &&
ok(u,v,dis[u],cons[e],pri[e])&& inq[v]==0){
dis[v]=dis[u]+pri[e];
tt.id=v;
tt.dis=dis[v];
q.push(tt);
}
e=jie[e];
}
}
}
int main(){
scanf("%d",&cas);
while(cas--){
scanf("%d%d",&n,&m);
edgenum=0;
for(int i=1;i<=n;i++)head[i]=0,inq[i]=0;
for(int i=1;i<=m;i++){
int u,v,a,b;
scanf("%d%d%d%d",&u,&v,&a,&b);
add(u,v,a,b);
}
spfa();
double tt=log((double)dis[n])/log(2.0);
if(dis[n]==dis[0])printf("-1\n");else{
printf("%lld\n",(ll)(tt+eps));
}
}
}
【寻宝游戏】
这个dp[i][j][x][y],表示走到(i,j)位置时考虑了 前面(i-1)*m+j个 点,这是被我们没想到的最主要因素。没有其它的状态了,如果想要正确地处理一些交换操作的话)
类似精妙的dp想到了寒假只做出了miner一个人的那道状压题,都是别人用力写写很套路自己却找不到套路怎么走的题。。
至于敲代码。。冷静的力量是无限的。。~
【quality算法】
一道被我们队从开场开到结尾的题。题解长得特别虚张声势,还特地跑去查了“拟阵”是什么意思。。
不过这个题,不用看出来那些东西也是可以强做的,但我猜出题人的过程可能是从最小基环生成树出发然后加了一点包装上来造出了这个题)
考虑对所有的点权按位分治,必定高位相同的权值需要尽可能的组内配对。然后发现跨组配对时,点数大于等于三的两个组必然GG。。然后就没有了。。
【代码party】
当时通读完题目以后我是这么说的,“回文数和对称数应该不太好开”,“异或点对和矩阵题应该好开一点”,无奈赛场上没有什么好好开题的条件。
对矩阵做二维前缀和,再用一些精妙的容斥来计算方案,最后每个三元组都只会在重复区域块的左上角被记录一次><
【对称数】
呈上一些链接
https://files-cdn.cnblogs.com/files/clrs97/2018CCPC%E5%A5%B3%E7%94%9F%E8%B5%9Bsol.pdf
还有DZY loves Chinese2