noip2013一开始打了260分。。。感觉最后一题实在太麻烦考试的时候也写不出来,,,打个60分也不错了。
1.一开始被数据吓了一跳。。。然后想到肯定是快速幂。但我有点脑残的先用gcd求了一个循环节。其实只要mod m就可以了脑残脑残。还要提醒一下能不用mod尽量不用。mod的常数真的很大。
代码:#include<iostream>
#include<cmath>
#include<cstdio>
using namespace std;
long long n,m,k,x;
long long p;
long long gcd(long long a,long long b)
{ if(b==0) return a;
else return gcd(b,a%b);
}
long long mi(long long w)
{ if(w==1) return 10;
if(w%2==0)
{long long r=mi(w/2)%p;
return (r*r)%p;
}
else
{ long long r=mi(w/2)%p;
return (r*r*10)%p;
}
}
int main()
{
freopen("circle.in","r",stdin);
freopen("circle.out","w",stdout);
cin>>n>>m>>k>>x;
p=gcd(n,m);
p=n*m/p;
long long ans=mi(k);
int ou=x;
for(int i=1;i<=ans;i++)
{ ou=ou+m;
if(ou>=n) ou=ou-n; //本来这里写的是取模。。结果就炸了两个点。。。
}
cout<<ou;
return 0;
}
2.直接告诉我们肯定是两边都从小到大匹配。。。不想证明了、、然后我们假设原来的a数组相对不懂。。根据排序后得出的编号匹配得出b数组的目标状态。然后逆序对数就是最少交换次数。。还是不想证明。。。还有一点有点意思。match也有匹配的意思。所以就要往匹配上想。。。
代码:#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;
int d[1000001];
long long ans=0;
void stablesort(int *a,int l,int mid,int r)
{
int i=l,j=mid+1,k=l;
while (i<=mid && j<=r)
if (a[i]<=a[j]) d[k++]=a[i++];
else { d[k++]=a[j++]; ans+=mid+1-i;ans=ans%99999997;}
while (i<=mid)
d[k++]=a[i++];
while (j<=r)
d[k++]=a[j++];
for (int q=l;q<=r;q++) a[q]=d[q];
}
void merge(int *a,int l,int r)
{
if (l<r)
{
int mid=(l+r)/2;
merge(a,l,mid);
merge(a,mid+1,r);
stablesort(a,l,mid,r);
}
}
int n,c[1000001];
struct data{
int xx,yy;}a[1000001],b[1000001];
bool cmp(data a,data b)
{
return a.xx<b.xx;
}
int main()
{
freopen("match.in","r",stdin);
freopen("match.out","w",stdout);
scanf("%d",&n);
for (int q=0;q<n;q++) {scanf("%d",&a[q].xx);a[q].yy=q;}
for (int q=0;q<n;q++) {scanf("%d",&b[q].xx);b[q].yy=q;}
sort(a,a+n,cmp);
sort(b,b+n,cmp);
for (int q=0;q<n;q++)
c[b[q].yy]=a[q].yy;
merge(c,0,n-1);
printf("%d",ans);
return 0;
}
第三题。。最大生成树肯定是的、、问题就是在怎么做。。。我想到了求最近公共祖先。。求出每个点的深度,x,y,两个谁深度较大就把谁向上找爸爸。。。这样拿了60分。。。真心不错了!!!#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<vector>
#include<set>
#include<iterator>
#include<algorithm>
#include<climits>
using namespace std;
int n,m,q;
int f[10001];
int fa[10001];
int lenfa[10001];
int deep[1001];
bool flag[1001];
struct edge
{ int from,to,len;
};
vector<edge> e;
vector<edge> g[10001];
bool comp(edge x,edge y)
{ return x.len>y.len;
}
int get(int x)
{if(f[x]==x) return x;
f[x]=get(f[x]);
return f[x];
}
void dfs(int x,int tot)
{
for(int i=0;i<g[x].size();i++)
{ int u=g[x][i].to;
if(flag[u]==0)
{flag[u]=1;
fa[u]=x;
lenfa[u]=g[x][i].len;
deep[u]=tot+1;
dfs(u,tot+1);
}
}
}
int main()
{
freopen("truck.in","r",stdin);
freopen("truck.out","w",stdout);
cin>>n>>m;
for(int i=1;i<=n;i++) {f[i]=i;fa[i]=i;}
for(int i=1;i<=m;i++)
{ int x,y,z; scanf("%d",&x); scanf("%d",&y);scanf("%d",&z);
e.push_back((edge){x,y,z});
}
sort(e.begin(),e.end(),comp);
int sum=0;int ans=0;
for(int i=0;i<=m-1;i++)
{ int a=get(e[i].to);int b=get(e[i].from);
if(a!=b)
{sum++;
f[a]=b;
g[e[i].to].push_back((edge){e[i].to,e[i].from,e[i].len});
g[e[i].from].push_back((edge){e[i].from,e[i].to,e[i].len});
}
if(sum==n-1) break;
}
for(int i=1;i<=n;i++)
{ if(f[i]==i)
{memset(flag,0,sizeof(flag));
flag[i]=1;deep[i]=0;
dfs(i,0);
}
}
cin>>q;
for(int i=1;i<=q;i++)
{int x,y;scanf("%d %d",&x,&y);
int cmin=INT_MAX/10;
while(x!=y&&(deep[x]!=0||deep[y]!=0))
{if(deep[x]>deep[y])
{cmin=min(cmin,lenfa[x]);
x=fa[x];
}
else
{cmin=min(cmin,lenfa[y]);
y=fa[y];
}
}
if(cmin==INT_MAX/10)
printf("-1\n");
else
printf("%d\n",cmin);
}
return 0;
}
好了给个满分程序。。。反正我是写不出来的树上倍增
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#define inf 0x7fffffff
using namespace std;
int n,m,q,cnt,tot,deep[10001],head[10001],f[10001],fa[10001][17],d[10001][17];
bool vis[10001];
struct edge{int x,y,v;}a[50001];
struct e{int next,to,v;}e[20001];
void ins(int u,int v,int w)
{e[++cnt].to=v;e[cnt].next=head[u];head[u]=cnt;e[cnt].v=w;}
void insert(int u,int v,int w)
{ins(u,v,w);ins(v,u,w);}
int find(int x)
{return f[x]==x?x:f[x]=find(f[x]);}
bool cmp(edge a,edge b)
{return a.v>b.v;}
void dfs(int x)
{
vis[x]=1;
for(int i=1;i<=16;i++)
{
if(deep[x]<(1<<i))break;
fa[x][i]=fa[fa[x][i-1]][i-1];
d[x][i]=min(d[x][i-1],d[fa[x][i-1]][i-1]);
}
for(int i=head[x];i;i=e[i].next)
{
if(vis[e[i].to])continue;
fa[e[i].to][0]=x;
d[e[i].to][0]=e[i].v;
deep[e[i].to]=deep[x]+1;
dfs(e[i].to);
}
}
int lca(int x,int y)
{
if(deep[x]<deep[y])swap(x,y);
int t=deep[x]-deep[y];
for(int i=0;i<=16;i++)
if((1<<i)&t)x=fa[x][i];
for(int i=16;i>=0;i--)
{
if(fa[x][i]!=fa[y][i])
{x=fa[x][i];y=fa[y][i];}
}
if(x==y)return x;
return fa[x][0];
}
int ask(int x,int f)
{
int mn=inf;
int t=deep[x]-deep[f];
for(int i=0;i<=16;i++)
{
if(t&(1<<i))
{
mn=min(mn,d[x][i]);
x=fa[x][i];
}
}
return mn;
}
int main()
{
memset(d,127/3,sizeof(d));
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)f[i]=i;
for(int i=1;i<=m;i++)
scanf("%d%d%d",&a[i].x,&a[i].y,&a[i].v);
sort(a+1,a+m+1,cmp);
for(int i=1;i<=m;i++)
{
int x=a[i].x,y=a[i].y;
int p=find(a[i].x),q=find(a[i].y);
if(p!=q)
{
f[p]=q;
insert(x,y,a[i].v);
tot++;if(tot==n-1)break;
}
}
for(int i=1;i<=n;i++)if(!vis[i])dfs(i);
scanf("%d",&q);
for(int i=1;i<=q;i++)
{
int x,y;scanf("%d%d",&x,&y);
if(find(x)!=find(y)){printf("-1\n");continue;}
else
{
int t=lca(x,y);
printf("%d\n",min(ask(x,t),ask(y,t)));
}
}
return 0;
}