T1:逗比三角形
【题目描述】
小J是一名OI退役滚粗文化课选手,他十分喜欢做题,尤其是裸题。他现在有一个二维盒子和一些二维三角形,这个盒子拥有无限的高度和L的宽度。而且他的三角形也都是一些锐角三角形或者是直角三角形。现在小J想把这些三角形放入盒子里,由于小J从txt大神犇那里学会了魔法,所以小J的三角形既可以无视盒子边界又可以重叠放置,但是必须有一条边紧贴盒子底面所在的直线。
现在小J想要最大化在盒子中的被三角形覆盖的区域的面积(即三角形间的重叠部分只算一遍),请问这个最大值应该是多少?
【输入格式】
一行一个整数T,代表数据组数。下面T部分,每部分第一行两个整数N,L分别代表三角形数量与盒子的宽度。下面N行每行三个整数ai,bi,ci表示三角形i的三条边长。
【输出格式】
T行,每行一个实数代表盒子内部被三角形覆盖的区域的面积的最大值。
T2:数三角形
Description
给定一个nxm的网格,请计算三点都在格点上的三角形共有多少个。下图为4×4的网格上的一个三角形。
注意三角形的三点不能共线。
Input
输入一行,包含两个空格分隔的正整数m和n。
Output
Sample Input
Sample Output
76
数据范围
1<=m,n<=1000
T3:树上三角形
Description
Input
Output
Sample Input
1 2 3 4 5
1 2
2 3
3 4
1 5
0 1 3
0 4 5
1 1 4
0 2 5
0 2 3
Sample Output
Y
Y
N
HINT
对于100%的数据,n,q<=100000,点权范围[1,2^31-1]
附加题:同名"数三角形"
Description
Input
Output
Sample Input
-5 0
0 2
11 2
-11 -6
11 -5
Sample Output
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#include<cmath>
#define MAXN 100005
#define EPS 1e-5
using namespace std;
double a[MAXN],b[MAXN],c[MAXN];
double S[MAXN],H[MAXN];
int n;
double LL;
bool cal(double h){
double ret=;
for(int i=;i<=n;i++){
if(H[i]>h){
ret+=(H[i]-h)*a[i]/H[i];
}
}
return (ret<=LL);
}
void solve(){
scanf("%d%lf",&n,&LL);
for(int i=;i<=n;i++){
scanf("%lf%lf%lf",&a[i],&b[i],&c[i]);
if(b[i]>c[i]){
swap(b[i],c[i]);
}
if(a[i]>b[i]){
swap(a[i],b[i]);
}
double p=(a[i]+b[i]+c[i])/2.0;
S[i]=sqrt(p*(p-a[i])*(p-b[i])*(p-c[i]));
H[i]=S[i]*/a[i];
}
double L=,R=;
double mid;
while(L+EPS<R){
mid=(L+R)/2.0;
if(cal(mid)){
R=mid;
}
else{
L=mid;
}
}
double ans=;
for(int i=;i<=n;i++){
if(H[i]>L){
ans+=(H[i]-L)*a[i]/H[i]*(H[i]-L)/2.0;
}
}
ans+=L*LL;
printf("%.6f\n",ans); }
int main()
{
freopen("sbtg10.in","r",stdin);
freopen("sbtg.out","w",stdout);
int T;
scanf("%d",&T);
for(int i=;i<=T;i++){
solve();
}
return ;
}
Code1
T2:
用组合数学的知识,先把所有的情况算出来,然后减去共线的情况即可
注意斜的用gcd算,不能用组合数了
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#define MAXN 1005005
#define ll long long
using namespace std;
ll C[MAXN][];
int g[][];
int Gcd(int x,int y){
if(g[x][y]){
return g[x][y];
}
return g[x][y]=((y==)?x:Gcd(y,x%y));
}
void make(){
for(int i=;i<MAXN;i++){
C[i][]=;
}
for(int i=;i<MAXN;i++){
for(int j=;j<=;j++){
C[i][j]=C[i-][j]+C[i-][j-];
}
}
}
int n,m;
int main()
{
make();
ll ans=;
scanf("%d%d",&n,&m);
n++;m++;
ans+=C[n*m][];
ans-=C[n][]*m;
ans-=C[m][]*n;
///!!!
for(int i=;i<=m;i++){
for(int j=;j<=n;j++){
int t=Gcd(i-,j-)-;
ans-=t*(m-i+)*(n-j+)*;
}
}
printf("%lld\n",ans);
return ;
}
Code2
T3:
这题好坑啊
设树链上的节点为a1,a2,a3,……,an
排一下序发现:
如果存在ai+ai+1>ai+2 那么就可以构成三角形了
即较小的两条边大于最大的边
然后考虑一下,如果出题目的人很变态,故意想卡你,使得这个数列尽可能多,但偏偏构成不了三角形(QAQ)
那么肯定是1 1 2 3 5 8 13 ……
斐波拉契数列,然后第50项就炸int了,由于数据又在int范围内,所以当节点数目超过50个时,再想卡你也卡不了,一定可以构成三角形
对于小于50的,直接暴力即可
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#define MAXN 100005
#define ll long long
using namespace std;
int v[MAXN];
int first[MAXN],Next[MAXN*],to[MAXN*],cnt;
int dep[MAXN],fa[MAXN];
int q[MAXN],top;
//double edge
int read(){
int ret=;
char c=getchar();
do{
ret=ret*+c-'';
c=getchar();
}while(''<=c&&c<='');
return ret;
}
void Add(int x,int y){
Next[++cnt]=first[x];first[x]=cnt;to[cnt]=y;
Next[++cnt]=first[y];first[y]=cnt;to[cnt]=x;
}
int n;
void dfs(int x){
for(int e=first[x];e;e=Next[e]){
int y=to[e];
if(y==fa[x]){
continue;
}
fa[y]=x;
dep[y]=dep[x]+;
dfs(y);
}
}
int main()
{
// freopen("sdtg2.in","r",stdin);
// freopen("my.out","w",stdout);
int T;
scanf("%d%d",&n,&T);
for(int i=;i<=n;i++){
scanf("%d",&v[i]);
}
for(int i=;i<n;i++){
int x,y;
scanf("%d%d",&x,&y);
Add(x,y);
}
dfs();
for(int i=;i<=T;i++){
int k,x,y;
scanf("%d%d%d",&k,&x,&y);
if(!k){
top=;
while(top<=){
if(dep[x]<dep[y]){
swap(x,y);
}
q[++top]=v[x];
if(x==y){
break;
}
x=fa[x];
}
if(top>){
printf("Y\n");
}
else{
int ok=;
sort(q+,q+top+);
for(int i=;i<=top-;i++){
if((ll)q[i]+q[i+]>q[i+]){
ok=;
printf("Y\n");
break;
}
}
if(!ok){
printf("N\n");
}
}
}
else{
v[x]=y;
}
}
return ;
}
Code3
附加题:
这是一道平面扫描的题目
先正难则反,找非黄金三角形
按角度排序之后,对于每个点i,与原点做一条直线,然后同在直线一侧的肯定是非黄金三角形,这样用组合数就可以计算了
我们只需要考虑左侧,因为后面的点考虑左侧时会正好可以补上,做到不重不漏
顺便说一下,向量a,b的叉积大于等于0(即a.x*b.y-a.y*b.x>=0)表示b在a的左边
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#include<cmath>
#define MAXN 100005
#define ll long long
using namespace std;
struct Point{
ll x,y;
double angle;
Point(){
x=y=;
angle=;
}
friend bool operator < (const Point &p1,const Point &p2){
return (p1.angle<p2.angle);
}
friend bool operator > (const Point &p1,const Point &p2){
return !(p1<p2);
}
friend ll operator * (const Point &p1,const Point &p2){
return (p1.x*p2.y-p2.x*p1.y);
}
}s[MAXN];
int n;
int main()
{
// freopen("data.in","r",stdin);
scanf("%d",&n);
for(int i=;i<=n;i++){
scanf("%lld%lld",&s[i].x,&s[i].y);
//atan222!!!
s[i].angle=atan2((double)s[i].y,(double)s[i].x);
}
sort(s+,s+n+);
int r=,t=;
ll ans=;
for(int i=;i<=n;i++){
while(r%n+!=i&&s[i]*s[r%n+]>=){
r++; t++;
}
ans+=(ll)t*(t-)/;
t--;
}
printf("%lld\n",(ll)n*(n-)*(n-)/-ans);
return ;
}
Code4
总结:还模版题,模版个鬼啊,这么难