BUPT 2014新生暑假个人排位赛11
BOJ 488 小妹妹个数
恶意题
用桶分法去做，先用第一段代码打表至某txt文件，然后再复制到第二代码的数组部分（0（1）分桶）
#include <cstdio>
const long long MAX = 1000000100;
const long long N = 1000000000;
const long long BASE = 10000;
bool used[MAX];
int main(void)
{
freopen("out.txt","w",stdout);
int cnt = 0;
for(long long i = 2; i <= N; ++i){
if(!used[i]){
cnt++;
for(long long j = i *i; j <= N; j += i)
used[j] = true;
}
if(i % BASE == 0)
printf("%d,",cnt);
}
return 0;
}
#include <algorithm>
#include <iostream>
#include <iomanip>
#include <cstring>
#include <climits>
#include <complex>
#include <fstream>
#include <cassert>
#include <cstdio>
#include <bitset>
#include <vector>
#include <deque>
#include <queue>
#include <stack>
#include <ctime>
#include <set>
#include <map>
#include <cmath>
#define eps 1e-9
#define INF 0x3f3f3f3f
using namespace std;
typedef long long ll;
typedef long double ld;
typedef pair<ll, ll> pll;
typedef complex<ld> point;
typedef pair<int, int> pii;
typedef pair<pii, int> piii;
template<class T>
inline bool read(T &n)
{
T x = 0, tmp = 1;
char c = getchar();
while((c < '0' || c > '9') && c != '-' && c != EOF) c = getchar();
if(c == EOF) return false;
if(c == '-') c = getchar(), tmp = -1;
while(c >= '0' && c <= '9') x *= 10, x += (c - '0'),c = getchar();
n = x*tmp;
return true;
}
template <class T>
inline void write(T n)
{
if(n < 0)
{
putchar('-');
n = -n;
}
int len = 0,data[20];
while(n)
{
data[len++] = n%10;
n /= 10;
}
if(!len) data[len++] = 0;
while(len--) putchar(data[len]+48);
}
//-----------------------------------
bool Prime(int x)
{
if (x == 1)
return false;
else if (x <= 0)
return false;
else if (x == 2)
return true;
else
{
for(ll i = 2; i*i <= x; i++)
if (x%i == 0)
return false;
return true;
}
}
const long long MAX = 1000000100;
const long long N = 1000000000;
const long long BASE = 10000;
int a[]= {};
int main()
{
int n;
// freopen("data.txt","r",stdin);
while(read(n))
{
int temp=n/BASE;
int ans;
if(temp==0)
ans=0;
else
ans=a[temp-1];
for(int i=BASE*temp; i<=n; i++)
if(Prime(i))
ans++;
printf("%d\n",ans);
}
return 0;
}
BOJ 480 田田背课文
最长重复子串的长度+1
可以作为模板题来解，反正方法特别多吧
#include <algorithm>
#include <iostream>
#include <iomanip>
#include <cstring>
#include <climits>
#include <complex>
#include <fstream>
#include <cassert>
#include <cstdio>
#include <bitset>
#include <vector>
#include <deque>
#include <queue>
#include <stack>
#include <ctime>
#include <set>
#include <map>
#include <cmath>
#define EPS 1e-10
#define INF 0x3f3f3f3f
#define ll long long
using namespace std;
template<class T>
inline bool read(T &n)
{
T x = 0, tmp = 1; char c = getchar();
while((c < '0' || c > '9') && c != '-' && c != EOF) c = getchar();
if(c == EOF) return false;
if(c == '-') c = getchar(), tmp = -1;
while(c >= '0' && c <= '9') x *= 10, x += (c - '0'),c = getchar();
n = x*tmp;
return true;
}
template <class T>
inline void write(T n) {
if(n < 0) {
putchar('-');
n = -n;
}
int len = 0,data[20];
while(n) {
data[len++] = n%10;
n /= 10;
}
if(!len) data[len++] = 0;
while(len--) putchar(data[len]+48);
}
//-----------------------------------------------------------------------
const int MAXN=1000010;
void radix(int *str,int *a,int *b,int n,int m)
{
static int count[MAXN];
memset(count,0,sizeof(count));
for(int i=0;i<n;++i)
++count[str[a[i]]];
for(int i=1;i<=m;++i)
count[i]+=count[i-1];
for(int i=n-1;i>=0;--i)
b[--count[str[a[i]]]]=a[i];
}
void suffix_array(int *str,int *sa,int n,int m)
{
static int rank[MAXN],a[MAXN],b[MAXN];
for(int i=0;i<n;i++)
rank[i]=i;
radix(str,rank,sa,n,m);
rank[sa[0]]=0;
for(int i=1;i<n;i++)
rank[sa[i]]=rank[sa[i-1]]+(str[sa[i]]!=str[sa[i-1]]);
for(int i=0;1<<i <n;++i)
{
for(int j=0;j<n;++j)
{
a[j]=rank[j]+1;
b[j]=j +(1<<i)>=n? 0:rank[j+(1<<i)] +1;
sa[j]=j;
}
radix(b,sa,rank,n,n);
radix(a,rank,sa,n,n);
rank[sa[0]]=0;
for(int j=1;j<n;j++)
{
rank[sa[j]]=rank[sa[j-1]]+(a[sa[j-1]]!=a[sa[j]]||b[sa[j-1]]!=b[sa[j]]);
}
}
}
string duplicate_substr(string str)
{
string rev;
static int s[MAXN],sa[MAXN],rank[MAXN],h[MAXN];
int n=str.length();
copy(str.begin(),str.end(),s);
suffix_array(s,sa,n,256);
for(int i=0;i<n;++i)
rank[sa[i]]=i;
int k=0;
int ans1=0,pos1=0;
for(int i=0;i<n;i++)
{
k= k==0? 0:k-1;
while(rank[i]>0&&s[i+k]==s[sa[rank[i]-1]+k])
++k;
h[rank[i]]=k;
if(h[rank[i]]>ans1)
{
ans1=h[rank[i]];
pos1=i;
}
}
return str.substr(pos1,ans1);
}
int main()
{
string s;
while(cin>>s)
{
s=duplicate_substr(s);
printf("%d\n",s.length()+1);
}
return 0;
}
BOJ 478 小妹妹采蘑菇
数学题利用期望的几何概型解题：第一个取到的蘑菇必定为新蘑菇，第二个取到的有 ( n - 1 ) / n 的概率取到新蘑菇，则可以认为取一次取到了 ( n - 1 ) / n 的新蘑菇，则取到一个新蘑菇的次数是 n / ( n - 1 ) 。则第 S 次取到新蘑菇的期望是 n / ( n - S )
利用组数+高数+数列知识：假设现在我们已经取到了 S 个蘑菇 ，则 下一次取到新蘑菇的概率为是 ( n - S ) / n ，取到老蘑菇的概率为 S / n ：
第一次取到新蘑菇的期望是 1 * (n - S) / n ；
第二次取到新蘑菇的期望是 2 * (n - S) / n * ( S / n )^1 ；
第三次取到新蘑菇的期望是 3 * (n - S) / n * ( S / n )^2 ；
....
趋于无穷大的时候（用数列的错位相减，以及取极限的方法可得出，第 S 次取到新蘑菇的期望是 n / ( n - S )
#include <algorithm>
#include <iostream>
#include <iomanip>
#include <cstring>
#include <climits>
#include <complex>
#include <fstream>
#include <cassert>
#include <cstdio>
#include <bitset>
#include <vector>
#include <deque>
#include <queue>
#include <stack>
#include <ctime>
#include <set>
#include <map>
#include <cmath>
#define EPS 1e-10
#define INF 0x3f3f3f3f
#define ll long long
using namespace std;
template<class T>
inline bool read(T &n)
{
T x = 0, tmp = 1; char c = getchar();
while((c < '0' || c > '9') && c != '-' && c != EOF) c = getchar();
if(c == EOF) return false;
if(c == '-') c = getchar(), tmp = -1;
while(c >= '0' && c <= '9') x *= 10, x += (c - '0'),c = getchar();
n = x*tmp;
return true;
}
template <class T>
inline void write(T n) {
if(n < 0) {
putchar('-');
n = -n;
}
int len = 0,data[20];
while(n) {
data[len++] = n%10;
n /= 10;
}
if(!len) data[len++] = 0;
while(len--) putchar(data[len]+48);
}
//-----------------------------------------------------------------------
double dp[111];
int main()
{
int n;
while(read(n))
{
memset(dp,0.0,sizeof(dp));
for(int i=1;i<=n;i++)
{
dp[i]=dp[i-1]+(double)(n)/(double)(n-i+1);
}
printf("%.6f\n",dp[n]);
}
return 0;
}
BOJ 482 焦级长的激光炮
二分图的最有匹配
将 n 个激光炮拆成 n * m 个点，代表着每个点第 i 次开炮消耗的时间（最多开m次
以第一个点拆后的点为例建图：第一次开炮与目标点的权值为——T1（化成分钟的单位）+直线距离 / 速度
第二次开炮与目标点的权值为——2 * T1（化成分钟的单位）+T2 + 直线距离 / 速度（以此类推
然后二分时间，再建图+增广路（要剪枝，否则会TLE
#include <algorithm>
#include <iostream>
#include <iomanip>
#include <cstring>
#include <climits>
#include <complex>
#include <fstream>
#include <cassert>
#include <cstdio>
#include <bitset>
#include <vector>
#include <deque>
#include <queue>
#include <stack>
#include <ctime>
#include <set>
#include <map>
#include <cmath>
#define CLR(x,y) memset(x,y,sizeof(x))
#define eps 1e-9
#define INF 0x3f3f3f3f
using namespace std;
typedef long long ll;
typedef long double ld;
typedef pair<ll, ll> pll;
typedef complex<ld> point;
typedef pair<int, int> pii;
typedef pair<pii, int> piii;
template<class T>
inline bool read(T &n)
{
T x = 0, tmp = 1; char c = getchar();
while((c < '0' || c > '9') && c != '-' && c != EOF) c = getchar();
if(c == EOF) return false;
if(c == '-') c = getchar(), tmp = -1;
while(c >= '0' && c <= '9') x *= 10, x += (c - '0'),c = getchar();
n = x*tmp;
return true;
}
template <class T>
inline void write(T n)
{
if(n < 0)
{
putchar('-');
n = -n;
}
int len = 0,data[20];
while(n)
{
data[len++] = n%10;
n /= 10;
}
if(!len) data[len++] = 0;
while(len--) putchar(data[len]+48);
}
//-----------------------------------
const int MAXN=60;
struct Point
{
int x,y;
}a[MAXN],b[MAXN];
struct Edge
{
int to,next;
int val;
}e[MAXN*MAXN*MAXN];
int n,m,t,t2,v;
double w[MAXN*MAXN][MAXN];
int head[MAXN*MAXN],pre[MAXN*MAXN],idx;
bool vis[MAXN*MAXN];
void add_edge(int u,int v)
{
e[idx].to=v;
e[idx].next=head[u];
head[u]=idx++;
}
void make_graph(double limit)
{
idx=0;
CLR(head,-1);
CLR(pre,-1);
for(int i=0;i<n*m;i++)
for(int j=0;j<m;j++)
if(w[i][j]<=limit)
add_edge(i,j);
}
bool dfs(int u)
{
for(int i=head[u];i!=-1;i=e[i].next)
{
int v=e[i].to;
if(!vis[v])
{
vis[v]=true;
if(pre[v]==-1 || dfs(pre[v]))
{
pre[v]=u;
return true;
}
}
}
return false;
}
double dist(int x,int y,int xx,int yy)
{
return sqrt((x-xx)*(x-xx)+(y-yy)*(y-yy))/(double)v;
}
int main()
{
// freopen("data.txt","r",stdin);
while(read(n)&&read(m)&&read(t)&&read(t2)&&read(v))
{
double t1=(double)t/60.0;
for(int i=0;i<m;i++)
read(b[i].x),read(b[i].y);
for(int i=0;i<n;i++)
read(a[i].x),read(a[i].y);
for(int i=0;i<n;i++) //第 I 个炮
{
for(int j=0;j<m;j++) //第 J 个敌人
{
int temp=i*m;
w[temp][j]=dist(a[i].x,a[i].y,b[j].x,b[j].y)+t1; //打到敌人的时间
// cout<<i<<" "<<j<<" : "<<dis<<endl;
for(int k=1;k<m;k++) //第 K 次开炮
w[temp+k][j]=w[temp+k-1][j]+t1+t2;
}
}
/*
for(int i=0;i<n*m;i++)
{
for(int j=0;j<m;j++)
cout<<w[i][j]<<" ";
cout<<endl;
}
*/
double b_low=0,b_high=(1<<31)-1,b_mid;
double ans=0;
int time=100;
while(time--)
{
b_mid=b_low+(b_high-b_low)/2;
make_graph(b_mid);
int temp=0;
for(int i=0;i<n*m&&temp!=m;i++)
{
CLR(vis,0);
if(dfs(i))
temp++;
}
if(temp==m)
{
ans=b_mid;
b_high=b_mid;
}
else
b_low=b_mid;
}
printf("%.6f\n",ans+1e-7);
}
return 0;
}
BOJ 489 小妹妹去划船
模拟题
用终点减去起点，根据差值对风向进行决策
#include <algorithm>
#include <iostream>
#include <iomanip>
#include <cstring>
#include <climits>
#include <complex>
#include <fstream>
#include <cassert>
#include <cstdio>
#include <bitset>
#include <vector>
#include <deque>
#include <queue>
#include <stack>
#include <ctime>
#include <set>
#include <map>
#include <cmath>
#define EPS 1e-10
#define INF 0x3f3f3f3f
#define ll long long
using namespace std;
template<class T>
inline bool read(T &n)
{
T x = 0, tmp = 1; char c = getchar();
while((c < '0' || c > '9') && c != '-' && c != EOF) c = getchar();
if(c == EOF) return false;
if(c == '-') c = getchar(), tmp = -1;
while(c >= '0' && c <= '9') x *= 10, x += (c - '0'),c = getchar();
n = x*tmp;
return true;
}
//-----------------------------------------------------------------------
const int MAXN=100010;
ll t,sx,sy,ex,ey;
char str[MAXN];
int main()
{
// freopen("data.txt","r",stdin);
while(read(t)&&read(sx)&&read(sy)&&read(ex)&&read(ey))
{
ll tx=ex-sx,ty=ey-sy;
int time=0;
gets(str);
for(int i=0;i<t;i++)
{
if(str[i]=='S'&&ty<0)
ty++;
else if(str[i]=='E'&&tx>0)
tx--;
else if(str[i]=='N'&&ty>0)
ty--;
else if(str[i]=='W'&&tx<0)
tx++;
time++;
if(!tx&&!ty)
break;
}
if(!tx&&!ty)
printf("%d\n",time);
else
printf("-1\n");
}
return 0;
}