CQOI2016滚粗记
day0
上去试机发现鼠标是坏的，旁边是巴蜀的大牛和教师机…
随手敲了一个树剖和Pollard_rho，顺路用批处理对拍了一下，然后想去交一发意外的发现竟然木有网….
晚上复习了一下板子…然后写得有点兴奋了就睡不着辣…
day1
上歌乐山的时候还在下雨…
然后我7点就悠悠的到了门口…然后就在外面晃荡了40多分钟…
8点过就进了考场，发现解压密码一点都不好玩…
看到T1的时候心里一阵窃喜，嘿嘿我做过bzoj2229的，不就是一道分治最小割的水题吗…30分钟敲完就写了一个对拍，拍了一个小时木有错就默默的水过辣…
#include <iostream>
#include <cstdio>
#include <algorithm>
#define MAXN 855
#define MAXM 20005
#define inf 0x3f3f3f3f
#define min(a,b) ((a)<(b)?(a):(b))
using namespace std;
struct node
{
int v, w, next;
}edge[MAXM];
int adj[MAXN], pos, n, m, a[MAXN], tmp[MAXN], T, S;
int d[MAXN], vd[MAXN], road[MAXN*MAXN], cnt;
inline void add(int u,int v,int w)
{
edge[pos].v=v, edge[pos].w=w, edge[pos].next=adj[u];
adj[u]=pos;
++pos;
edge[pos].v=u, edge[pos].w=w, edge[pos].next=adj[v];
adj[v]=pos;
++pos;
}
int aug(int u,int augco)
{
if(u==T)return augco;
int augc=augco, mind=n, delta, v;
for(int p=adj[u];~p;p=edge[p].next)
if(edge[p].w>0)
{
v=edge[p].v;
if(d[v]+1==d[u])
{
delta=aug(v,min(augc,edge[p].w));
edge[p].w-=delta, augc-=delta, edge[p^1].w+=delta;
if(!augc||d[S]>=n)return augco-augc;
}
mind=min(mind,d[v]);
}
if(augco==augc)
{
if(!--vd[d[u]])d[S]=n;
else ++vd[d[u]=mind+1];
}
return augco-augc;
}
int q[MAXM], head, tail;
bool vis[MAXN], in[MAXN];
void bfs()
{
head=tail=0;
q[++tail]=T, d[T]=0;
int u, v;
while(head<tail)
{
u=q[++head], in[u]=0;
for(int p=adj[u];~p;p=edge[p].next)
if(edge[p^1].w>0&&d[v=edge[p].v]>d[u]+1)
{
d[v]=d[u]+1;
if(!in[v])q[++tail]=v, in[v]=1;
}
}
}
int isap()
{
for(int p=0;p<pos;p+=2)
edge[p].w=edge[p^1].w=(edge[p].w+edge[p^1].w)/2;
for(int i=1;i<=n;++i)d[i]=inf, vd[i]=0, vis[i]=0;
bfs();
for(int i=1;i<=n;++i)
if(d[i]!=inf)++vd[d[i]];
int ans=0;
while(d[S]<n)
ans+=aug(S,inf);
return ans;
}
void dfs(int u)
{
vis[u]=1;
for(int p=adj[u], v;~p;p=edge[p].next)
if(edge[p].w>0&&!vis[v=edge[p].v])
dfs(v);
}
void solve(int l,int r)
{
if(l==r)return;
S=a[l], T=a[r];
int ans=isap();
road[++cnt]=ans;
dfs(S);
int p1=l, p2=r;
for(int i=l;i<=r;++i)
if(vis[a[i]])tmp[p1++]=a[i];
else tmp[p2--]=a[i];
for(int i=l;i<=r;++i)
a[i]=tmp[i];
solve(l,p1-1), solve(p2+1,r);
}
int main()
{
scanf("%d%d",&n,&m);
for(int i=1;i<=n;++i)adj[i]=-1, a[i]=i;
int u, v, w;
while(m--)
{
scanf("%d%d%d",&u,&v,&w);
add(u,v,w);
}
solve(1,n);
sort(road+1,road+cnt+1);
int ans=1;
for(int i=2;i<=cnt;++i)
if(road[i]!=road[i-1])++ans;
printf("%d\n",ans);
return 0;
}T2 不会的说…就默默滴写了个大暴力…外加一个分治乱写一共骗了40分…
附上之后写的题解… 题解链接
T3 一眼看穿是一个数位DP，然后…我TM搞了两个多小时都木有搞出来，真是太垃圾辣…
day1就这样160滚粗辣…
晚上就默默滴敲了几个计算几何和数学的板子…
day2
早上继续雨蒙蒙…
怀着对滚粗的美好向往，我默默滴进入机房…
然后…
T1 TM不就是想考Pollard_rho吗…不虚…
30分钟敲完，又写了一个对拍，没问题…
恩，100分稳稳滴…
#include <iostream>
#include <cstdio>
#include <algorithm>
#define abs(a) ((a)>0?(a):-(a))
#define LL long long int
using namespace std;
LL mul(LL a,LL pos,LL mod)
{
LL ans=0;
for(;pos;pos>>=1)
{
if(pos&1){ans+=a;if(ans>=mod)ans-=mod;}
a+=a;if(a>=mod)a-=mod;
}
return ans;
}
LL power(LL a,LL pos,LL mod)
{
LL ans=1;
for(;pos;pos>>=1, a=mul(a,a,mod))
if(pos&1)ans=mul(ans,a,mod);
return ans;
}
bool Miller_Rabin(LL n)
{
if(n==2||n==3)return 1;
if(!(n&1))return 0;
LL m=n-1, a, tmp, ans;
int cnt=0;
while(!(m&1))m>>=1, ++cnt;
for(int i=0;i<20;++i)
{
a=rand()%(n-1)+1;
tmp=power(a,m,n);
for(int j=0;j<cnt;++j)
{
ans=mul(tmp,tmp,n);
if(ans==1)
{
if(tmp!=1&&tmp!=n-1)return 0;
break;
}
tmp=ans;
}
if(ans!=1)return 0;
}
return 1;
}
LL gcd(LL a,LL b)
{
if(!b)return a;
return gcd(b,a%b);
}
LL work(LL a,LL b)
{
LL cnt=0, k=2, x0=rand()%a, y=x0, d;
while(1)
{
++cnt;
x0=(mul(x0,x0,a)+b)%a;
d=gcd(abs(y-x0),a);
if(d!=1&&d!=a)return d;
if(y==x0)return a;
if(cnt==k)k+=k, y=x0;
}
}
LL fac[105];
int tot;
void dfs(LL n)
{
if(Miller_Rabin(n)){fac[++tot]=n;return;}
LL p=n;
while(p>=n)p=work(p,rand()%(n-1)+1);
dfs(p), dfs(n/p);
}
void exgcd(LL a,LL b,LL &d,LL &x,LL &y)
{
if(!b)d=a, x=1, y=0;
else
{
exgcd(b,a%b,d,y,x);
y-=a/b*x;
}
}
LL ask_ni(LL a,LL b)
{
LL d, x, y;
exgcd(a,b,d,x,y);
LL mod=b/d;
return (x%mod+mod)%mod;
}
LL e, N, c;
int main()
{
cin>>e>>N>>c;
dfs(N);
LL r=(fac[1]-1)*(fac[2]-1);
LL d=ask_ni(e,r);
cout<<d<<' '<<power(c,d,N);
return 0;
}T2 一开始木有看懂题啊…
然后磨了半个小时，然后加了一点补充说明就看懂了。
不就是一个水水的01Trie吗，不虚200稳稳滴…
#include <iostream>
#include <cstdio>
#include <algorithm>
#define MAXN 1000005
#define MAXM 15050000
using namespace std;
char w;
inline void GET(int &t)
{
t=0;
do w=getchar();while(w<'0'||w>'9');
do{t=t*10+w-'0';w=getchar();}while(w>='0'&&w<='9');
}
int tmp[MAXN], temp, len[MAXN];
struct Trie
{
int ch[MAXM][2], pos, id[MAXM];
inline void insert(char w[],int len,int d)
{
int p=0;
for(int i=1;i<=len;++i)
if(!ch[p][w[i]])p=ch[p][w[i]]=++pos;
else p=ch[p][w[i]];
if(!id[p])id[p]=d;
}
inline void query(char w[])
{
int p=0;
for(int i=1;i<33;++i)
{
if(ch[p][w[i]])p=ch[p][w[i]];
else return;
if(id[p])tmp[++temp]=id[p];
}
}
}t;
char s[36];
int main()
{
int n, u, a, b, cnt=0, ans, mv, tot;
char ask[10];
GET(n);
while(n--)
{
scanf("%s",ask);
if(ask[0]=='A')
{
++cnt, tot=0;
for(int i=1;i<5;++i)
{
GET(u);
for(int j=7;~j;--j)
if(u&(1<<j))s[++tot]=1;
else s[++tot]=0;
}
GET(len[cnt]);
t.insert(s,len[cnt],cnt);
}
else
{
tot=0;
for(int i=1;i<5;++i)
{
GET(u);
for(int j=7;~j;--j)
if(u&(1<<j))s[++tot]=1;
else s[++tot]=0;
}
temp=0, t.query(s);
sort(tmp+1,tmp+temp+1);
GET(a), GET(b), ans=mv=0;
for(int i=1;i<=temp;++i)
if(mv<len[tmp[i]])
{
mv=len[tmp[i]];
if(tmp[i]>b)break;
if(a<=tmp[i])++ans;
}
printf("%d\n",ans);
}
}
return 0;
}T3 一看，(⊙v⊙)嗯好像是一道赤裸裸数学题，然后就默默滴推了2个多小时，然后木有推出来，然后就写了一个大暴力，弃疗…
正解：大力出奇迹…我真是…
把128内的质数打出来，随便搞搞这个搞出来这个质数最多能被选多少次，然后枚举一个要选择的最大的质数，这样就可以枚举这个质数被选次数内的比这个质数小的所有质数了。就是个集合的无序拆分。
剪枝：10^14的第800000个伪光滑数一定比10^15的第800000个伪光滑数小，这样假设我们知道10^14的第800000个伪光滑数，这个数可以拿来剪枝，然后10^15的第800000个伪光滑数算起来就很快了，可以拿去剪10^16的第800000个伪光滑数。
by quack…
我就%%%%%%
day2 230彻底滚粗…
加上noip成绩默默滴排到了第8左右…
幸好noip考得好…
然后女生保护政策就默默滴拿了一个A类…
感觉不是很爽，毕竟D1T2出题人面向巴蜀选手出数据，居然让旋转卡壳和三分就这样水过辣，表示很不爽…
SB出题人！！！SB出题人！！！SB出题人！！！
重要的事情说三遍…